






THE JOURNAL 


OF 


CHEMICAL PHYSICS 





VotumE 16, Numser 1 


January, 1948 





Dielectric Properties of Aqueous Ionic Solutions. 
Parts I and II. 


J. B. Hastep, D. M. Ritson, AND C. H. COLLie 
The Clarendon Laboratory, Oxford, England 


(Received June 20, 1947) 


The dielectric constants and loss angles of a series of concentrated aqueous ionic solutions 
have been measured at wave-lengths of 10 cm, 3 cm, and 1.25 cm. From these results the values 
of the static dielectric constant and relaxation time for these solutions have been calculated on 
the basis of the Debye formula, which appears to hold accurately. All salts show a lowering of 
the dielectric constant and a shift in the relaxation time of water. It is found that the dielectric 
constant € can be represented by a formula e=€ +26c, where €» is the dielectric constant of 
water, c is the concentration in moles per liter, and 6 has values between —7 and —15 for 


various salts in concentrations of up to 2 M. 


In Part I the measurements are described and the results discussed in relation to the 
structure of ionic solutions. In Part II the validity of the Debye-Sack saturation theory of the 
dielectric constant and the effects of the fall of dielectric constant on the electrolytic properties 


of concentrated solutions are discussed. 





PART I 
Introduction 


LARGE variety of methods have been used 

to determine the dielectric properties of 
salt solutions. The difficulty of the determination 
lies in the high value of the loss angle 6 which is 
given by the relation, tanéd = 2¢/fe where a is the 
conductivity, in e.s.u., f the frequency, and, e the 
dielectric constant. 

At the frequencies available when the early 
measurements were carried out, it was necessary 
either to work with extremely dilute solutions, 
the deviations of which from the dielectric 
properties of water are correspondingly small, or 
to work with very high values of tané, which 
greatly limits the accuracy of measurement. It is 
even difficult from a study of the literature to 
decide whether the dielectric constant of an 
ionic solution increases or decreases with in- 





creasing concentration.! It was shown by 
Hiickel? that a variation of dielectric constant 
would have a significant, if not dominating, effect 
on the properties of concentrated salt solutions. 
Further, Sack made a theoretical estimate of the 
lowering of dielectric constant that would be 
expected because of the saturation of the dielec- 
tric in the neighborhood of an ion. Both these 
studies are of great theoretical interest. It is 
therefore important to study the dielectric 
properties of ionic solutions in a region where 
tané is comparatively low, i.e., at frequencies of 
the order of 10'° cycles a second. We have pre- 
viously described methods used to determine the 
dielectric properties of water in the range \=1 
to 10 cm. The methods used depended only on 
the validity of Maxwell’s equations, and the use 

1See for summary C. H. Douglas Clark, The Fine 
Structure of Matter (Chapman and Hall, Ltd., London, 


1941), p. 232. 
2 E. Hiickel, Phys. Zeits. 26, 93 (1925). 





















































of standard attenuator technique.’ The limits to 
the accuracy of measurement were entirely tech- 
nical and did not arise from any uncertainty in 
the theory. It was found possible to work with 
reasonable accuracy with values of tané of the 
order of 1.5. 

In the present paper these methods are applied 
to ionic solutions. There are two disadvantages 
of measurements at these frequencies: 

(1) Correction must be made for anomalous 
dispersion so that measurements at a number of 
frequencies are necessary. It can be shown that 
the results obtained may be interpreted on the 
basis of a single time of relaxation, slightly dif- 
ferent from that of water. This enables accurate 
calculation of the static dielectric constant to be 
made, and further conclusions to be drawn from 
the behavior of the relaxation time. An apparent 
disadvantage may thus be turned to good 
account. 

(2) With the 2-4 percent accuracy in e, ob- 
tainable in solutions of tané=1.5, deviations 
from the static dielectric constant of water are 
only noticeable in solutions of 0.5- to 2-M con- 
centration; dilute solutions can therefore not be 
studied by these methods. 


Methods of Measurement 


Of the methods developed for the measurement 
of the dielectric properties of water and heavy 
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3R. J. Clayton, J. E. Houldin, H. R. L. Lamont, and 
W. E. Willshaw, J. Inst. Elec. Eng. 93, 97 (1946). 
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water, which are fully described by the authors,‘ 
the three most suitable have been selected for 
study of ionic solutions. These are: (1) the direct 
measurement of the absorption coefficient x in a 
rectangular wave guide; (2) the measurement of 
absorption coefficient in two cylindrical wave 
guides full of solution, the two diameters being 
selected to be just above and just below cut-off 
for the wave guide used. By this method both , 
the refractive index, and x, the absorption coef- 
ficient, are determined; (3) measurement of the 
change in Q of an Ho resonator on the intro- 
duction of an axial capillary containing liquid. 

(1) The rectangular wave-guide cell used for 
measurements of x at 3 cm and 1.25 cm is shown 
diagrammatically in Fig. 14. A pure Ho; mode is 
propagated in the liquid-filled guide, and a 
smaller guide picks up the wave after a known 
distance. The variation of response as the smaller 
guide is moved through the liquid is compared 
with the response variation of a standard piston 
attenuator. Since the attenuation is exponential, 
a plot of absolute attenuation, as determined by 
the distance moved through the liquid against 
the piston attenuator, should be a straight line, 
from the slope of which x may be calculated. The 
achievement of a straight-line plot confirms that 
a single pure mode is being propagated in the 
liquid. With a movement through 30 db of 
attenuation, with a dead space at the beginning 
of 15 db, which is possible if a superheterodyne 
method of response measurement is used, an 
accuracy of 1 percent in «x is obtained. The 
method was found to give results at 3 cm in 
agreement with those obtained by the second 
method described here. 

(2) Absorption coefficients were measured at 
both 3 cms and 10 cm by similar technique using 
the Hy, mode in filled cylindrical wave guides 
near the cut-off diameter, fed from a rectangular 
guide as shown in Fig. 1B. A loop on the end of 
concentric line is used to pick up the response. 
The propagation constants p; and 2 in the 
cylindrical guides are given by the relations: 


p2=(my—je)? =e —je’ —kyc?/w, 
po? = (n2—jke)? = €’ —je”’ — Roc? /w’, 
where ky=1.831/r; and ko=1.831/re, 7; and 72 


4C. H. Collie, J. B. Hasted and D. M. Ritson, Proc. Phys. 
Soc. 73 (1948). 
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TABLE I. Dielectric properties of NaCl solutions. Observed values of e’ and e’’, the latter corrected for conductivity, 


are compared with values calculated from the corresponding values of €, and ),. 











a 10 cm 3 cm 1.27 cm 
Tem- conduc- P €" obs. ; €" obs. P € obs. A» in 
pera- Nor- tivity €obs. corr. €obs. corr. fobs. corr. €s cm 
ture malcy  (I.C.T.) a3 +3 € calc. €' calc +1.3 +1 €calc. © calc. +1 +1 fcaic, (calc, +2 +04 
21°C 0.5 85.6 72.8 12.6 71.9 11.3 58.6 28.3 58.9 28.4 — — — — 73.9 1.71 
1.0 78.7 67.4 13.1 67.5 10.1 56.2 25.0 56.0 25.7 30.1 304 29.6 30.8 69.1 1.63 
1S 73.0 61.5 12.7 63.1 o> S25 225 S383 23.1 —- — — — 64.3 1.55 
2.0 68.6 56.2 11.8 58.0 7.9 49.6 20.7 49.4 20.6 27.8 26.0 276 26.3 59.0 1.51 
3.0 59.7 -— — — — 45.8 18.5 45.5 18.4 — — -—— — 54.0 1.47 
5.0 42.4 -—- — — -= 41.8 15.7 41.9 15.8 — —_— — — 48.4 1.39 
13° @5 54.3 --- = = -— 44.7 36.1 44.5 35.9 —— --- — -= 77.6 2.95 
1.0 49.7 — —- — — 43.0 32.8 42.7 32.4 -~ — — — 71.0 2.80 
2.0 43.3 -= -- -— — 39.4 264 39.6 26.5 — — — — 60.0 2.49 
3.0 38.0 — — — — gee 234 333 a3 oo —_ — _— 51.0 2.26 
5.0 ones — od — —_ 28.8 15.1 28.6 14.9 — — — — 38.3 2.07 
0° .66 50.1 70.7 22.9 71.2 194 45.5 36.5 45.3 36.7 — — — 79.1 2.96 
10° .66 63.2 704 19.1 72.2 145 51.9 32.1 53.4 32.4 — -— — — (EB ee & & 
20° .66 81.2 69.1 15.1 70.1 11.0 56.3 27.2 57.5 27.6 — — — 72.0 1.70 
30° .66 98.8 67.7 12.0 65.8 8.1 55.7 21.8 56.8 21.6 —- _- = 66.6 1.35 
40° .66 117.4 62.8 10.0 62.4 63 54.3 16.9 544 16.8 — — — 63.0 


1.10 








being the radii of the guides. We thus have the 
equations: 


: (c?/w?) (R?k1? — Re? ke”) 





€ 


os (k1?+ ke”), 


kK — Ko 


¢” = 2kiles’ —_ (Ry2c?/w?) +x? ]}, 


representing both dielectric constants e’ and ¢’’ 
in terms of the two absorption coefficients «; and 
xo. The latitude in the choice of r; and re is suf- 
ficient to allow the two measurements of the 
attenuation in the guide to be approximately 
independent measurements of and x, the optical 
constants in free space. 

This method is used for the real dielectric 
constant at 3 cm and for both dielectric con- 
stants at 10 cm. It was found to give results for 
water at 10 cm in agreement with those obtained 
by a similar method in rectangular guide, and 
also using an Eo19 resonator with axial capillary. 

In both methods constancy of power during 
the experiments was insured by checking back 
to the initial position of liquid cell or attenuator 
movement and noting the (unchanged) value of 
the response. Detailed descriptions of the experi- 
mental methods are given in the éarlier papers. 

(3) The Ho: resonator method for the deter- 
mination of ” at 1.25 cm is described in an earlier 
paper.‘ 

In all methods the circulation of water pumped 


from a thermostat through jackets on the ap- 
paratus provided temperatures constant to 
+0.2°C. 


Conductivity Measurements 


It is necessary to make a correction for the 
conductivity before the correct value of ¢’’ can 
be found. On solutions of which the low fre- 
quency conductivity was not available in the 
literature, rough measurements were made in a 
cell with platinum electrodes of area 0.2 sq cm, 
8 cm apart, using standard a.c. bridge technique. 
The accuracy was about 2 percent. The values 
obtained are given in the same table as the 
dielectric constant results. 


Results 


The experimental values of e’ and e’’ at 10 cm 
and 3 cm for sodium chloride solutions are given 
in Table I and those of other solutions in Table 
II. The measurements in Table II are all at 
25°+0.2°C, but those on NaCl in Table I were 
taken at various temperatures in order to study 
the effects more fully. For the same reason some 
readings at 1.25 cm were taken on the NaCl 
solutions. Corrections for conductivity are made 
according to the equation: 


€’ corsect _ €' checrved + (20/f) ’ 


where ¢ is the specific conductivity in e.s.u. and 


4 


f the frequency. For the purpose of these cor- 
rections there is assumed to be no change in this 
conductivity with frequency, an assumption 
which is discussed later. It is possible to analyse 
the corrected results in terms of a static dielectric 
constant e,, falling to an atomic dielectric con- 
stant ¢. If there is a single relaxation time the 
relations below are generally true, and not de- 
pendent on the dielectric model. 


e’ = n* — x? = (€,—€0) /(1+w?r?) +e 
= (€,—€0)/(1+[A./d ]?) +0, 
(€,—€0)wT 3 (€,—€0)[As/A ] 


ée’ =2nx= = . 
1+w?r? 1+[A./d}? 








\,=22Ccr, 


where 7 is a single time of relaxation, corre- 
sponding to a wave-length A, (‘‘Sprungwellen- 
lange’) and X, w, and c have their usual sig- 
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nificance. In the tables the experimental values 
of e’ and e” at each frequency are compared with 
the values calculated from the corresponding 
values of ¢«, and \,. These values of e, and , are 
chosen so as to give the best fit to the experi- 
mental data when the above equations are used. 

It is clear that the accuracy of measurement 
depends on the value of tané, being greatest 
where tané is smallest, that is, for solutions of 
lowest conductivity. The estimated accuracies 
given below are those for a typical case, a normal 


solution of KCI, NaCl, or RbCI. 


Accuracy in e at10cm +3 
e’ 10cm +3 
é’ 3cm +1.3 
€ 3 cm +1 
é’ 1.25 cm +1 
” tid 1.25cm +1 
Es +2 
As +0.04 


TABLE II, Dielectric properties of solutions at 3 and 10 cm, at 25°C. Observed values of e’ and e”, the latter corrected 
for conductivity, are compared with values calculated from the corresponding values of ¢, and \,. Conductivity values 
are denoted ‘‘7”’ when taken from the International Critical Tables, ‘‘E’’ when found experimentally. 








10 cm a 
i € obs. rs in 
os € obs. corr. és cm 
€ calc. +3.0 €calc. obs. +3.0 "calc, +2 +.04 





3 cm 
las €" obs. 
os : obs. corr. 
Salt fans <a alee +13 cate, €obs, +£1.0 
LiCl 0.5 81.0T 58.3 58.9 33.6 26.0 
1.0 72.9T 53.9 53.6 37.4 23.3 
1.5 67.0T 49.0 48.1 39.6 20.2 
2.0 62.0T 44.2 44.1 41.1 17.0 
RbCl 0.5 125.5E 59.5 60.8 37.0 25.0 
1.0 114.7E 56.1 57.1 45.1 23.2 
1.5 107.5E 52.5 53.4 52.8 22.0 
2.0 104.5E 49.2 49.7 59.9 20.2 
KCl 0.5 117.2T 60.0 60.8 ky 25.9 
1.0 111.97 Sie. Stok 45.2 23.6 
25 107.67 54.4 53.4 53.0 21.8 
2.0 104.27 49.0 49.7 58.4 18.4 
BaCl: 1.0 80.5T 52.3 53.4 37.7 22.2 
2.0 68.5T 44.1 43.6 43.3 16.9 
HCl 0.25 377.0T 58.3 59.3 44.2 26.0 
0.5 360.0T 55.7 56.5 63.6 26.3 
NaOH 0.25 206.0T 58.9 60.0 36.5 26.6 
0.5 196.07 57.0 55.7 43.6 24.7 
KI 0.198 126.4T 61.0 61.5 32.0 27.2 
0.396 122.0T 58.7 59.8 34.4 25.1 
NaF 0.415 83.7E 58.3 60.4 31.7 25.0 
0.830 70.2E 56.2 57.3 35.2 24.0 


LaCls 0.52 WS5E 57.3 584 33.0 25.2 
1.04 69.2E 52.6 53.3 35.9 22.0 


KF 0.242 108.9EF 61.5 61.5 32.0 27.9 
0.484 984E 58.9 59.3 34.9 25.7 
Nal 0.428 100.87 594 59.3 33.4 25.1 
0.856 944E 54.5 54.2 364 20.8 
MgCl, 0.468 84.9T 574 58.5 32.9 25.2 
0.935 744.3E 53.8 53.55 36.3 23.0 


0.5 70.0T 57.3 57.9 33.2 26.4 
1.0 59.07 55.0 55.9 36.5 25.1 
2.0 


Na2SO, 
45.57 51.2 51.7 384 20.8 





25.7 69.0 706 346 10.7 10.0 71.2 1.55 


22.55 606 63.0 51.9 8.0 8.6 64.2 1.52 
99.7 3325 35 GS 5.0 74 57.0 1.49 
17.5 51.0 1.45 
26.5 72.1 72.0 46.2 8.9 10.2 73.5 1.54 
24.4 71.0 67.0 77.5 9.4 9.2 685 1.50 
21.9 63.5 1.47 
19.7 58.5 1.43 
26.5 70.3 72.00 45.7 10.5 10.2 73.5 1.54 
24.1 643 67.0 73.9 6.8 9.8 68.5 1.50 
21.9 63.5 1.47 
19.7 58.5 1.43 
22.4 64.1 62.8 57.1 8.8 8.6 64.0 1.50 
16.8 51.0 1.42 
26.6 75.9 71.0 64.4 72 02 825 15 
25.3 69.0 1.59 
26.8 75.7 71.0 42.6 11.7 10.3 73.0 1.58 
244 665 67.0 564 4.9 9.5 68.0 1.56 
27.2 73.1 %740 25.5 104 10.7 75.0 1.56 
25.7 733 710 40.1 111 100 72.0 1.52 
26.6 71.3 72.0 28.2 74 10.3 73.0 1.55 
24.5 724 68.0 43.4 8.4 [5 @SA i132 
25.6 709 70.0 33.1 8.9 99 71.0 1.54 
22.4 66.1 63.0 51.0 7.8 8.8 64.0 1.50 
27.55 73.5 73.5 24.6 8.8 10.7 75.0 1.58 
26.1 72.1 70.5 35.9 7.3 10.1 72.0 1.56 
25.3 71.5 70.0 34.2 8.3 95 71.0 1.50 
21.5 64.0 1.42 
25.7 70.0 70.0 369 130 10.1 71.0 1.56 
22.9 64.9 63.5 49.8 8.1 8.6 64.5 1.53 
25.0 73.0 1.54 
23.4 67.0 1.49 
20.1 60.5 1.40 
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The accuracy of the calculated values of «, and 
\, depends on the correct choice of €9. This has 
been taken as 5.5, the value which has been found 
most suitable in the interpretation of our 
previous results on water. Additional weight is 
lent to this value for water by the infra-red 
measurements of Rubens® and Hawley Cart- 
wright.® It is found that , which is not much 


’ greater than 2 at 152y, rises as the wave-length 


increases. It is unlikely that a very large de- 
pression or elevation of €9 for water occurs when 
ions are added, but a variation of +1 would only 
alter the values of e, by +1 percent. 

The possibility of fitting the data to the above 
equations employing a single time of relaxation 
adds additional weight to the results; but the 
accuracy of some measurements at 10 cm is 
unsatisfactory when compared with our similar 
measurements on water, even when allowing for 
the higher loss angle. The experimental method 
was not pushed to its limit of accuracy, there 
being three principal defects: 

(1) Selection of too large a diameter of smaller 
wave guide for the absorption measurement of e’. 

(2) Unevenness in the cylindrical wave-guide 
walls in absorption measurements. 

(3) Mechanical defects of the pick-up probe 
movement. 

The first two defects arose owing to the dif- 
ficulty of drilling 5-mm diameter holes through 
distances greater than about 7 cm if commercial 
drills are used. Therefore the wave guides were 
made of extruded copper tube, which is only 
available in certain diameters, and which shows 
defects when it is strained by the soldering of 
water jackets. The movement of the pick-up 
probe was designed for measurements on water, 
through which greater distances are traveled; 
it was therefore not so accurate for the small 
movements required in measurements on ionic 
solutions. 

At 3 cm the cylindrical wave guides of diam- 
eters of a few millimeters were machined, and a 
standard micrometer used for the movement. 
There is every prospect of being able to increase 
the accuracy of these methods beyond even that 
at 3 cm, since the errors are not theoretical but 
technical in origin. 


5H. Rubens, V. D. Physik. Ges. 17, 315 (1915). 
6 C, Hawley Cartwright, J. Chem. Phys. 4, 413 (1936). 
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Fic. 2. Static dielectric constant of NaCl solution plotted 
against normality N. Temperature 21.0°C. 


Previous Measurements 


Previous measurements on ionic solutions at 
low frequencies suffer from the drawback of a 
high loss angle mentioned above. 

The small increase in the dielectric constant of 
dilute solutions observed by Wien is not in con- 
flict with the present results, since the variations 
in question are not of great magnitude, and 
might very well take place before the fall in the 
dielectric constant of concentrated solutions 
comes into play. This can be seen from Fig. 2 in 
which the static dielectric constant of NaCl 
solution is plotted against the concentration. It 
is clear, however, that the considerable fall of 
dielectric constant observed by Sack is in dis- 
agreement with the present results. Some of the 
most reliable work, that of Wyman,’ and Drake, 
Pierce and Dow* on dilute solutions at 60 cm and 
4 meters, shows that the lowering of dielectric 
constant is of such small magnitude as would be 
in keeping with our results on concentrated solu- 
tions. The only available results on concentrated 
solutions at centimetric wave-lengths are those 
of Elle’ at 4 cm using a free-wave method with 
damped waves from a Herzian oscillator, and 
those of Cooper!® on NaCl solutions at wave- 
lengths down to 7 cm using a transmission-line 
reflection method. The latter are in rough agree- 
ment with the results of the present authors, 
though not interpreted as such. (In his Fig. 4 
much better agreement would be obtained by 
employing our values of the dielectric constant.) 

7J. Wyman, Phys. Rev. 35, 623 (1930). 

8F, H. Drake, G. W. Pierce, and M. T. Dow, Phys. 
Rev. 35, 613 (1930). 


®D. Elle, Ann. d. Physik 30, 355 (1937). 
1 R, Cooper, J. Inst. Elec. Eng. 93, 69 (1946). 
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re) 
A B 
Fic. 3. Methods of hydration round an ion—A, cation; 
B, anion. 
Discussion 


(a) Conductivity 


The consistency of the results with the above 
equations, assuming the conductivity to be the 
same as the d.c. value, is in confirmation not only 
of the accuracy of the results but also of the cor- 
rectness of choice of value for the conductivity. 
It is interesting to note that the increase of con- 
ductivity of a dilute solution at a sufficiently 
high frequency predicted on the Debye-Falken- 
hagen theory" cannot be of significant size in 
these concentrated solutions. It can be seen from 
the experimental results that an increase in con- 
ductivity of more than 3 percent at 10 cm or 10 
percent at 3 cm is not likely. The Debye-Falken- 
hagen theory does not apply quantitatively to 
solutions of these concentrations, but for a 
weaker solution, N/2 NaCl, the increase would 
be of the order of 20 percent at 3 cm and 5 per- 
cent at 10 cm. It is therefore probable that the 
effect is smaller, or occurs at a higher frequency, 


11H, Falkenhagen, Electrolytes (Oxford University Press, 
London, 1934), p. 73. 
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than would be predicted for these concentrations 
by use of the theory in its present form. 

There is also the possibility of a fall of con- 
ductivity in solutions containing ions of ab- 
normal mobilities such as Ht. The abnormal 
mobility of the hydrogen is believed to be due to 
its passage through water by an exchange 
mechanism which is written as follows: 


lm tend 
| | 
H—O—H +0—H—H—0+H—O—H 
~ ~ 


It was predicted by Hiickel” that at wave-lengths 
of the order of 8 cm the conductivity of hydrogen 
ions should fall to the normal value for monova- 
lent ions. Hiickel’s treatment, like that of other 
authors, assumes a Lorentz field 


E+[(e+2)E/3] 


acting on the ions, but this should perhaps be 
replaced by an Onsager field 3E/2. With this 
modification we should not expect the fall of 
conductivity to occur until about 5 mm. 


TABLE III. Values of 5 and da,. 











51 5d\s 2 
HCl —10 0 
LiCl — 7 — 3.5 
NaCl — 5.5 —4 
KCl — 5 — 4 
RbCl — 5 — 4.5 
NaF — 6 — 4 
KF — 6.5 — 3.5 
Nal — 7.5 — 9.5 
KI — 8 — 9.5 
MgCl, —15 — 6 
BaCl. —14 —_ 8.5 
LaCl; —22 —13.5 
NaOH —10.5 — 3 
Na2SO, —11 — 9.5 

Values of 5+, 6-, 6A,*, 5A.7 
6+, d-+1 5A,t, 6Ag +2 

H* —17 + 4 
Eat —i1 — 3 
Nat — 8 —4 
Kt — 8 — 4 
Rb* — 7 — 5 
Mg** —24 — 4 
Batt —22 — 9 
Latt++ —35 —15 
: a — § — 4 
i — 3 —4 
I- — 7 —15 
OH- —13 — 2 
SO. — 7 —il1 











12 E. Hiickel, Zeits. f. Elektrochemie 34, 546 (1928). 
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TABLE IV. Minimum hydration of positive ions. 








H+ 10 
Lit 6 
Nat 4 
Kt 4 
Rbt 4 
Mg*t 14 
Batt 14 
Lat 22 








Now the maximum decrease of conductivity is 
estimated as above from the experimental results 
on HCI in Table II to be 3 percent at 10 cm and 
10 percent at 3 cm. This would indicate either 
that the replacement of the Lorentz field by the 
Onsager field in Hiickel’s treatment is justified, 
or that the exchange process occurs at a still 
greater rate. Measurements were also made of 
the absorption coefficient of 2N and 5N HCI at 
1.25 cm, and the experimental values are suf- 
ficiently large (2N, «=4.74, 5N, x=5.04) to 
indicate that at this wave-length there is no 
considerable decrease in conductivity. 


(b) Dielectric Constant 


We shall now consider qualitatively the experi- 
mentally measured depression of the dielectric 
constant of water brought about by the addition 
of ions. The most acceptable model is one of each 
ion attracting a sheath of several layers of water 
molecules around it; the dielectric properties of 
these molecules are naturally different from those 
of the molecules not in the neighborhood of an 
ion, and the effects fall off as we go further from 
the ion itself. The types of hydration believed to 
exist round a positive and a negative ion are 
shown in Fig. 3, A and B. 

In solutions of concentrations less than 2M, 
a linear relationship between the dielectric con- 
stant and the concentration is found: 


€=€,+2éc, 


where ¢, is the static dielectric constant of water, 
c is the concentration in moles per liter, and 
5=(6++6-)/2, 6+ and 6- being the (negative) 
contributions from the two ions, assumed uni- 
valent in this case. Values of 5, 5+ and 6 are 
given in Table III, it being necessary to assume 
the relative contributions of the two ions in a 
single case. We choose the case of sodium 
chloride, and 6+y,+ is taken as —8 and d-¢;- as —3 
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on the following basis; a negative ion, attracting 
water molecules as in Fig. 3B, still leaves them 
partially free to rotate, and so saturates them 
less than would a cation which attracts them as 
in Fig. 3A. Therefore, the depression of dielectric 
constant is much larger for cations than anions; 
quantitative treatment in Part II of this paper 
gives 6+= —10, 6-= —1; on the other hand, the 
differences between the saturation are probably 
not so sharply defined as this; we have therefore 
estimated that 


Stwat=—841 and sor =—3+1. 


It is possible to estimate the minimum hydra- 
tion possible round a cation from the values of 
5+, on the following basis: the model adopted in 
Part II of this paper is one in which short-range 
order dominates the orientation of the inner 
hydration sheath, effectively removing a fraction 
of the water molecules from playing any part in 
the dielectric process; outside this there is a 
region in which a small amount of saturation 
takes place, so that these molecules make a 
small contribution to the fall of dielectric con- 
stant; it is reasonable to assume that this 
accounts for about a quarter of 6+y,+, that is, —2. 
Subtracting this contribution, we have for the 
minimum number of water molecules in the inner 
sheath, N, 

(—6+—2)/e,=NM/1000, 


where M is the molecular weight of water. Values 
of this minimum hydration number WN are given 
in Table IV. 








Fic. 4. Possible hydration round a large anion. 


13D. H. Everett and C. H. Coulson, Trans. Faraday Soc. 
36, 633 (1940). 
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These results, in particular the value of N=4 
for ions such as Na* are of the correct order so 
far as the inner sheath is concerned, but it should 
be emphasized that this method can only give a 
lower limit to the hydration number; molecules 
in the first sheath are only incompletely pre- 
vented form taking part in the dielectric process, 
and therefore the hydration number will be cor- 
respondingly greater. The freedom of rotation 
of the water molecules in the sheath round an 
anion makes it impossible to treat these ions in 
the same way. 

A quantitative treatment of the dielectric 
constant depression is given in Part II of this 
paper, but the following points may be noted 
here. 

1. The cation contribution is dependent on its 
charge. 


2. Cations of small radius make a larger con- 
tribution than those of large radius. Both these 
facts are in accordance with the present ideas on 
the hydration of cations. 

3. A number of effects operate in determining 
the anion contribution, which is on the whole 
smaller than that of the cation. It is possibie 
that large anions such as I- bind the water 
molecules by both hydrogen atoms, as in Fig. 4, 
and so have a greater saturating effect than the 
small anions, which leave their hydration mole- 
cules free to rotate as in Fig. 3b. The large con- 
tribution of OH™ is explained on the basis that 
it can attract and saturate water molecules very 
easily, in a pattern similar to that of the water 
itself, without undue disturbance of the water 
lattice. 

4. As will be recalled later, values of 6+ and 
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5- are of the same order as those assumed by 
Hiickel in his theory of concentrated solutions; 
but the relative values for different ions are not 
quantitatively the same. 

It is of interest to compare our values of 6+ 
and 6~ with ionic entropy data. In a recent 
paper,'* Frank and Evans give the molal entropy 
of vaporization of ions at 25°, on changing from 
a molar solution to a gas at one-atmosphere 
pressure. This entropy change AS, is to be 
regarded as being a measure of the ordering of 
water molecules in the vicinity of an ion (com- 
pare with Latimer’) comparable to the lining 
up of the water molecules. According to Frank 
and Evans part may be due to the restriction of 
an ion into a free volume box of water molecules. 
Since the solutions are of molar concentration in 
each case, there should be a certain correlation 
between AS, and 6+, and probably a different 
correlation with 6~-. These are shown in Fig. 5. 
It is significant that 6,~-lies on the cation cor- 
relation line, indicating the rigid dipole type of 
hydration. 

The “structure-breaking entropy” of Frank 
and Evans is considered under the heading of 
relaxation-time depression. 

As can be seen from Fig. 2, the linear de- 
pendence of dielectric-constant depression on 
concentration no longer holds at concentrations 
above 2M. It is possible to regard solutions in 
this region as ionic lattices with interspersed 
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water molecules. At least we may be certain that 
the saturation effects of adjacent ions overlap. 
On this type of model the prediction of dielectric 
constant behavior is much more difficult. 

It may be seen from Table I that the dielectric- 
constant depression is greater at high tem- 
peratures than at low. This is more clearly 
shown in Fig. 6, in which the percentage de- 
pression of the dielectric constant of water by 
0.66N NaCl is plotted against temperature. This 
effect could be explained on the assumption that 
the short-range order round an ion increases 
with the temperature, as the competing order of 
the water decreases. This would imply that the 
hydration increases with the temperature, a 
problem which requires further experimental 
elucidation. 


(c) Relaxation Time 


The linear depression of the relaxation wave- 
lengths of water by the addition of ions at con- 
centrations of less than 1.5M, which is shown in 
Fig. 7, may be discussed on the same basis as the 
dielectric constant depression. The authors‘ dis- 
cuss the mechanism of re-orientation in water 
and heavy water in the light of experimental 
results. Comparison with the viscosity data of 
water and heavy water over a temperature range 
leads us to believe that re-orientation takes place 
by molecules undergoing the same process as 
that of viscous flow, or self-diffusion, that is, a 
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Fic. 7. Plot of relaxation wave-length of NaCl solution against normality N. 
4H. S. Frank and M. Evans, J. Chem. Phys. 13, 479 (1945). 


15 W. M. Latimer, Chem. Rev. 18, 349 (1935). 
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Fic. 8. Speculative correlation of fluidity elevation with 
sum of short-range and long-range effects. Anions, repre- 
sented as blackened circles, lie on shorter line. Errors of 
+6 in 6+—6d,* not marked in. 


replacement of empty lattice spaces by molecules 
which themselves leave further spaces for re- 
placement. Thus re-orientation takes place at 
the moving boundaries of the broken pieces of 
“lattice.” 

Most ions when added to water rupture the 
lattice, creating more boundaries at which re- 
orientation can take place; they have a so-called 
“structure-breaking”’ effect, of which the de- 
pression of relaxation time or wave-length is a 
reasonable index. 

Using a similar notation to that of dielectric- 
constant depression, we have at concentration 
lower than 1.5M; \,=Asa+2d,c, where Asw is 
the relaxation wave-length of water, and 4,, 
=(6A,++6A,-)/2, 5A,+ and 6A, being the 
(negative) contributions from the two ions, both 
taken as univalent in this case. In Table III are 
given values of 6\,+, 5,-, and 6A,, it being 
assumed that in a solution of ions of equal radius, 
such as KF, there are equal contributions from 
K and F (compare with Bingham"). The error 
of +2 in 6d, is such that these values can only 
be regarded as semiquantitative. 

It is apparent that 6\,+ is dependent on the 
charge of the cation, and is greater for large ions 
than for small, much-hydrated ions of the type 
of Lit, Mgt+, and H+. The concept of a large 
ion breaking the structure more than a small 


16 E, C. Bingham, J. Phys. Chem. 45, 885 (1941). 
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ion, because it is difficult to fit into the structure 
without splintering, is reasonable. Moreover 
these qualitative differences are in agreement 
with the generally accepted views on the inter- 
pretation of x-ray diffraction, Raman and 
infra-red absorption spectra, and on the de- 
pression of the temperature of maximum density. 

It is usual to correlate relaxation-time data 
with macroscopic viscosity data, as was done by 
the authors in the case of water and heavy water. 
In dilute solution the increase of viscosity is satis- 
factorily explained on the theory of Falkenhagen 
and Dole,!? which employs the concept of the 
ionic atmosphere, but at high concentrations 
there is a negative viscosity effect which requires 
a different explanation. Finkelstein'* developed a 
theory of the re-orientation of solvent molecules, 
giving a positive viscosity effect varying with the 
concentration; but we must be able to explain 
the effects of various ions in concentrated solu- 
tion in both elevating and depressing the vis- 
cosity. In this connection the ‘“‘ionic elevation of 
fluidity” of Bingham" is of value. Let us first 
examine the relaxation-time data purely from 
the point of view of “‘structure-breaking”’ effect. 
Frank and Evans" tabulate values of a quantity 
AS,: which can be regarded as the entropy change 
due to the disturbing effect of the ion and hy- 
dration sheath on the surrounding water; but it 
contains other factors also. We should expect this 
to be correlated with 6\,+ and with Bingham’s 
elevation of fluidity. The values of 6A,+ are not 
sufficiently accurate to correlate with AS,,, but 
for monovalent ions there is no violent disagree- 
ment. For polyvalent ions, however, AS,; is not 
calculated, and large negative values of fluidity 
elevation appear. Now polyvalent ions are shown 
on the basis of relaxation time to have a strong 
structure-breaking effect and they should there- 
fore give large positive elevations of fluidity. 
There is, therefore, a second factor operating in 
the cases of viscosity or fluid flow. This is the 
short-range order effect, by which water mole- 
cules in the inner hydration sheath, by reason 
of their strong orientation, are made more vis- 
cous and depress the fluidity. This effect, 
analogous to the depression of dielectric con- 
stant, will be greater for polyvalent ions, and 


17 See H. Falkenhagen, reference 11, p. 245 for summary. 
18 B, N. Finkelstein, Phys. Zeits. 31, 130, 165 (1930). 
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will be correlated to 6+. A speculative treatment 
involving both effects, that is, correlating 
fluidity elevation, with values of 6+—6,,*, is 
shown in Fig. 8, but the present level of accuracy 
is such that further measurements of relaxation 
wave-length are necessary before this picture of 
viscosity can be regarded as fully justified by 
relaxation-time measurements. 

The variation of depression of relaxation 
wave-length with temperature is shown in Fig. 6, 
in which the percentage depression by 0.66M 
NaCl is plotted against temperature. The effect 
of the ions is greater at low temperatures than at 
high, as can be seen from the following picture; 
the low values of the relaxation wave-length of 
water at low temperatures indicate that the 
boundaries at which re-orientation takes place 
are much more rare; yet the effect of the ions 
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in introducing new boundaries is of the same 
order as at higher temperatures, since the size 
and number of the hydration sheaths is of the 
same order; therefore, the introduction of the 
ions into the water will have a proportionately 
greater effect on the frequency of occurrence of 
the boundary, and so on the relaxation wave- 
length. 


Conclusion 


We therefore find that it is possible to give a 
reasonable qualitative explanation of our results 
on depression of dielectric constant and relaxa- 
tion time on the basis of the model outlined 
above. In Part II of this paper we shall consider 
the depression of dielectric constant from the 
standpoint of dielectric theory, and the bearing 
of these results on electrolyte theory. 


PART II 


D. M. Ritson AND J. B. HAsTED 


In Part II of this paper the theory of dielectric saturation in the neighborhood of an ion is 
shown to give values for the depression of dielectric constant of concentrated ionic solutions 
of the same order as those found experimentally in Part I. Modern theories explaining the de- 
viations of concentrated solutions from the Debye-Hiickel theory are considered, and semi- 
quantitative estimations made of the various contributing factors—viz.: Dielectric constant 
variation, Bjerrum-type association, co-volume correction, and structure-breaking effect. 


Theoretical Treatment of the Dielectric 
Constant Depression 


It was predicted by Sack! that the dielectric 
constant of an aqueous ionic solution would be 
considerably lower than that of water. On the 
basis of the Debye theory of dielectrics he cal- 
culated a depression which was about a hundred 
times greater than the experimental values re- 
corded in Table III, Part I of this paper. 

In general we should expect that any theory 
based on the assumption that an ion is sur- 
rounded by a continuous dielectric medium of 
macroscopic dielectric constant, should give a 
lower limit to the actual dielectric saturation in 
the neighborhood of an ion; short-range order 
round the ion, which is ignored when we treat the 
surrounding medium as continuous, must further 
increase the saturation. 


1H. Sack, Phys. Zeits. 27, 206 (1926); ibid. 28, 199 
(1927). 


The apparent discrepancy between calculation 
and experiment arises from the use of a Debye 
model for the dielectric constant, which assumes 
an improbably small value for the dipole moment 
in water. We have attempted in this paper to 
bring Sack’s original conception of the electric 
field of an ion effectively saturating the sur- 
rounding dielectric into line with more recent 
views on the structure of water. 

There are two fundamental difficulties en- 
countered in any such treatment. There is still 
considerable uncertainty as to the correct model 
for the dielectric constant of water, and there is 
no method of treating the short-range order 
round the ion. One is therefore limited to ob- 
taining a qualitative estimate of the effect, and 
showing that this estimate is in agreement with 
experiment. 

We consider first a point charge placed in a 
continuous dielectric medium. For this dielectric 
medium we take two models. 
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Model 1 is a dielectric obeying Onsager’s 
equation exactly. At 25°C we can approximate 
to water by taking values for the dipole moment 
wo Of 1.88 Debye units and refractive index 
n®? = 4.13, leading to a dielectric constant of 78.7. 
‘There is one grave drawback to this model, 
namely, that it gives an effective moment (the 
dipole moment plus the dipole moment induced 
by the reaction field) of 3.77 Debye units, which 
is rather improbable. 

We obtain for the dielectric constant « 





An (n?+2)(2e+1) fe(m?+2) po? 
E=— 1 —E| 

V = 3(2e-+n?) (2e+n?) kT 
3e(n?—1)E 


2e+n? 








where L is the Langevin function, and all the 
other symbols have their usual significance. 


We can make the approximation 


An (n?+2) /(m?+2) po? 3 
E=— 1( “z) +-(n?—1)E, 
V3 2 kT 2 





€ 


and on substituting the values given above we 
find at 25°C that 


eH =5.69E+1.58-10°L(1.39-10£). 


In addition to this treatment, which is open 
to the criticism given above, we have taken a 
second model. This is an empirical modification 
of Kirkwood’s? formula for small fields, so that 
it applies to large fields. Kirkwood’s formula 
gives 


cE = (49/V)LaE+(ua/3kT)E |, 
where zu is the effective dipole moment, and @ is 


the moment plus the moment induced in the 
surroundings. 
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Fic. 1. Plot of dielectric constants ¢,, ¢g against distance 7 in angstroms from an ion, 


2 J. Kirkwood, J. Chem. Phys. 4, 592 (1936). 
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We have made the empirical modification 
eE = (42/V)aE+uL(u/kT)E. 


Putting u=2.45 Debye units, which is Ver- 
wey’s® estimate of its magnitude, we obtain 


cE =4E+1.03-10°L(2.18-10-4E). 


Although this second formula is certainly not 
entirely correct, our main interest is to show that 
our value for the saturation effect is not markedly 
dependent on the model we adopt. 

Now consider the point charge as the center 
of coordinates 7; and 6. We know that the field 
E, round the ion must obey the equation 


div(E+4P) =0, 
or 
E+4nP=e/r’. (1) 
We shall write this equation in the form 
AE,+BL(CE,) =e?/r’, (1a) 
where for, model 1, 
A=5.69, B=1.58X10®, C=2.18X10-, 
and for model OF 
A=4, B=1.03X10®, C=2.1810~. 
We now define two dielectric constants 
€-=OP/0E,, «=OP/dEs. 


The significance of these quantities lies in the 

fact that they represent the effective dielectric 

constants when a small additional field is applied. 
From Eq. 1(a) we see that 


€,=A +BC(dL(CE,)/d(CE,)], 3(a) 
and 


These are the effective dielectric constants when 
a small additional field is applied either along or 
perpendicularly to the radius vector. 

We have plotted in Fig. 1 values of €, and e 
against those calculated from Eqs. 3(a), 3(b), 
and (1). 

When r>4A, the values of e€ do not differ 
appreciably from e, and can be expressed ac- 
curately by the approximations 


€r = €s|.1 — (3b/r*) ], (4a) 
€9 =e, 1—(b/r*) ], (4b) 
7 E. J. W. Verwey, Rec. Trav. Chim, 60,8127 (1941). 


where in model 1, }=57.8X10-*, and in 2, 
b=163X10-*. These results are obtained by 
expanding L(x) as (x/3)—(x*/45)+---, and 
neglecting terms above the second. 

Returning now to the original problem of the 
effect of an ion of finite size on the surrounding 
dielectric, it is apparent that in the region below 
4A, which is of the order of the radius of the 
first hydration sheath, the representation of the 
dielectric as continuous will not be a valid approx- 
imation. Therefore, for simplicity, we have made 
the following assumptions: 

(a) That we can represent a positive ion as 
surrounded by a completely saturated first shell, 
and beyond that by continuous dielectric of 


‘dielectric constant given by (4a) and (4b). We 


assume that a saturated shell has a radius r,+2r,, 
where r is the radius of the ion, 7, that of the 
water. 

(b) That the negative ion can be’ represented 
by a completely unsaturated first shell of dielec- 
tric constant €,, surrounded outside a radius of 
r_+2r., by continuous dielectric. 

We shall now consider the application of a 
small field to the region round a point charge. 
We can represent the applied field E by a poten- 
tial y. We shall only briefly indicate the argument 
which can be found in full in Debye’s Polar 
Molecules, Chapter VI. As r-+* we know that 
yY must be of the form 


v= —([E,+(u/r*)] cosd. (5) 
We assume that 


¥= —LEf(r)+(u/r*)g(r) ] cosé 
= E-R(r) cos, 
for 
r>ra+2r,. 


Now it is a fundamental equation of electro- 
magnetic theory that 


dive-grady =0. 


This gives, on rearrangement, 





dj dR(r) 
|e" |- 2eR(r) =0. (6) 
dr dr 
The solution of this equation is 
f(r) = —r[1—(2b/r')], (6a) 


g(r) = (1/r*) [1+ (6/r')], (6b) 
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if 
r<ra+2r,. 

The solution must be of the form ~=Ar cos@, 
where A is a constant, as the dielectric constant 
is now independent of r and 0. 

Putting in the boundary condition that y is 
everywhere continuous and that e,(dy/dr) is also 
continuous, we obtain 


p=(e9—1)/(2e+1)(r++2r.)*. 


The total polarization P, outside r=r++2r, is 


given by 
rank (cet 


where dr is any volume element. Therefore, 


dy 
P=af fr(S+E sin@-drd@. (7) 
T+ +2ry 


We require to obtain the decrease in polariza- 
tion AP of the medium caused by the presence 
of the central ion. If we write P» for the polariza- 
tion of the medium if no ion were present, P, 
for the polarization in the first sheath round the 
* jon, we obtain 


AP,=P,)— 


Ss 


—Pr,—Ps,. (8) 


The dielectric constant e, for a normal solution 
is, therefore, ¢,=€.+25c, where §6=(6++6-/2) 
and 

5+ = NAP, -4x-10-%, (9) 


Substituting the values we have obtained for 
€r, €, and y, and ignoring terms of the order of 
(b/r*)?, we find 


AP=[—Db/4(r+2r..)4|—P:. (10) 


It also follows from our initial assumptions that 
for negative ions P,=0, for positive ions 
P,=r'*e,/3. This leads directly to the following 
values for 6+ and 6~ for ions of radii r, =r_=1A, 
and r. equal to 1.36A: 


Model 1, 6+= —11.5, 6-= —1.5, and b= —6.5. 
Model 2, 6t+=—13, 6-=-—3, and 6=—8. 


The calculated values of 6 are thus of the right 
order of magnitude. The main contribution to 
the dielectric lowering comes almost entirely 
from the first sheath. The large increases for 


Lit, H+, and OH™- ions are undoubtedly due to 
the extent of short-range order round the ions. 
For divalent and trivalent ions the above calcu- 
lation is unsatisfactory, as the short-range order 
probably extends further than the first shell. 


Theory of Concentrated Solutions of 
Strong Electrolytes 


In 1926 Hiickel* suggested that the deviations 
from the Debye-Hiickel theory of strong elec- 
trolytes at moderate concentrations could be 
explained by supposing that the dielectric con- 
stant was a function of the concentration. He 
suggested the law for the dielectric constant D 
of a solution of molecular concentration c, 
D=D,—28, where Dy is the dielectric constant 
of the solvent; he showed that this led to a 
formula for the activity of an electrolyte of the 
following type: 

logiof=D+Be, (11) 


where D is the usual Debye-Hiickel term of the 
type [0.43432,2_e?x/2DokT(1+«a) ], where 


[= y_z_?) } 
k= , 
DokT 





(11a) 


n is the number of molecules of salt per cc, vy 
and v_ the numbers of positive and negative ions 
of valencies z; and z_, corresponding to one 
molecule of salt, e the electronic charge in e.s.u., 
k the Boltzmann constant, and a the ionic radius 
measured in angstrom units. 

Hiickel’s term Bc, due to the variation of the 
dielectric constant, is given by the relation: 





Be= 


0.4343e2 sv,6+ y_d- 
( ). (11b) 


nDP?kT a~ 

Now the type of relation logif=D+Be is a 
well-known explanation of the variation of 
activity with concentration, found experimen- 
tally. Values of Bexperimenta for a large number 
of salts are given by Harned and Owen.*® Using 
§ and a as adjustable parameters, Hiickel ob- 
tained remarkedly good agreement with these 
values Of Bexperimentar: We were interested to find 


4E. Hiickel, Phys. Zeits. 26, 93 (1925). 

5 Harned and Owen, The Physical Chemistry of Electro- 
lytic Solutions (Reinhold Publishing Corporation, New 
York, 1943). 
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that our experimental values of 6+ and 6-, given 
in Part I of this paper, were of the same order of 
magnitude and showed the same qualitative 
trends as those used by Hiickel, although quan- 
titatively the agreement is not good. 

Later authors rejected the view that the devi- 
ations were due to changing dielectric constant, 
and attributed them to other causes. It is likely 
that a number of different factors contribute, 
each one involving a term Bc; we shall consider 
each in turn and attempt to assess semiquanti- 
tatively its contribution to B. They are as 
follows: 

1. Van der Waals co-volume correction. 

2. Correction due to the breaking up of the 
water structure by the ions. 

3. Correction due to the type of ionic associ- 
ation envisaged by Bjerrum. 

4, Varying dielectric constant correction, out- 
lined above. 


Co-Volume Effect 


Since the distance of closest approach of two 
ions is fairly large, it would seem likely that a 
van der Waals co-volume correction must be 
made to the activity. Van Rysselberghe and 
Eisenberg*® have calculated values of B on this 
basis, using Ursell’s treatment. We have inves- 
tigated whether it is legitimate to apply such a 
correction by assuming a highly simplified lattice 
model. 

If we consider the distribution of ~ molecules 
on WN lattice points, such that around each 
molecule the neighboring positions must be 
unoccupied, we find that the number of dis- 
tinguishable configurations Q is given by 


‘iol | N! 
(N—2n)!N! (N—n)!n! 
for N>>n. (12) 





This result is obtained by using a very general 
method of enumerating lattice configurations, 
which we shall describe elsewhere. The equation 
for Q follows from the fact that we have x axes 
along which each center must be separated by at 
least two lattice points. The probability that any 
random configuration satisfies this condition 


6 P. van Rysselberghe and S. Eisenberg, J. Am. Chem. 
Soc. 61, 3030 (1939). 


along any axis is given by the expression in the 
brackets. The power x gives the simultaneous 
probability that we can satisfy the condition 
along the axes simultaneously. N!/(N—n)!n! 
represents the total number of random dis- 
tinguishable configurations. 

We now consider a model approximating to a 
number of solvated molecules B in a solvent A. 
A two-dimensional model is shown diagrammati- 
cally in Fig. 2. The B molecules and their solva- 
tion sheath occupy one unit cell of the lattice. 
This type of model is chosen as being similar to 
water, where the, coordinated sheath of water 
molecules fits into the general lattice structure. 
For simplicity we| consider a cubic lattice. 

Let there be NV, solvent molecules A, and Ng 
solute molecules B. 8Nz of the solvent molecules 
A are held by the Nz solute molecules. There are 
13 axes along which the B’s must be separated 
by at least one cell. Four axes are shown in the 
two-dimensional figure. The number of positions 
that it is possible for B molecules to occupy is Na. 
Making the assumption that the solution is 
“regular,” i.e., that the potential and molecular 
energies are not altered for the different con- 
figurations, we can write down the configura- 
tional partition function with the aid of Eq. (12) as 





Sy 





(Na—2Nz)!Na! 
Na 
. » 3s) 
(Na—Nz)!Nez! 
The ‘‘partial potentials’ are given by 
wa=pw0a—kT(d InQ/dNa), (14a) 
p= u0p—kT(d InQ/dNz). (14b) 
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If we consider the case for N4 >Nz so that we 
can ignore term in (Ng/Na)* we find that 


ba poa+kT In(Na/Na+Na), (15a) 


Mp=pop+kT In(Ne/Nat+Na) 
X[1+28(Ne/Na)]. (15b) 


If we had assumed van der Waal’s law we should 
have written the equation as 


wa=poat+kT In(Na/Na+Ns), (16a) 


bp=popt+kT In(Ng/Niat+Nzs) 
<[1+32(Ne/Na)]. (16b) 


We see, therefore, that in the particular case we 
have considered above, the expressions reduce 
to a form very near to that of van der Waals law; 
it follows that the treatment of ions in solution 
may be made along similar lines to that of 
molecules in a gas, and it is likely that the appli- 
cation of Ursell’s method may be made with 
success. Harned and Owen’s® revision of van 
Rysselberghe and Eisenberg’s® treatment gives 
a value of B=2.2063 X10-* Xa‘ for a single 1:1 
electrolyte, or 32.2063 X10-* Xa for a single 
2:2 electrolyte, where a represents the distance 
of closest approach of two ions measured in 
angstrom units. For consistency we have used 
the empirical values of a which are obtained by 
fitting the experimental activity coefficients to 
Eq. (11). These values are given by Harned and 
Owen? and are greater than the crystallographic 
radii. 


Structure-Breaking Effect 


The interactions of the ion with the water 
lattice is also an important contributing factor; 
it has been brought out clearly in a recent paper 
of Frank and Evans.’ The authors point out that 
for ‘‘structure-breaking” salts the water lattice 
becomes increasingly random with increasing 
concentration. Accordingly, dilution should be 
accompanied by a decrease of entropy. The 
change of entropy of dilution with concentration 
is a measure of the structure-breaking effect of 
an ion, and values obtained by this method show 
correlation with those obtained from Raman 
and infra-red absorption spectra, viscosity data, 
etc. There is, in addition, some sort of correlation 


7H. S. Frank and M. Evans, J. Chem. Phys. 13, 479 
(1945). 
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with the shift of relaxation time of the water by 
the ion, experimental values of which were given 
in Part I of the present paper. When the relaxa- 
tion-time measurements are considered in con- 
junction with the previous study by the authors 
of the re-orientation process in water and heavy 
water, it appears likely that relaxation-time shift 
can be regarded as a fairly reliable measure of 
structure-breaking effect. As was demonstrated 
by Falkenhagen,® it is impossible to obtain the 
correct ‘heats of dilution from Hiickel’s theory, 
and it would seem equally difficult for modern 
theories to account fully for the activity coef- 
ficient deviations without taking lattice effects 
into account. It is, however, the free energy and 
not the entropy change with which we are con- 
cerned, and this is, of course, smaller. We found 
in the first part of this paper that the effect 
was much greater at low temperatures than at 
high, and this will lead us to predict that activity 
deviations are likely to be much greater near the 
freezing point than at high temperatures. 

However, it is not yet possible to calculate 
even approximately the contribution made by 
this effect to the value of B, and we have, there- 
fore, made an estimate simply by subtracting the 
other contributions from the experimental value. 

It will be seen from Table III where the con- 
tributions to B are tabulated, that the structure- 
breaking contributions so obtained are in ac- 
cordance with the usual comparisons between 
various ions, i.e., 


Bat+>Rbt++ >Kt>Nat>Ht, 


and are, moreover, fairly small. 


Bjerrum-Type of Ionic Association 


Bjerrum® attempted to account for anomalies 
at high concentrations by assuming that ionic 
association occurs. He calculates the electro- 
static and other forces acting between two ions, 
and equates the resulting potential energy with 
mean kinetic energy at various inter-ionic dis- 
tances. Within a certain distance the ions will 
“associate.” The association constant (1—a) 
depends on concentration in a complicated, but 


8H. Falkenhagen, Electrolytes (Clarendon Press, Oxford, 
England, 1934), p. 282. 

9N. Bjerrum, D. Kgl. Dansk. Vid. Selsk. Medd. 7, No. 
9 (1926); Ergeb. der_exakt. Naturwiss. 6, 125 (1926). 
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approximately linear, manner, but may be shown 
to be of the order of 0.14 in molar solution of 1:1 


-electrolyte of ionic radii sum equal to 2.82A. For 


higher concentrations, smaller ions, ions of 


_ higher valency, and solvents of low dielectric 


constant the association will be higher. 

It is unfortunately impossible to accept the 
calculated values of (1—a) on Bjerrum’s theory, 
or to calculate back from the heats of dilution, 
since a number of approximations are made 
which are not permissible in more concentrated 
solutions; moreover, at distances of less than 4A 
from an ion the dieclectric constant is a rapidly 
varying quantity, as is shown above. However, 
it appears very likely that some sort of associa- 
tion of this type occurs in concentrated solutions 
because there certainly are strong forces acting 
between ions at very short distances. 

The only method of estimating Bjerrum’s 
association that remains is by treatment of the 
conductivity data. Let us consider in detail the 
simplest case, a 1:1 strong electrolyte at 25°C, 
in which there can be little question of incom- 
plete dissociation in the usual sense of the term, 
but for which accurate conductivity data over a 





wide concentration range exist. Such a case is 
potassium bromide, for which there are the 
measurements of Jones and Bickford.'!® Between 
concentrations of 1.5—6N many such salts, 
including the present case, show a linear de- 
pendence of equivalent conductivity on concen- 
tration. This is illustrated for a number of cases 
in Fig. 3, for which the data have been taken 
from the International Critical Tables. It is true 
that many salts show a non-linear dependence in 
such concentrated solution, for instance, salts of 
lithium. But, nevertheless, the dependence in a 
large number of the simplest cases is linear, and 
this fact must be explained. Now it is with solu- 
tions of these concentrations that we are con- 
cerned when investigating Bjerrum association, 
because the effect is more significant in this 
range. It is well known that for the conductivity 
of dilute solutions, many equations with satis- 
factory theoretical basis have been proposed 
(Onsager), but all explain the dependence on the 
square root of the concentration which is found 

10G, Jones and C. F. Bickford, J. Am. Chem. Soc. 56, 


602 (1934). 
L. Onsager, Phys. Zeits. 28, 277 (1927). 
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TABLE I. Approximate percentage Bjerrum association 
100(1 —a) calculated from conductivities of concentrated 
solutions. 








Approximate 





percentage 
association 
Electro- estimated for Temperature 
lyte normal solution 7 Data 
HCl 11 18 I.C.T. 
NaCl 10 0 LCT. 
NaCl 10 18 fACT. 
KCl 5 18 iB oy 
Nal 9 25 iC oe 
KBr 5 25 Jones and 
Bickford * 
KI 5 18 LR tg ip 
RbCl 4 18 iS aay 
CsCl 4 18 A.T. 
BaCl, 12 25 Jones and 
Dole ** 








* See reference 10. 
** See reference 12. 


in dilute solution; modifications are made to 
Onsager’s equation to explain deviations from 
this. It is very unlikely that such modifications 
will transform the dependence into a linear one 
in concentrated solution, by the addition of a 
term which counter-balances the non-linear 
dependence, nor is it likely that a viscosity cor- 
rection will transform the relationship into a 
linear one; in other words, a single equation of 
the accepted type will not hold over the entire 
concentration range unless the Onsager term 
reduces to a constant value and a linear term, 
— K,, is introduced at high concentrations only. 
The empirical equation of Jones and Dole” is an 
example of this type of behavior: 


A=Aw—Av/c/(1+Ga/c) —Ke, (17) 


but in order to fit the dilute solution data, the 
constants of this equation are so chosen that it 
does not reduce to a linear equation until very 
high concentrations (~10N) are reached. Al- 
though this equation is capable of explaining the 
experimental data over the widest concentration 
range, and is of a theoretically satisfactory form, 
the parameters are chosen entirely empirically, 
and a single equation of this type would not be 
capable of explaining both the dilute solution 
data and the linear dependence in concentrated 
solution. In particular, the values of K are such 
that it is difficult to find a theoretical explanation 
of them. It is necessary to find a theoretical 


122G. Jones and M. Dole, J. Am. Chem. Soc. 52, 2245 
(1930). 
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equation which holds in dilute solution and 
reduces to a constant term by the time a con- 
centration of about 2N is reached; at greater 
concentrations, a new effect comes into play to 
explain the linear dependence. We shall assume 
that this effect is the Bjerrum association. Its 
linear dependence on concentration cannot be 
predicted theoretically, in fact, as we have seen 
above, the dependence on concentration is 
believed to be complicated, but only approxi- 
mately linear. But in view of the criticisms of 
the theory given above, the original treatment is 
not unsatisfactory. 

We have thus two equations involving the con- 
centration on which the conducivity depends: 

(1) An equation which we will call C, tending 
to a constant conductivity value at concentra- 
tions greater than 1.5M. The form of this equa- 
tion will be discussed later. 

(2) A simple equation D, holding over the 
linear range; A=aA’ where 1—a is the degree of 
association and A’ is the constant value to which 
equation C has tended. We can estimate the 
association approximately from the slope of the 
conductivity-concentration curve. The accuracy 
is limited by the following facts: 

(1) It is not possible to obtain a good value 
for A’; we have used what seems the most reason- 
able value, namely, that obtained by producing 
the linear section back to zero concentration. 

(2) It is doubtful whether a viscosity correc- 
tion should be applied and if so what form it 
should take. There seems to be insufficient jus- 
tification for simply multiplying the conductivity 
by the macroscopic viscosity, as, for instance, 
Jones and Dole have done. We have omitted the 
viscosity correction entirely. 

The approximate values of percentage associ- 
ation in normal solution estimated in this way 
for a number of electrolytes are given in Table I, 
to the nearest whole number. The qualitative 
dependence on ionic radius and valency are 
according to theory. 

It remains to discuss the form of equation C. 
It is possible to explain the data below 0.02N on 
a theoretical basis, and at the same time obtain 
a relation which is capable of tending to a con- 
stant value. The equation 


A=A,—Av/c/(1+G/c) 
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which was first suggested by Lattey'® as an 
empirical relationship, and was used with the 
addition of a Ke term by Jones and Dole, has 
actually a theoretical basis provided the con- 
stants are chosen correctly. Kaneko“ suggested 
that the introduction of a factor of 1+xa in the 
denominator would correct the Onsager theory 
for the finite size of the ion. Although this sug- 
gestion has not been given much consideration, 
and although the treatment has been criticized, 
it is apparent that this factor corresponds to the 
1+.xa which appears in the denominator of the 
Debye-Hiickel expression for the logarithm of 
the activity coefficient of a concentrated solu- 
tion; it seems perhaps a plausible correction. 
For dilute solutions we have the equation (C’): 


A=A,—Ov/c/(1+ xa), (C’) 


where O,v/c is the theoretical Onsager term, a 
modification of the original Debye-Hiickel ex- 
pression. 
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This type of treatment is supported by our 
experimental results on Debye-Falkenhagen con- 
ductivity increase at high frequency. We found 
the increase in concentrated solutions to be much 
smaller than would be expected on the original 
theory. If the Debye-Hiickel conductivity term 
is diminished by the factor 1+ xa, this becomes 
intelligible, since it is only this term which is 
annulled at high frequencies. An independent 
term reducing the equivalent conducivity of con- 
centrated solutions below its value at infinite 
dilution, the type of term we introduce as being 
due to Bjerrum association, would not be 
affected : so that the high frequency conductivity 
would not show such a large increase above the 
d.c. conductivity. 

In the above equation for dilute solutions 
there are no empirical factors, but the equation 
will not hold except in dilute solutions, although 
its range is greater than that of the Onsager 
equation. The exact position for potassium 


TABLE II. Experimental data of Jones and Bickford for KBr equivalent conductivity at 25°C 
compared with various equations. 





























Ao(Onsager) 
calculated 
back from AJonesand AAJones and APresent AAPresent 
Cc Aobs. obs. AShedlovsky AAShedlovsky Dole Dole equation equation 
0.00025 150.16 151.65 150.16 0.00 150.24 —0.08 150.16 0.00 
.00036 149.87 151.66 149.88 — 01 149.93 — .06 149.89 ,— 02 
.0005 149.55 151.66 149.57 — .02 149.61 — .06 149.57 — .02 
.00075 149.12 151.70 149.12 .00 149.13 — 01 149.12 .00 
.001 148.78 151.76 148.75 .03 148.74 .04 148.75 .03 
.0016 148.02 151.79 148.02 .00 147.98 .04 148.02 .00 
.002 147.64 151.86 147.62 .02 147.56 .08 147.61 .03 
.005 145.47 152.14 145.46 01 145.36 ll 145.44 .03 
01 143.15 143.19 — .04 143.11 .04 143.12 .03 
.02 140.26 140.25 01 140.25 01 140.07 _-15 
.05 135.44 135.36 .08 135.51 — .07 134.74 .70 
1 131.19 131.19 .00 131.27 — .08 
2 126.59 127.61 —1.02 126.62 — .03 
a 120.35 127.33 —6.98 120.23 AZ 
1.0 115.46 115.50 — .04 
2.0 109.37 110.98 —1.61 
3.0 103.55 108.39 —4,.84 
3.75 98.70 106.98 — 8.26 
Onsager Ao=(A+59.790+/c)/(1 —0.2274/c) 
Shedlovsky A= 151.63 —94.27+/c+ 100.96c(1 —0.2274) /c 


Jones and Dole 


A= 151.84 —103.77+/c/(1+1.8953+/c) —0.492c 


Present equation A=151.63 —94.274/c/(1+1.084/c) 
Horizontal lines indicate limits beyond which equations do not fit the observed data. 








%R. T. Lattey, Phil. Mag. 4, 831 (1927). 


4S. Kaneko, J. Chem. Soc. Japan 56, 793, 1320 (1935); ibid. 58, 985 (1937). 
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Fic. 4. Qualitative variation of equivalent conductivity 
with concentration according to various equations. E, 
experimental; C, C’, D, equations as in text; O, Onsager; 
S, Shedlovsky; JD, Jones and Dole. 


bromide at 25° is given in Table IJ,* and the 
type of behavior of the Onsager, Shedlovsky,'® 
and present equations is shown qualitatively 
in Fig. 4. Our values of a are from the ionic 
radius in the crystal. 

There is no doubt that the equation of Jones 
and Dole fits the data over the widest concen- 
tration range, but it fails at high concentration 
and is moreover entirely empirical. The viscosity 
correction which the authors apply widens the 
range, but will not reduce the equation to a 
linear variation; moreover, it is of doubtful 
theoretical validity. The Onsager equation 
deviates steadily at concentrations greater than 
0.001N, but when modified, as suggested here, 
will hold to 0.02N, whereas the Shedlovsky equa- 
tion with one empirical factor deviates from 
0.1N. We have tried our modification of the 
Onsager equation on Shedlovsky’s data for HCI, 
LiCl, NaCl, KCl,!* for which the agreement is 
satisfactory up to about 0.01N. 

The form of the equation is therefore promising, 
but it cannot be regarded as even approximately 
correct in concentrated solution, since the as- 
sumptions underlying the Onsager term do not 
apply. We cannot therefore offer any detailed 
explanation of how the equation of type C is 
made up, and must be content simply to offer 
this as a possible method of treatment of the 


* The calculations are taken from the paper of Jones 
and Bickford. 
18 T, Shedlovsky, J. Am. Chem. Soc. 54, 1405 (1932). 
16 T, Shedlovsky, J. Am. Chem. Soc. 54, 1411 (1932). 
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linear variation of conductivity with concen- 
tration in concentrated solution. 

It will be seen that this method of estimating 
association is somewhat similar to that used by 
Davies and his collaborators!? in more dilute 
solutions of incompletely dissociated salts. The 
cases considered here, would, however, be re- 
garded as completely dissociated on Davies’ 
treatment. 

Our estimated values of 1—a, given in Table I, 
are a direct contribution to the value of B in the 
activity equation. Although the contributions to 
B are estimated in other cases at 25°, there are 
no sufficient experimental data available at this 
temperature. In these cases, estimations of (1—«a) 
at 18°C have been used. As will be seen by com- 
paring the values of NaCl at 0° and 18°, there is 
very small temperature dependence. But it 
appears that the amount of Bjerrum association 
is an extremely difficult term to estimate. 


Variable Dielectric-Constant Correction 


If we accept the fact that this correction is 
only a contribution to the value of B and need 
not be equated to the experimentally determined 
values of this constant, it is no longer necessary 
to treat a or 6* as adjustable parameters in 
Hiickel’s theory. We have, therefore, calculated 
contributions to B on the basis of Eq. (11b), 
using Harned and Owen’s® experimental values 
of a and our experimental values of 6*. The jus- 
tification for using these values of a is that the 
dielectric constant will only vary outside the 
distance of closest approach of the ions. 


Conclusions 


In Table III are given the contributions to B 
from the various effects considered above; the 
contribution from the structure-breaking effect 


TABLE III. Valves of contribution to activity coefficient 
correction factor B. 








Electro- B B B B B 
lyte Experimental Bjerrum Hiickel Co-volume Structure 


HCl +0.13 —0.11 +009 +018  -—0.03 
NaCl +0.04 —0.10 +005 +016  —0.07 





. KCl +0.02 —0.05 +0.05 -+0.09 —0.07 


RbCl +0.02 —0.04 +0.06 +0.07 —0.08 
BaCl. +0.03 —0.12 +0.11 +0.18 —0.14 
Bexperimental = Bpjerrum+ Buiickel + Beo-volume+ Bstructure. 








17C, W. Davies, Trans. Faraday Soc. 23, 351 (1927). 
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is not calculated theoretically, but estimated by 
subtracting the other contributions from the 
experimental value of B. It should be possible 
eventually to correlate the values estimated in 
this way with data obtained by other methods, 
although only entropies, and not free energy 
data are at present available. As has been men- 
tioned above, the structure-breaking contribu- 
tions are qualitatively correct for different ions, 
i.¢., 
Bat*+ > Rbt+>K*t>Nat>Ht. 


This order would be more strongly accentuated 
if the crystallographic radii, with an allowance 
for the first hydration sheath, had been used 
instead of Harned and Owen’s empirical values 
in calculating the various contributions. 

Unfortunately the Bjerrum association for 
lithium cannot be easily estimated, as the con- 
ductivity curve is not linear with concentration 
over a wide range. It is clear that considerable 
progress must be made in study of the various 
effects before a quantitative verification of the 
behavior of activity coefficients of concentrated 
solutions can be made. Our table is extremely 
tentative and little weight should be attached to 
the numerical values. 

Many of the effects described would be capable 
of direct experimental verification. 

Accurate measurements of the Falkenhagen 
effect over a concentration range 0—1N and 
over a frequency range would be capable of 
showing whether the high frequency conductivity 
is a linear function of concentration, and, hence, 
whether Bjerrum association can be evaluated 
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from conductivity data, as has been done in this 
paper. Further experimental work on this 
problem is in the course of preparation. 

The importance of structure-breaking effects 
of the ions could be evaluated by measuring ac- 
tivities and heats of dilution at low temperatures 
where this effect should be very marked, and at 
high temperatures when it should be of much 
less importance. Further comparison of the dilu- 
tion of salts with those of rare gases, to explain 
the anomalies mentioned by Frank, are also 
indicated. 

Significant data could be obtained by carefully 
selected experiments and analyses of the existing 
data on mixtures of electrolytes where the co- 
volume effect and dielectric effect of the salts 
present would probably be the deterinining 
factors. 

It appears that future progress in electrolyte 
theory of concentrated solutions must be in an 
experimental evaluation of the various factors 
involved and by a vigorous theoretical attack on 
the basic problems of solute-solvent interaction. 
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In the usual theory of absolute reaction rates it is assumed that the molecules in the activated 
complex are in statistical equilibrium with the molecules in the initial state. This assumption 
is not justified in the case of unimolecular decompositions because the molecules in the initial 
state must gain the energy of activation by a stepwise process, usually by the transfer of quanta 
of vibrational energy from one normal mode to another. The size of these energy jumps and 
their frequency inevitably affect the over-all reaction rate. The problem is considered in con- 
siderable detail for a rough idealized mode. In this case, the actual reaction rate should be 
smaller by a factor of 0.5 to 0.387 than the value computed from the usual theory of absolute 
rates. A rigorous development would require further development of the original Ramsperger- 


Rice-Kassel treatment. 





N the usual theory of absolute reaction rates 
it is assumed that the molecules in the ac- 

tivated complex are in statistical equilibrium 
with the molecules in the initial state. This 
assumption is not justified in the case of uni- 
molecular decompositions because the molecules 
in the initial state must gain an energy of ac- 
tivation by a stepwise process, usually by the 
transfer of quanta of vibrational energy from 
one normal mode to another. The size of these 
energy jumps and their frequency inevitably 
affect the over-all reaction rate. The treatment 
given here is a generalization of the recent work 
of Zwolinski and Eyring.! 

In the case of normal chemical reactions where 
neither the initial nor the final state is uni- 
molecular, the assumption of statistical equi- 
librium is probably valid. However, the equilib- 
rium assumption may require modification if the 
activation energy is not large compared to kT, 
as in the case of some problems of viscous flow. 

Unimolecular reactions involve a special con- 
sideration which is not important in other cases. 
The way in which the systems obtain their 
energy of activation may determine the rate at 
which the decomposition takes place. The usual 
theory of absolute reaction rates must be modi- 
fied to take this into account, otherwise it would 
lead to erroneous conclusions. For example, the 
usual theory of absolute reaction rates leads to 
absurd results under the following conditions. 
Consider the possibility of a diatomic molecule 


1 Zwolinski and Henry Eyring, J. Am. Chem. Soc. (in 
process offpublication). 


AB dissociating unimolecularly without transfer 
of energy by radiation. The potential energy is 
shown in Fig. 1 as a function of the separation 
between A and B. Physically it is clear that the 
atoms A and B cannot unite nor can AB dissoci- 
ate without either a collision with another 
molecule or the transfer of energy by radiation. 
However, we can examine the problem from the 
standpoint of the theory of absolute reaction 
rates. According to all of the rules, it is permis- 
sible to consider the top of the energy hump as 
the activated state. Expressions for the rate of 
dissociation and the rate of association can be 
written down in the usual manner. The impos- 
sibility of a collision between A and B resulting 
in a stable molecule AB is expressed by setting 
the transmission coefficient equal to zero for the 
association reaction. This is altogether reasonable 
and does not require any further comment. 
However, the transmission coefficient for the 
dissociation reaction is certainly unity since the 
potential energy curve of Fig. 1 indicates that 
A and B repel each other for separations greater 
than that of the activated state. Furthermore, all 
of the other factors in the expression for the rate 
of decomposition of AB are finite and non- 
vanishing. Thus the theory of absolute reaction 
rates leads to a non-vanishing rate of dissociation 
of a diatomic molecule where it should have given 
zero. This paradox disappears when we examine 
the way in which the molecule AB gains its 
energy of activation, i.e., only by collision with 
other molecules. Since the system cannot get 
from the initial to the activated states under the 
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hypothetical conditions which we imposed, the 
normal theory of absolute reaction rates does not 
apply. Or in the language of statistical me- 
chanics, that part of phase space corresponding 
to the activated state is inaccessible to the 
molecules in the initial state so that no uni- 
molecular decomposition is possible. The accessi- 
bility of the activated state is therefore an 
important factor in reaction kinetics. 

Polyatomic molecules, in contrast to diatomic 
molecules, can decompose unimolecularly with- 
out requiring energy transfer by radiation. The 
energy surges from one degree of freedom to 
another within the molecule until the energy of 
activation is accumulated in the degree of 
freedom along which the molecule decomposes. 
The original Ramsperger-Rice-Kassel* treatment 
is still, after twenty years, the most fundamental 
approach to this problem. Ramsperger, Rice, and 
Kassel developed equations for the rate of 
decompositions in terms of the specific transition 
probability for the transfer of vibrational energy. 
These transition probabilities depend on the 
exact shape of the potential energy surface (or the 
anharmonicity of the vibrations). The lower 
vibrational states of the normal modes are almost 
harmonic, and energy transfers are infrequent. 
The higher vibrational states are no longer inde- 
pendent, energy transfers are frequent, and the 
system behaves classically. Unfortunately, no 
one has ever developed a satisfactory theory for 
the specific transition probabilities, so that it is 
not possible to use the Ramsperger-Rice-Kassel 
equations for the solution of practical rate 
problems. Instead, the Eyring theory of absolute 
reaction rates is usually applied. 

The use of the theory of absolute reaction 
rates for the unimolecular decomposition of a 
polyatomic molecule implies two assumptions: 
First, the molecule is sufficiently large that it 
may be considered as a heat reservoir possessing 
a characteristic temperature independent of the 
other molecules in the vessel. For any given 
molecule this approximation becomes better the 
higher the pressure. Thus the experimental vari- 
ation of the unimolecular rate constant with 
pressure is the best test for the validity of this 
assumption. The assumption is justified at 


? Louis Kassel, Kinetics of Homogeneous Gas Reactions 
(Chemical Catalog Company, 1932), Chapter V. 


pressures sufficiently high so that the uni- 
molecular rate constant is no longer dependent 
on pressure. The second assumption is that the 
molecules in the activated state are in equi- 
librium with the molecules in the initial state. 
This is the condition which we wish to examine 
in this present paper. For this purpose we fix 
our attention on the vibrational quantum state 
of the degree of freedom along which the molecule 
decomposes. In going from the initial to the 
activated state the molecule passes successively 
through each of the vibrational states. Sometimes 
energy is lost to other degrees of freedom within 
the molecule and this quantum number de- 
creases; at other times, energy is gained from 
the other degrees of freedom and this quantum 
number increases. Thus the process of gaining 
the energy of activation is quite devious in 
contrast to bimolecular reactions where the 
activation energy is attained in a single violent 
collision. The following calculation shows that 
the number of molecules in the activated state 
is somewhat less than would be expected from 
statistical equilibrium with the initial state if 
this energy of activation is accumulated in a 
number of steps. 

Consider a hypothetical polyatomic molecule 
in which the vibrational quantum number in the 
degree of freedom along which the decomposition 
procedes is indicated by J. Let the initial state 
be called “‘O” and the activated state ‘‘N.’’ The 
successive reactions by which the molecule goes 
from the initial to the activated state are indi- 
cated by 


k r1 k 72 TN—1 TN 


Here kr is the rate constant for molecules in the 
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activated state proceeding to the decomposition 
products. For simplicity of presentation we 
assume here that once a molecule has decom- 
posed past the activated state it does not return. 
Then we can write the set of simultaneous equa- 
tions for the rate of change of the concentrations, 
Cy, in each of the vibrational states. After a 
steady state is reached, the rate of change of 
each of these concentrations (except the Co) is 
zero so that: 





dC4 
Y =kjyiCni—(kr +krv)Ca=0, (2) 
dCny-1 
ae =kjy_2Cnv_2— (Rinatkrn1)Cy-1 
+krvnCa=0, (3) 
dCy 
— =kRypyiCyi— (kts t+krs) Cr 
thryyiCy41 =0, (4) 
dC; 
— =k aCe— (kath) Crt bnCe=0. (5) 
( 


These equations suffice to determine C4 in terms 
of Cy. Let us examine the simplest possible case 
where 


kp=ka=kRp= sada =ksy_41=k,y, (6) 
kr =kre=kr3= ++ -Rrvn=k,, (7) 

and 
k;/k-=exp(—e/kT). (8) 


This situation would occur if the energy jump 
between successive vibrational states is e if all 
vibrational states have the same degeneracy, and 
if the transmission coefficients are in all cases 
unity. The solution to Eqs. (2) through (5) is 
then :* 


c. 1 
(Ca)en 2 el 
k,\ 1—exp(—e/kT) 








Here, (C4)xq is the concentration of molecules in 
the activated state which would exist if statis- 
tical equilibrium were attained between the 
molecules in the initial state and the molecules in 


*The author wishes to thank the Journal of Chemical 
Physics reviewer for pointing out this simple form for the 
solution. 
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the activated state. For the model which we are 
assuming, it is clear that 


(Ca)Eq = Co exp(— Ne/kT). (10) 


Furthermore, with this model we can determine 
the ratio kr/k,. At complete statistical equi- 
librium, including back reactions, the same 
number of systems are moving backwards as 
forwards through the activated state. Since we 
have assumed that all transmission coefficients 
are unity, it follows that kr/k,=1. If, in addi- 
tion, the energy of activation for the reaction, 
Ne is large compared to kT, we obtain: 


Ca 1—exp(—e/kT) 
(Ca)eq 2—exp(—e/kT) 





(11) 


Here if ¢ is large compared to RT, 
Ca/(Ca)Eq=?- (12) 
If € is equal to kT, 
C4/(Ca) ey = (€—1)/(2e —1) = 0.3873. (13) 


And if ¢ is small compared to kT, Ca/(Ca)xq 
becomes correspondingly small. However, it is 
likely that if the energy between successive vibra- 
tional states were small compared to kT, the 
motion would become classical and in the degree 
of freedom along which the molecule decomposes, 
the average energy change on successive surges 
would remain of the order of kT. We conclude 
that for this simple model of a hypothetical 
polyatomic molecule the concentration of mole- 
cules in the activated state is roughly between 
0.5 and 0.38 times the value expected on the 
basis of statistical equilibrium. Thus, for this 
case, the rate of the unimolecular decomposition 
is smaller by this same factor than would be 
expected on the basis of the usual theory of 
absolute reaction rates. 

The model considered in this paper is ad- 
mittedly crude, but it is gratifying in that it 
leads to a concentration of molecules in the 
activated state only slightly smaller than the 
equilibrium concentration assumed in the usual 
theory of absolute reaction rates. This superficial 
examination does not uncover any evidence that 
the theory of absolute reaction rates would lead 
to any large errors in the rate of unimolecular 
decompositions of polyatomic molecules. How- 
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ever, the use of this theory leaves much to be 
desired from the logical standpoint. Real 
progress in the theory of unimolecular reactions 
requires the evaluation of the specific transition 
probabilities in the Ramsperger-Rice-Kassel 
treatment. Until this is done, the theory of uni- 
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molecular decompositions must remain on a 
semiquantitative basis. 

The author wishes to express his thanks to the 
Navy Bureau of Ordnance for financial assistance 
during the course of this work through NOrd 
Contract 9938 with the University of Wisconsin. 
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‘ 


The statistical method of treating the sorption of vapors, begun by Fowler and Guggenheim 
and extended to multilayer sorption by Cassie and Hill, is further generalized to include vari- 
able heats of sorption in different layers. It is shown under what conditions multilayer sorption 
can lead to a linear isotherm and to various modifications of the BET equation. The treatment 
is also extended to the case of interaction between sorbed molecules in the same layer. 


N a previous paper! methods of calculating the 

free energy and heats of sorption of water 
vapor by proteins and high polymers were dis- 
cussed for the cases of linear, Langmuir, and 
BET? sorption isotherms. At first sight it might 
appear unsound to apply equations derived on 
the basis of a kinetic picture of equilibrium 
between gas molecules colliding with and evapo- 
rating from a solid surface to the sorption of 
water by solids such as proteins and high 
polymers where the water molecules must diffuse 
through the solid to reach the points of sorption, 
but when it is realized that Fowler and Guggen- 
heim* and Cassie‘ and Hill’ have derived both 
the Langmuir monolayer adsorption and the 
BET multilayer adsorption equations on the 
basis of statistical reasoning, the theoretical 
justification for this procedure becomes apparent. 

At this point we wish to emphasize that appli- 
cation of the BET equation to the sorption 


1M. Dole and A. D. McLaren, J. Am. Chem. Soc. 69, 
651 (1947). 

2S. Brunauer, P. H. Emmett, and E. Teller, J. Am. 
Chem. Soc. 60, 309 (1938). 

’R. H. Fowler and E. A. Guggenheim, Statistical 
a (University Press, Cambridge, 1939), p. 
*A. B. D. Cassie, Trans. Faraday Soc. 41, 450 (1945). 

5T. L. Hill, J. Chem. Phys. 14, 263 (1946). The treat- 
ree of Hill is a correction to and extension of that of 
assie, 


isotherms of water sorbed by high polymers in 
no way implies any estimate of a surface. In 
fact, by a careful study of the statistical deri- 
vation it is realized that the BET equation gives 
as one of its parameters the number of adsorp- 
tion sites, not the surface. The surface can only 
be obtained by a further calculation in which a 
value for the area per site is assumed. Thus the 
area is a derived quantity. By adsorption we 
mean adsorption on a free surface, where the gas 
molecule can collide directly with the adsorbing 
site, while by sorption we mean the increase of 
weight of the solid as the result of the penetration 
and sorption of water vapor within the solid, at 
sorption sites where there is no free surface. 
Justification of the application of the BET 
equation to the sorption of water vapor by high 
polymers depends upon the existence of a real 
reversible equilibrium, no hysteresis, upon the 
negligible contribution of change of volume to 
the free energy of the sorption process, upon the 
sorption sites of the polymer remaining fixed in 
space during sorption (otherwise the process 
would be that of mixing two liquids), and upon 
the absence of interaction between sorbed 
molecules. In a later paper giving new data on 
the sorption of water by certain experimental 
polymers we shall discuss the applicability of the 











BET equation more fully and compare the 
theory presented here with that of Barrer.® 

In this paper we extend the statistical treat- 
ment to consider cases of varying heats of sorp- 
tion between layers and of multilayer sorption 
with interaction between molecules in the same 
layer. The nomenclature of Fowler and Guggen- 
heim will be used for the most part. 


GENERAL STATISTICAL TREATMENT 
Let 


So=number of sites having no sorbed molecules, 
S;,=number of sites having one sorbed molecule, 
S:=number of sites having two sorbed mole- 
cules, etc., 
N,=total number of sorbed molecules in the 
first layer, 
N.2=total number of sorbed molecules in the 
second layer, etc., 
Ngs=total number of sorption sites, 
N=total number of sorbed molecules, 
a,=internal partition function for sorbed mole- 
cules in the first layer, 
a2= internal partition function for sorbed mole- 
cule in the second layer, etc., 
a,=internal partition function for molecule in 
the liquid state, 
x=p/po, partial pressure ratio of sorbate in 
the vapor phase, 
Q=total partition function of all sorbed mole- 
cules, and 
n =total number of layers. 


As a result of these definitions, 


So=Ns—Mi, S;=Ni—Nsz, etc., 
> N;=N. 
i=1 


The total partition function is given by the 
equation 


Ns! » 
Q=—— TI (a)™ (1) 
Os! 
0 
or 
Ns! 
Q 





~ (Ns—Ni) (Ni— Ns) !---(Naa—N,) IN! 


een -a,Ni.- . “An. (2) 


6 R. M. Barrer, Trans. Faraday Soc. 43, 3 (1947). 
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The Helmholtz free energy, F, then follows from 
(2), after approximating the factorials by Stir- 
ling’s theorem. 


F 
——=InQ=Nsg InNs—(Ns—N;) In(Ns—™) 
RT 
—-+++(Ny-1— Np) In(Nn-1— Nn) 


—N,1InN,+> N,lna;. (3) 


i=1 


For the conditions of equilibrium between the 
first and second layers, similarly to Hill,? we 
obtain 


Ns—N; Ni—N2e ao 
Ni-N2 No—N; a, 





(4) 


and for equilibrium between the second and 
third layers, etc., 


Ni-—N2 No—Ns as 
No—N; Ns—N, do 





etc. (5) 


In Hill’s treatment 
a3 = a2, 


with a consequent simplification of (5) and 
similar relationships. 

Finding the variation of the free energy with 
respect to NV and setting the result equal to the 
chemical potential of the pure liquid phase, we 
deduce that 

Na Qn 


——_ez.—. (6) 
Nr-si— Nn ay 


Substituting (6) into (4),-(5) and similar rela- 
tions, we find, also similarly to Hill, 


a 
Ni—Ne=(Ns—Nj)x-—, 





ay 
x2 
N2—N3=(Ns—Ni1)— 41:2, 
az? 
xn 
Nnr=(Ns—M1) Q1°* "An, (7) 
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except that in Hill’s case (the BET model) 
Qe=A3=—°*** =An=AQL. 


Making use of (7), the generalized sorption 
equation results for the model used here, that is, 
localized sorption on definite sites with no inter- 
action between molecules on adjacent sites, 
namely, 


N  y(ait2yai-de+3y7a1-d2-a3+---) 





aie ee ’ 8) 
Ns 1+9a1+3*d1-d2+y*d1*d2-d3+°-* ; 
where 
x 
ja, 
ay 
Letting 


@=1+yaityaid2+---, 
it is readily seen that 
N ¢ ; 
Wee (9) 
SPECIAL SOLUTIONS 
1. Partition Functions in All ‘Layers Equal 


The simplest assumption is that the a,’s are 
all equal, and are some fraction, k, of the pure 
liquid-state partition function @,. For an infinite 
number of layers the solution is 





N kx 
Ns 1—kx 
and for a finite number n, 
N 1—(n+1)(kx)"+n(kx)"*! 
aii J] (11) 





Ns (1—kx)[1—(kx) #7] 


In case k is equal to unity, which would be true 
for an ideal solution, then (10) reduces to 


N=x, (12) 


which is Raoult’s law (N is the mole fraction of 
the vapor in the solid). 


2. Partition Functions All Equal Beyond the 
First Layer 
If 


a1 
—=42=—-a3= eee =Q,=-<aL1, 
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the BET equations result. 
If 


d2=ka,=a3=a4=--- =A,¥F aL, 
and if m is limited to a finite number, we obtain 


N y cx1—(n+1)(ckx)*+n(ckx)"**] 








Ns [1-—chkx][1+cx{1—k—(ckx)"}] 
where ¢=4a;/az,. When n goes to infinity 
N XC 
(14) 


Ng [1—ckx][1+cex(1—&)] 


which becomes identical to the BET equations 
on setting k=az/a;. Equation (14) is closely 
similar to that deduced by Anderson’ (Ander- 
son’s Eq. (6)) which in our terminology can be 
written 


N cxk 
en i (15) 
Ns [1—kx][1+(c—1)kx] 





When Egs. (14) and (15) are put into a linear 
form and Solved for Ns and c, it can be seen that 
in Anderson’s case 


1 
Ns=——, 
Tk+S 
where J and S are the intercept at x equal to 
zero and the slope, respectively, and in our case 





1 
Ns = . (16) 
Ickh+S 
As the functions plotted are 
x x : 
—————- against x and -———— against x 
N(1—ckx) N(1—kx) 


in the two cases, it is obvious that both Ander- 
son’s and our treatment will yield equal values 
of Ns. However, if the constant c is solved for 








Ick+S 
C= 
I 
in our case, and 
Ik+S 
C= 
Ik 


7R. B. Anderson, J. Am. Chem. Soc. 68, 686 (1946). 











in Anderson’s case, it is evident that a difference 
results. As Anderson derived his equation for a 
heat of adsorption in the second to the ninth 
layer different from that in the first layer as well 
as different from the heat of condensation in the 
pure liquid, one would not expect his and our 
equation to be identical (Anderson did not give 
details of derivation). 

It should be noted that our equation reduces 
to the Langmuir isotherm on setting k equal to 
zero, whereas Anderson’s equation does not. 


3. Partition Functions Decreasing Harmonically 


Let us make the arbitrary assumption that 








d2=1/2a,, a3=1/3a, etc., (17) 
then 
ya ya; 
ta FF -*, 
2! 3! 
= eye 
and 
N a, 
—=x-—=Cx. (18) 
Ns ay 


This demonstrates that given the proper decrease 
in the internal partition function with increasing 
number of sorbed layers, a linear sorption 
isotherm can be obtained. In the case of the 
sorption of water by high polymers there are a 
number of examples of the sorption following a 
linear relationship with the pressure up to 50 to 
60 percent relative humidity. There is no a priori 
reason in these cases why the partition function 
should not fall to zero in the limit. In fact, the 
Langmuir isotherm can be obtained from (8) by 
setting all partition functions except a: equal to 
zero. 


MULTILAYER SORPTION WITH INTERACTION 


Fowler and Guggenheim® have shown how the 
Langmuir treatment is to be modified if inter- 
action between adjacent molecules forming a 
number of pairs, Ni, in any layer is taken into 
account. If the contribution to the internal par- 
tition function for such pairs is a,;, then the 
total partition function for multilayer sorption 


8 See reference 3, p. 429. 
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becomes 


Ns! 1 
Q= ona (ai) ¥*(ayi)¥® (19) 


i=1 


II S;! 
i=0 


where N;;’ is a statistical average number of 
pairs when N; molecules are adsorbed. Making 
the assumption of a completely random arrange- 





ment, Fowler and Guggenheim deduce the 
equation 
N?? 
N;,/ =42 ’ (20) 
2 Ng 


where z is the number of nearest neighbors in any 
layer. The contribution of the interaction 
energy, W, to the total energy can be expressed 
by the equation, 

N 2 


a 


Fint. _ 





W. (21) 


Ns 


(This follows from the assumption that the 
excess energy resulting from interaction can be 
expressed by the factor N;i2W/z].) 

The free energy of the adsorbed material and 
other relations follow immediately from (19) 


F 
—=-—NsgInNs+(Ns—N;j) In(Vs— Mi) 
RT 
+(Ni—WN-2) In(Ni—N2)+-::: 
+N, Inv,.—-¥ N; Ina; 


i=1 














n N?? Pa 
—> gina Ind;;, (22) 
i=1 Ns 2 
Bax Ni-—N:; de Ag9X2/N 8 
- pi (23) 
Ni—Nz Ne—Nsa au%¥!/"s 
Ns—Ni Na-1— Na Gafan®*7'%5 
_ “| | (24) 
Ni—N2 Nw aylaiy™) 
Na an 
on =X—OnnNa- (IN 8), (25) 
Nn-1—Na aL 
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2 
ai x ; ; 
vaste) +2| =] 4°81 G/N S) © @g-ogh 2" G/N s)4 ae 





Ns a, 


(26) 
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Fowler and Guggenheim® point out that the 
assumption of interaction leads at certain values 
of the temperature to the splitting of the ad- 
sorbed layer into two phases. Hill®® suggests 
that in the case of multilayer sorption this 
region of instability, corresponding to a sudden 
jump in the value of N/Ns, is to be conceived 
physically as capillary condensation, and he 
concludes that the multilayer adsorption and 
capillary condensation points of view are there- 
fore really not in conflict. 

In the case of the sorption of water vapor 
by solids, the present author prefers to consider 
the sudden increase in N/Nsg at constant vapor 
pressure as analogous to compound formation 
such as occurs in the hydration of crystals when 
the percentage of water in the solid (water of 
hydration) rises from a low to a high value 
without change of vapor pressure. Such hydra- 
tion may occur occasionally with very little 
change in the x-ray diagram of the solid.'” 

Let us consider the application of Eq. (26) to 
the phenomena of hydration of crystals. When a 
higher and lower hydrated form of the crystal 
are simultaneously in equilibrium with a certain 
vapor pressure of water, then the percentage of 
water in the crystal will rise without change in 
the chemical potential, u, of the water, or 


Ou 


— (27) 


It is interesting to find the value of N, corre- 


*In a paper in press Dr. Hill has gone considerably 
deeper into the mathematics of the unstable region than 
he did in the paper quoted in reference (5). We are greatly 
indebted to him for the opportunity of seeing his paper in 
advance of publication. 

10See, for example, the interesting paper by W. O. 
Milligan and H. B. Weiser, J. Phys. Chem. 41, 1029 
(1937), where x-ray diffraction patterns are given for a 
number of hydrated zeolites. In their case they heated 
the hydrated crystals at constant vapor pressure and 
observed a stepwise decrease in the percentage of water in 
the crystal. 








sponding to condition (27). Letting N,,° represent 
this value, we obtain 





An 

1 x—agg 0" INS) 
ay 

N,° = ’ (28) 
q 
where 

Z 

q oo Ina, ne 

Ns 


Differentiating N,° with respect to x, we find 


ON n? 


An*Annr™ GINS) 


= ——_________—_.,__ (29) 


Ox An | : 
q:a L 1 = i—B on ad n:(z/Ng) 
ay 


so that for dN,°/dx to be equal to infinity as it 
must in the case of formation of hydrates at 
constant vapor pressure, 


an 
NX> —OAnyN 2’ @/NS) =1., 


ay 


(30) 


Hence 
2 


N,°=- (31) 


For the layer immediately under the outer layer 
we obtain, on rearranging (25), 


Qn 
1+-x—ann® @/NS) 
aL 
Na-=Nr ’ 
an 
X*—Ona®™’ (z/N s) 


ay 





which for V,=N,° becomes 
Na-1=2N,’. 


If N,° equals 1/2N 5, which it will at the critical 
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temperature when 


W = —2kT 
or 
4 
q Ne’ 
then 
N,-1=Ns. 


In other words, the sorption in the next layer 
below the outer layer is complete when the top 
layer is half-covered, as would be expected when 
the hydration of a crystal proceeds in steps. It 
would appear, then, that the physical fact which 
distinguishes between sorption of the types 
represented by the Langmuir, BET, or other 
equations and hydrate formation is the inter- 
action energy between the sorbed molecules. 





BERNSTEIN AND G. HERZBERG 


When the sorbed molecules are so far apart that 
no or insufficient interaction energy exists, no 
region of increase of water sorbed at constant 
vapor pressure can exist, no separation of the 
solid into two phases will occur; on the other 
hand, a gradual increase of the amount sorbed 
with increase of vapor pressure will be observed. 
Application of the equations contained in 
this paper to new data on the sorption of water 
by certain high polymers will be carried out at 
a later date. The work given here is part of a 
general program of study of the thermodynamics 
of high polymers for which grants have been 
gratefully received from the du Pont Corpora- 
tion, the Research Corporation," the Richardson 
Company, and the Visking Corporation. 


1 Frederick Gardner Cottrell grant. 





THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 16, 


NUMBER 1 JANUARY, 1948 


Rotation-Vibration Spectra of Diatomic and Simple Polyatomic Molecules with Long 
Absorbing Paths. 


I. The Spectrum of Fluoroform (CHF;) from 2.4u to 0.7u 


H. J. BERNSTEIN* AND G. HERZBERG 
Yerkes Observatory of the University of Chicago, Williams Bay, Wisconsin 


(Received August 25, 1947) 


The spectrum of fluoroform has been investigated under high resolution (21-ft. grating) in 
the photographic infra-red with an absorbing path of up to 60 meters, obtained by multiple 
reflection according to the method of J. U. White. In addition, with the same path length the 
spectrum from 1.2 to 2.44 was obtained under low resolution with a photoelectric infra-red 
spectrometer. A large number of overtone and combination bands, 47 in all, were found. Two 
parallel bands in the photographic region were well resolved and analyzed, yielding, for the mo- 
ment of inertia about an axis perpendicular to the symmetry axis, 


Tp =81.08 X 10- g cm? 


From this value, assuming tetrahedral angles and a C—H distance as in methane, the C—F 
distance in fluoroform is found to be 1.329A, a value that is appreciably lower than that found 
in methyl fluoride, 1.380A. A provisional vibrational assignment of all the bands is given. 


A. INTRODUCTION 


HE investigation of the rotation-vibration 
spectra of diatomic and polyatomic mole- 

cules in the photographic infra-red by means of 
large diffraction gratings makes possible a deter- 
mination of molecular constants (moments of 
inertia, internuclear distances, etc.) of very high 


* Division of Chemistry, National Research Laboratories, 
Ottawa, Canada. 





precision. Indeed it was only during the last 
year that this precision has been surpassed by 
the data obtained from microwave-absorption 
spectra. The microwave technique is applicable 
only to a certain group of molecules and to 
certain molecular constants. Therefore, it re- 
mains imperative to continue the study of photo- 
graphic infra-red spectra. It appears that there 
are still a number of diatomic and simple poly- 
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atomic molecules which ought to be studied in 
the photographic region, and others for which 
the information available at present could be 
greatly extended with improved technique. The 
present series of investigations is directed at 
this aim. 

The study of molecular absorption spectra in 
the photographic infra-red entails the use of 
long absorbing paths. Particularly long paths 
are required for molecules not containing hydro- 
gen, since their vibrational frequencies are small 
and therefore only rather high overtones fall into 
the photographic region. Indeed up to now only 
for one molecule of this type (CO2) have photo- 
graphic infra-red spectra been obtained.! 

In most infra-red work interest is centered on 
the fundamentals, and only comparatively short 
path lengths are used. If long paths are used 
many overtone and combination bands may be 
expected in the near infra-red as well as in the 
photographic region. It seems important to 
locate as many of these overtone and combina- 
tion bands as possible in order to obtain a reliable 
representation of all the vibrational levels of the 
molecule under consideration and in order to 
study its potential function. Toward this end, 
in the present series low dispersion spectra up to 
2.4u have been taken with the same long path as 
for the photographic high resolution spectra. It 





(a) 


is hoped that later on it will be possible to study 
also this spectral region with much _ higher 
resolution. 

A precise knowledge of the higher vibrational 
levels, as well as of the rotational constants, is 
not only of interest in the study of molecular 
structure, but also for the purpose of calculating 
thermodynamic functions particularly at higher 
temperatures. This need gives added stimulus 
to the investigations of infra-red spectra with 
long absorbing paths. 

The present first paper in this series gives the 
results of the investigation of the spectrum of 
fluoroform. Spectra of ketene and carbon dioxide 
have already been taken and are being evaluated. 
The study of other molecules is being prepared. 


B. EXPERIMENTAL 


A long absorbing path was obtained by 
multiple reflection in a comparatively short tube 
according to the ingenious method suggested by 
White,” with slight modifications. Figure 1 shows 
the optical arrangement. Three concave mirrors 
of equal radius of curvature are used. Two of 
these, A and B, which were cut from one circular 
mirror as shown in Fig. 1b, are mounted at one 
end of the tube while the third, C, of a shape 
indicated in Fig. 1c, is mounted at the other 
end at a distance equal to the radius of curvature. 














(b) 











Fic. 1. Optical arrangement for obtaining a long absorbing path. 


Herzberg and Verleger [Phys. Rev. 48, 706 (1935) ]. 
* J. U. White, J. Opt. Soc. Am. 32, 285 (1942). 


‘In the spectrum of Venus by Adams and Dunham [Astron. Soc. Pac. 44, 423 (1932)] in the laboratory by 
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Fic. 2. Spectrogram of the CHF; bands at (a) 1.1370u and (b) 1.1639y. 


Light from an automobile headlight bulb is 
focused on an entrance slit at 0 in Fig. 1a and c. 
From there the light falls on the two mirrors, 
A and B, which form images of the slit at 1 and 
—1, respectively. The image —1 is discarded 
but is very useful in the adjusting procedure. 
The mirror C is so adjusted that the light from 
A received at 1 is reflected to B. With proper 
adjustment A is imaged on B and no light is 
lost except by reflection losses. The mirror B 
then forms an image of 1 at 2, whereupon the 
light is reflected to A, and so on, until the light 
emerges at 8 after having traversed the space 
between the mirrors 16 times. By turning mirror 
A by small amounts about an axis perpendicular 
to the plane of the paper of-Fig. 1a, the number 
of images on mirror C, that is, the number of 
traversals, can be readily changed. In order to 
accommodate as many images as possible on 
mirror C and for the purpose of simplifying the 
adjustment, the even and odd images are sepa- 
rated (by the proper adjustment of mirrors A 
and B) as shown in Fig. 1c. This diagram also 
makes it clear why the mirror C has been 
indented at the top half as shown. The separation 
of the axes (dot-dash lines) of the two mirrors 
A and B at C is half the separation of successive 
odd (or even) images. It is not necessary to 
have the odd images halfway between the even 
ones. In fact if the adjustment is made by 
turning mirror A only, such a symmetrical 
arrangement of the images will obtain only for a 





certain number of images depending on the 


orientation of mirror B. 

Thus far two long paths as just described have 
been set up: a small one in a glass tube of 50 mm 
diameter, and with 150 cm radius of curvature 
of the mirrors, and a large one in a galvanized 
iron pipe of 10 inches diameter and with mirrors 
of 22 m radius of curvature. With the small 
tube up to 40 traversals have been obtained, 
with the large one up to 144 traversals. The 
resulting path lengths are 60 m and 3168 m, 
respectively. The limit is set only by the reflec- 
tion losses and by the number of images that 
can be accommodated on mirror C. The spectrum 
of fluoroform has been taken entirely with the 
small absorption tube, since not enough of the 
gas was available for the large one. 

The photographic infra-red spectra were taken 
in the first order of a 21-ft. grating ruled by R. 
W. Wood. This grating has an excellent concen- 
tration in the first order of 11,000A and is 
thus particularly suitable for work with the 
Eastman I-Z plates. The grating is mounted 
in a Eagle mounting. A region of 2500A in first 
order can be covered in one exposure. The 
dispersion is 2.5 A/mm. While below 9000A the 
exposure times (on I-N plates) were of the order 
of a few minutes for 20 traversals, at 11,000A the 
exposure times were of the order of 24 hours. 

Spectra in the region from 1 to 2.4y were 
taken with the new infra-red spectrometer built 
by Dr. G. P. Kuiper for work with the McDonald 
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telescope.* In it a photo-conductive cell is used 
as the receiving element. The resolving power 
is 800. The spectrum is recorded on an Esterline- 
Angus recorder. Depending on the speed of 
recording, the dispersion on the record varies 
from 20 to 100 A/mm. 

The fluoroform was obtained through the 
kindness of Dr. A. F. Benning from the Jackson 
Laboratory of E. I. du Pont de Nemours and 
Company. It was specified as ‘‘practically pure’. 
Pressures up to 1 atmos. were used for the weaker 
bands. However, for a satisfactory resolution of 
the fine structure of the intense bands the use of a 
lower pressure (about 4 atmos.) proved of great 
advantage, since the pressure broadening is much 
reduced at these lower pressures. 


C. ROTATIONAL ANALYSIS 


A considerable number of photographic infra- 
red bands were found (see below) but only two 
of these, at 1.1639u and 1.1370y, were sufficiently 
strong so that their fine structures could be 
measured. In Fig. 2 spectrograms of these two 
bands are reproduced as obtained at a pressure 
of 12 cm and a path length of 18 meters. Both 
bands consist of a strong central Q branch with 
a sharp edge at the long wave-length side and a 
P and R branch on either side of it. In addition, 
in each band a secondary Q branch occurs which 
is obviously caused by a difference band forming 
the 1—1 band of a sequence starting with the 
main band (see below). The band 1.1370, is 
partially overlapped by the water band at 1.135. 

The fluoroform molecule is an oblate sym- 
metric top with the figure axis in the direction 
of the C—H bond. There are parallel and 
perpendicular bands (with AK=0 and +1, re- 
spectively). There can be no question that the 
two bands reproduced in Fig. 2 are parallel 
bands, since the central branch has a sharp edge 
and since there is only one P and one R branch 
in each band. To be sure, on account of the 
comparatively small difference between the mo- 
ments of inertia Iz and I¢, the Q branches of the 
sub-bands of a perpendicular band will also form 
a fairly pronounced central branch; but it will 
not have a sharp edge, and the P and R branches 
of different sub-bands will not form apparently 





*G. P. Kuiper, W. Wilson, and R. J. Cashman, Astro- 
phys. J. 106, 243 (1947). 





single branches of the whole band as they do 
for parallel bands. It should be emphasized that, 
on account of the variation of the rotational 
constants B and C with the vibrational quantum 
numbers, parallel bands in the photographic 
region do not necessarily have this simple struc- 
ture. Only when the effective moments of inertia 
differ very slightly in the upper and lower states, 
can such a simple structure be expected.‘ In the 


TABLE I. Wave numbers of the band lines. 











Band at 1.16394 Band at 1.1370 
Vhead(Q) =8589.30 Vhead(O) =8792.75 
P P(J) R(J) P(J) R(J) 
4 8592.72 8796.12 
5 93.33 96.73 
6 8585.00 93.97 8788.52 97.42 
7 84.37 94.66 87.86 98.18 
8 83.64 95.35 87.16 98.84 
9 82.94 96.02 86.41 8799.52 
10 82.24 96.64 85.72 8800.18 
11 81.52 97.31 85.03 00.90 
12 80.79 98.08 84.32 01.54 
13 80.10 98.76 83.60 02.14* 
14 79.45 8599.43 82.91 02.84 
15 78.79 8600.07 82.22 03.52 
16 78.02 00.73 81.50 04.20 
17 77.28 01.38 80.78 04.86 
18 76.60 02.06 80.09 05.53 
19 75.86 02.68 79.33 06.15 
20 75.13 03.30 78.62 06.83 
21 74.38 03.98 77.88 07.48 
22 73.69 04.64 77.15 08.15 
23 72.97 05.29 76.46 08.78 
24 72.26 05.94 75.73 09.38* 
25 71.52 06.53 75.00 10.09 
26 70.79 07.18 74.27 10.74 
27 70.02 07.81 73.56 11.43 
28 69.31 08.46 72.79 12.00* 
29 68.56 09.06 72.18* — * 
3 67.78 09.70 71.40* 13.38* 
31 67.06 10.32 70.60 13.95* 
32 66.31 10.90 69.86 14.48 
33 65.52 11.50 69.12 15.14 
34 64.79 12.16 68.34 15.78 
35 64.06 12.73 67.63 16.39 
36 63.28 13.37 66.86 16.98 
37 62.54 13.97 66.08 
38 61.73 14.54 65.37 
39 60.98 15.10 64.59 
40 60.20 15.74 63.89 
41 59.46 16.29 63.13 
42 58.71 16.91 62.23 
43 57.92 61.58 
44 57.11 60.76 
45 56.40 60.10* 
46 59.33 
47 58.55 
48 57.68 
49 56.96 
50 56.20 














* An asterisk indicates overlapping by H20O lines. 


4See. G. Herzberg, Infrared and Raman Spectra of 
Polyatomic Molecules (in future referred to as J.R.S.), (D. 
Van Nostrand Company, Inc., New York, 1945), p. 418. 
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present case this condition is fulfilled since the 
vibration of the H atom which is mainly involved 
(see Section D) changes the moment of inertia 
but slightly. 

In Table I the wave numbers of the lines of 
the P and R branches, as well as those of the Q 
heads of the two bands, are given. The num- 
bering of the lines was readily obtained since 
the head of the Q branch gives very nearly the 
position of the band origin. . 

Table II gives the combination differences 
AoF” (J) =R(J—1)—P(J+1) for the two bands. 
They agree for each J value within the accuracy 
of the measurements. It should be mentioned 
however, that an almost equally good agreement 
is obtained if the numbering in both bands is 
shifted by one unit. The reason for this insensi- 


TABLE II. Combination differences. 














; AoF’"(J) AoF" (J) 

from 1.13704 from 1.1639. AsF’"(J) 

P band band mean 
5 7.60 7.72 7.66 
6 si 8.87 8.96 8.92 
7 10.26 10.33 10.29 
8 11.77 11.72 11.75 
9 13.12 13.11 13.12 
10 14.49 14.50 14.50 
11 15.86 15.85 15.86 
12 17.30 17.21 17.26 
13 | 18.63 18.63 18.63 
14 | (19.92) 19.97 19.97 
15 | 21.34 21.41 21.38 
16 22.74 22.79 22.77 
17 | 24.11 24.13 24.12 
18 | 25.53 25.52 25.53 
19 26.91 26.93 26.92 
20 | 28.25 28.30 28.28 
21 29.68 29.61 29.65 
22 31.02 31.01 31.02 
23 | 32.42 32.38 32.40 
24 33.78 33.77 33.78 
23 | (35.11) 35.15 35.15 
—s | 36.53 36.51 36.52 
27 | 37.95 37.87 37.91 
28 =| ~— (39.25) 39.25 39.25 
29 =| ~— (40.60) 40.68 40.68 
30. on 42.00 42.00 
31 (43.52) 43.39 43.39 
32 | (44.83) 44.80 44.80 
33 46.14 46.11 46.13 
34 47.51 47.44 47.48 
35 48.92 48.88 48.90 
36 50.31 50.19 50.25 
37 51.61 51.64 51.63 
38 52.99 52.99 
39 54.34 54.34 
40 55.64 55.64 
57.03 57.03 

58.37 58.37 

59.80 59.80 














Values in brackets are from blended lines (* in Table I) and have 
not been used in forming the mean. 





BERNSTEIN AND G. HERZBERG 


tivity of the agreement of the combination 
differences to the numbering is due to the fact 
that the rotational constants B’ in the upper 
states of the two bands are almost the same. 
It appears, however, practically certain that the 
numbering as obtained from the positions of the 
Q branches is correct. If the effect of centrifugal 
distortion is included’ the rotational term values 
are given by 


F(J, K) =BJ(J+1)+(C-—B)K?—D,J*(J +1)? 
—DyxJ(J+1)K?-—DxK*‘, 
where 
h h 
B= — and C=—— 
8rrclp 8rclo 





and D,;, Dyx, and Dx are very much smaller 
than B and C. From this formula it follows 
immediately that 


AoF(J) = F(J+1, K)—F(J—1, K) 
= (4B—6D, —4D jx K*)(J+3) 
—8Dj(J +34)’. 


If the effect of centrifugal distortion is neglected 
A>F(J)/(J +4) would be a constant, 4B. Actu- 
ally, as shown by Fig. 3, the observed A2F’’(J)/ 
(J+) values show a slight decrease with (J+ 4)? 
indicating the effect of the term 8D;(J+})’. 
From the slope of the straight line in Fig. 3 one 
obtains D,’’= 3310-7 cm™ while the intercept 
on the ordinate axis gives 4B” —6D 7". —4D yx" K° 
= 1.3806, cm~!. The term —6D,""—4Dyx"'K? 
may be neglected in view of the magnitude of 
D,;" and the consideration that D;x’’ would be 
expected to be of the same order as D,’’. Thus 
we have 
BY, = 0.34515 cm—. 


The rotational constants B’ for the upper 
states can be more accurately determined from 
the B’—B” values obtained by plotting R(/J—1) 
+P(J) against J?. The same plots also give the 
band origins vo. The resulting B’ and yo values 
are given together with the other rotational 
constants in Table III. The band origins vo agree 
very well with the positions of the Q heads given 
in Table I. The latter contain, of course, the 
error inherent in measurements of band heads. 

It may be mentioned that the R(J—1)+P(J) 


versus J* curves show a slight curvature at high 


5 See G. Herzberg, I.R.S., p. 26. 








VI 


ve 











er 
mM) 


he 
es 
al 
ce 
en 


he 


J) 
gh 








SPECTRUM 








OF FLUOROFORM 





J J 





7 | 
- °5 500 


1000 1500 (+h) 


Fic. 3. Combination differences A2F’’(J)/(J+ 4). 


J values. This is due to the fact that each “‘line”’ 
is actually a superposition of J+1 lines corre- 
sponding to K=0, 1, 2, ---J which spread out 
more and more to longer wave-lengths as J 
increases. Therefore in determining B’—B” from 
the slope of these lines only the first part of these 
lines up to J/=30 was used. 

The rotational constants B depend on the 
vibrational quantum numbers 2; according to 


Bio) =B,-— dYai(v;+d;/2), 


where d;=1 or 2 depending on whether the 
vibration v; is non-degenerate or doubly degen- 
erate. According to the vibrational analysis in 
Section D, the two photographic bands 1.1370u 
and 1.16394 are 3v,; and 2»;+2v4, respectively. 
Therefore the B’—B” value of the former gives 
3a; while that of the latter gives 2a,:+2a,. 
In this way one obtains a,;=0.00025 and 
a4=0.00021. Since the other a,’s have not been 
determined, it is not possible to ebtain B,, the 
rotational constant for the equilibrium position. 

The Bio; value determined above gives the 
effective moment of inertia in the lowest vibra- 
tional level. One obtains, using Birge’s*® con- 
version factor 27.9865 x 10-*°, 


Tl = 81.08 10-*° ¢ cm?. 





°R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 


On account of the effect of Coriolis interaction 
a determination of the second moment of inertia 
I-® would require the rotational analysis not 
only of one but of all three perpendicular funda- 
mentals or of a perpendicular vibration in the 
Raman effect. Since such analyses have not been 
carried out, one must be satisfied with the 
determination of one unknown internuclear dis- 
tance from J;"!, assuming plausible values for 
the others. 

It appears very probable that the C—H dis- 
tance in CHF; is closely the same as in methane, 
that is 

r(C-H) = 1.0936 X 10° cm. 


But even if a value as found in ethylene (1.071) 
or acetylene (1.0597) would be assumed, the 
following conclusions would be substantially 
unaltered. More critical is the assumption made 


TABLE III. Rotational constants of CHF3. 








Boy 0.34516 cm! 

Bie for 2vi +25 0.34423 cm7! 

Bix) for 3 0.34449 cm 

a 0.00025 cm 

a4 0.00021 cm= 

Dy; 331077 cm™! 

vo for 2v; +24 8589.28 cm= 

Vo for 3r1 8792.70 cm— 

Tz) 81.08 X 10-*° ¢ cm? 


Tc \ calculated assuming tetrahedral {148.67 x 10-* g cm? 
Cio{ angles \0.1882; cm™ 
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for the F-C-F angle. If one assumes a tetrahedral 
angle the C—F distance remains the only un- 
known, and one obtains 


r(C—F) =1.329X 10-8 cm. 


For other values of the F-C-F valence angle 
Fig. 4 gives the corresponding r(C—F) values. 
It is seen that for F-C-F angles larger than the 
tetrahedral angle the r(C —F) distance is smaller 
than the value given above. _ 

It is interesting to compare the C—F distance 
in fluoroform with that in other related mole- 
cules. This is done in Table IV where both 
spectroscopic and electron-diffraction data are 
summarized. While the electron-diffraction values 
are rather larger than the spectroscopic values 
in all cases where both are available, for either 
set the C—F distance in CH;F is appreciably 
larger than in the molecules containing more than 
one fluorine atom. Of course by a proper choice 
of the F-C-F angle in CHF; it would be possible 
to obtain from the above J," value a C—F dis- 
tance equal to that in CH;F. However, the 
angle would have to be as low as 99°, which 
seems highly improbable. 


120° A F-C-F 


An unambiguous determination of the C—F 
distance as well as the F-C-F angle will be 
possible once the moment of inertia J, of 
“heavy” fluoroform (CDF3;) has been deter- 
mined from its parallel bands. 

The C—F distance given above and the as- 
sumption of tetrahedral angles may be used to 
calculate the moment of inertia J¢"! about the 
symmetry axis as well as the rotational constant 
Cio}. The values obtained are included in Table 


TABLE IV. C-F distances in different molecules. 











C-F distance (in 1078 cm) 
from from electron 

spectrum diffraction F-C-F-angle 
CH3F 1.380* 1.42¢ 
CHF. 1.320> 1.36° 106°50’> 
CHF; 1.329 109°28’ (assumed) 
CF, 1.36¢ 109°28’ (assumed) 
CF;CH; 1.374 107°32’ 














* Derived from the + value given in a recent paper by Yates and 
Nielsen [Phys. Rev. 71, 349 (1947)] assuming an H-C-H angle of 
111°48’ (see G. Herzberg, J.R.S., p. 439). The large deviation from 
the value given by Herzberg is due to the large change of Yates and 
Nielsen's Bro} as compared to the earlier one of Bennet and Meyer 
[Phys. Rev. 32, 888 (1928)]. 

b From a paper presented by W. H. Shaffer and H. Stewart at the 
Symposium on Molecular Structure and Spectroscopy at Ohio State 
University, June 1947. 

¢L. O. Brockway, J. Phys. Chem. 41, 185 (1937). 

4 W. Shand and R. A. Spurr, quoted in H. Russell, Jr.. D. R. V. 
Golding, and D. M. Yost, J. Am. Chem. Soc. 66, 16 (1944). 
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TABLE V. Fundamental frequencies of gaseous CHF; after 











Price.* 

Assignment Infra-red band 
vi1°8 (a1) | 3035.6 (s.) 
vo&F (a) } 1209 (m.) 
v3(a)) 703.2 (m.) 
vslF (e) 1351.5 (s.) 
vs°H(e) 1152.4 (v.s.) 
ve(e) 509.4 


| cm! 








* Note added in proof: New measurements by Price in the long wave- 
length region seem to indicate that the two low frequency fundamentals 
v3 and ve have to be exchanged since the band at 509.4 cm~! has the 
appearance of a parallel band while the band at 703.2 cm™ has an 
indifferent appearance. 


111. However, it must be remembered that these 
values are not a direct consequence of the 
measured wave numbers as are Jz"), By, and 
the other constants of Table III. 


D. VIBRATIONAL STRUCTURE 


The Raman spectrum of fluoroform in the 
liquid phase has been studied by Glockler and 
Edgell.?. The infra-red spectrum of the gas 
in the region of the fundamentals has been 
investigated by Price.* There is an appreci- 
able shift of the fundamental frequencies in 
going from the gas to the liquid. Therefore, for 
the interpretation of the overtone and combina- 
tion bands found here in the gas, we use the 
fundamentals determined by Price as given in 
Table V. The designation is that used by 
Herzberg. There are three totally symmetric 
fundamentals (a,) and three degenerate funda- 
mentals (e). 

In addition to the two photographic bands 
whose fine structure has been analyzed (Section 
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C), a considerable number of fainter photo- 
graphic bands were previously found by Mackay’® 
and have been confirmed in the present work. 
They extend as far as 7000A. Furthermore, 
many strong bands were found in the region 
from 2.4y to 1.24 with the photoelectric infra-red 
spectrometer referred to in Section B. Figure 5 
reproduces one of the recordings obtained. The 
wave numbers of all the observed overtone and 
combination bands are given in Table VI. The 
wave numbers measured with the infra-red 
instrument (marked I. in the second column) 
are accurate to about +5 cm~'; those measured 
in the photographic region (marked P.I.) are 
accurate to less than 1 cm™ except in the case 
of very diffuse bands. 

The dispersion of the photoelectric spectrom- 
eter is not sufficient to distinguish parallel and 
perpendicular bands. But even with the high 
resolution possible in the photographic infra-red 
it is often difficult to decide whether a given 
weak band is a parallel or a perpendicular band. 
If the band has a very broad central branch 
without a head it is fairly certain that it is a 
perpendicular band. In one case of such a band, 
at 10154.6 cm™, the Q branches of the sub-bands 
forming this central branch were partially re- 
solved into a series which even showed an 
indication of the proper intensity alternation 
(strong, weak, weak, strong) so that its identi- 
fication as a perpendicular band is beyond doubt. 
The completely resolved parallel bands, 8589 and 
8793 cm (Fig. 2), have a sharp edge of their 
central branches since B’ is somewhat smaller 
than B’’. Other bands showing such heads, even 
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Fic. 5. Infra-red spectrum of CHF; obtained with the photoelectric infra-red spectrometer. 


8 W. C, Price, unpublished. 
°G. Herzberg, I.R.S., p. 316. 


0 A, E. Mackay, Master’s thesis, University of Saskatchewan, 1942. 


7G. Glockler and W. F. Edgell, J. Chem. Phys. 9, 224 (1941). 
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TABLE VI. Overtone and combination bands of CHF;. 

















Yvw 
Vyac, | Cwict 
observed | bateseis 
(cm~!) | Band type* Assignment? (cta7)) 
4340 v.s. I. Le vitva(E) 4387 
4540 v.w. I. 3vs+ve(A1+3E) 4580 
yor = ; | -+ys-+r0(A1 +E) 4702 
4770 v.w. I. vitve+ve(E) 4759 
4872 v.w. I. vitvetve(A1 +E) 4901 
5047 v.w. I. vitvs+va(E) 5087 
5186 m. I. vi+2v4—ve[Ai1+2E] 5185 
5330 a. f. vit+2v5(Ait+E) 5340 
5530 ad vitvat+yvs(Ait+E) 5539 
oo : one) \n+20u(As +E) see text 
5938 re & vitve+u+ve(Ai1+E)? 5911 
5978 v.s. I. (doublet) | 271(A1) see text 
6290 v.w. I. vitve+v3+ra(E) 6296 
6435 v.w. I. 2r1+0(E) 6492 
6664 w. I. 2v1+v3(A1) 6678 
6853 m. I vit2vs+r5(A1+2E) 6851 
6906 =. 3. vi+2v4+v2(A1+E) 6908 
7023 s. I. vi+3va(A1+E) 7050 
7100 a A | 2x1 +5(E) 7130 
7318 v.s. I. 2v1+v4(E) 7329 
8135 v.w. I. 2vuit+vs+va(E) 8029 
8270 m. I, 2v1+25(Ai +E) 8283 
8551.0 v.w. P.I. 2v1 +2v4(E) see text 
8589.30 m. P.I. || head 2v1 +2v4(A1) 
8586.1 w. P.I. || head 2r1+2r4+v6 — v6 see text 
8791.1 w. P.I. || head 3vi+ve — 6 
8792.75 m. P.I. || head 311(A1) see text 
9540 v.w. I. 3v1+¥3(A1) 9493 
9730.0 va PEL 21 +2v4+05(A1+2E) 9742 
9881.6 w. P.I.4 2v1+3v4(A1 +E) 9941. 
9938.3 ow, PS. Lt 3r1+5(E) 9945 
10154.6 w. P.Le | 3rni+va(E) 10144 
11109.2 v.w. P.L. line 3v1+25(A1 +E) 11098 
> 
113002 | vw. Pilling [}S-tretn (Ante) 11297 
11346.2 v.w. P.I. line 3v1+2v4+v6 —ve[2A1+4E] 
> a. 
iwa970 | vow. Pa. Lt (}S"t2n(an te) cee text 
11555.5 v.w. P.I. || head | 4v1-+-v3 —v3[A1 +E] 
11559.2 v.w. P.I. || head | 4v1-+-v6 —ve[A1+£] 
11563.4 w. P.I. || head 41(A1) see text 
1264s | vows Padcline  |}S"+3(4i +2) 12699 
eee | vw Ty [}4et teeta te) see text 











® The intensities are indicated by v.s., s., m., w., and v.w. meaning 
very strong, strong, medium, weak, and very weak, respectively. I. 
and P.I. indicate bands measured with the photoelectric infra-red 
spectrometer and with the 21-ft. grating in the photographic region, 
respectively. || and | refer to bands that have been identified as 
parallel or perpendicular bands. Bands marked with “line” have a 
line-like but diffuse central branch whose center was measured. They 
may be either parallel or perpendicular bands. Parallel bands marked 
head show a pronounced head of the central branch to which the 
wave number of the first column refers, 

b The species of the allowed components of the transition moment 
are given in brackets. For all but the difference bands they coincide 
with the species of the upper states (omitting the Az substates). Note 
added in proof: There will be a few changes of species if v3 and ve are 
interchanged; see footnote (a) of Table V. 

© This band was also measured by Price with 10-cm path, at 4310 
cm~!, The identification of the band as | band is due to him. 

4 This band has a very peculiar structure. There is a strong violet- 
shaded head to which the wave number given in column 1 refers; 
in addition there is a diffuse band shaded apparently to longer wave- 
lengths. 

e This band shows more clearly than any other photographic band 
the expected structure of a perpendicular band. The series of 0 branches 
is partially resolved, and the appropriate intensity alternation (strong, 
weak, weak, strong) is just recognizable. 
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if the P and R branches are too weak to be 
observed, must also be identified as parallel 
bands. However, there is a considerable number 
of bands whose central branches (frequently the 
only feature observed) are like diffuse lines. 
They may either belong to parallel bands for 
which B’=B” or to perpendicular bands for 
which the spacing of the sub-bands, which is 
27 B—C(1—¢)], is small because the Coriolis 
interaction constant ¢ happens to be close to 
— 0.84. 

It was expected that the series of bands 
v1, 21, 31, 4v1, --+ (that is the overtones of the 
C—H vibration) would be the most prominent 
feature of the near infra-red spectrum of fluoro- 
form, as it is for chloroform. Actually there are 
in each case two neighboring parallel bands of 
similar intensity rather than one. The two bands 
analyzed in Section C (see Fig. 2) form one such 
pair, namely, the one corresponding to 3y, It is 
readily seen that these pairs of bands can be 
interpreted as v1y; and (v;—1)v1 +24, since 2v4 
has a magnitude similar to », particularly to 
the latter’s higher vibrational quanta, and since 
one substate of 2v, has the same species as 1, 
that is A;. The separation of the two components 
of a pair is larger than 200 cm~ and therefore 
Fermi resonance will affect the positions of the 
levels only slightly, even though it may affect 
the intensities to an appreciable extent. The 
shortward components of each pair which repre- 
sent v1, 271, 31, 4v1, are 3035.6, 5978, 8792.7, 
and 11563.4, respectively. The separations of 
successive vibrational levels of the CH vibration 
are, therefore, 3035.6, 2942.4, 2824.;, and 2770.7 
cm. The levels can be represented approxi- 
mately by 

3088.00; — 52.401’, 


although the deviation for v;=4 is much larger 
than compatible with the accuracy of the meas- 
urements, probably because of a non-negligible 
cubic term. The anharmonicity x:,;=52.4 is of 
the right order for a CH vibration (compare for 
example those of HCN and C;H.,). 

The separations of v1; +24, 2v1+2v4, --- from 
v1, 2”, ‘, respectively, are 2663.4, 2611.3, 
2554.5, and 2439.4 cm~'. The accelerated de- 





f Accompanied by a diffuse fine structure on both long and short 
wave-length sides. 

a were calculated with old values, 700 and 514 cm™, for » 
and ve. 
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crease again points to a cubic term in the energy 
expression. Neglecting this, the anharmonic term 
X14 comes out to be approximately 28 cm™. 
However, this anharmonic coefficient does not 
represent very well bands of the type 171+ 4 
and v17;+3v4. This is probably because the 
observed bands v1»4+ 2, are parallel bands, hav- 
ing the A; component level as upper state while 
the bands vy,+v, are perpendicular bands, the 
upper state having the species E only. For 
71¥1+3v, the situation is more complicated. 

The vibrational data here obtained are not 
sufficient to determine all the anharmonicity 
coefficients x;,. Since, in addition, the effect of 
cubic terms and of Fermi resonance is difficult to 
take into account accurately at the present 
stage, we give in the last column of Table VI, for 
all bands other than those just discussed, the 
wave numbers computed by simply adding to 
the observed wave numbers of the bands of type 
vv; and v\»1;+2v, the fundamental or sum of 
fundamentals involved according to the assign- 
ments given." On account of the effect of an- 
harmonicity and the other reasons mentioned 
above, one should not expect a very close agree- 
ment of calculated and observed wave numbers. 


'' See footnote (g) of Table VI. 


The only reason for listing the calculated values 
is as support for the assignments given in the 
third column. 

The central branches of the strong photo- 
graphic bands are accompanied by weaker com- 
panions which must clearly be interpreted as 
difference bands of the type v-+v;—y; in which a 
low frequency fundamental is excited in both 
upper and lower state of the main band. Their 
positions relative to the main band are deter- 
mined by the anharmonicity constants and can- 
not be predicted. 


We are greatly indebted to Dr. A. F. Benning 
of the Jackson Laboratory of E. |. du Pont de 
Nemours and Company for putting at our dis- 
posal the fluoroform used in this investigation. 
Grateful thanks are due Dr. G. P. Kuiper for his 
cooperation in taking the infra-red spectra of 
fluoroform with his new photoelectric infra-red 
spectrometer. We also wish to thank Dr. W. C. 
Price for letting us have his data about the 
fundamentals of CHF; in advance of publication. 
Finally, we would like to express our appreciation 
to Professor R. W. Wood for supplying us with 
the remarkable grating with which the high 
resolution work was carried out. 
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The reactions of ethylene oxide with Hg(?P;) atoms have been investigated at 25°C at 
pressures between 10- and 600-mm Hg, and at 200°C and 300°C at 100-mm Hg. The principal 
products are He, CO, CH;CHO and a polymer; smaller amounts of CH,, C2 hydrocarbons and 
formaldehyde are also formed. The nature of the products and the very small effect of pressure 
on the reaction are best explained by a mechanism based principally on the formation of the 
activated ethylene oxide molecule, with a minor portion of the reaction proceeding via a C—H 


split. 


INTRODUCTION 


HE reactions between Hg(?P;) atoms and 
the lower hydrocarbons have been studied 


* Contribution No. 1602 from the National Research 
Laboratories, Ottawa, Canada. 


in some detail.! It appears that, in general, the 
unsaturated hydrocarbons react principally by 
the transfer of energy from the excited mercury 


1 See E. W. R. Steacie, Atomic and Free Radical Reactions 
(Reinhold Publishing Corporation, New York, 1946). 
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atom to the hydrocarbon molecule, leading to 
the formation of an excited hydrocarbon from 
which the products arise; on the other hand, the 
saturated hydrocarbons appear to react by a 
direct C—H split, which may or may not involve 
the formation of HgH, leading to the formation 
of the corresponding alkyl radical and an H 
atom. It has recently been shown that diethyl- 
ether? probably reacts with Hg(*P,) atoms in a 
manner similar to the paraffins, whereas acetone*® 
appears to behave like the olefins, and to react 
by the transfer of energy. The mercury photo- 
sensitized reactions of ethylene oxide are of 
interest in connection with the photo-chemical 
reactions of the simple oxygen containing organic 
compounds, since the oxygen atom may form 
part of a ring structure.* - 


EXPERIMENTAL 


The ethylene oxide used was obtained from 
the Matheson Company, East Rutherford, New 
Jersey. It was used after repeated melting and 
freezing in vacuum followed by a simple bulb to 
bulb distillation. No impurities were detectable 
by the analytical methods used. The absorption 
spectrum was checked at atmospheric pressure 
on a Beckmann quartz spectrophotometer (10- 
cm path), and it was found to be completely 
transparent to 2200A. 

The reaction was investigated in a circulatory 
system in which all of the ethylene oxide was 
present in the gas phase, no attempt being made 
to remove the higher boiling products of the 
reaction by means of cold traps since very small 
conversions were used. The gases were circulated 
by means of a pump of the type described by 
Puddington® over warm mercury, through a 
mercury desaturator and a cylindrical quartz 
vessel with plane ends in which the reactants 
were illuminated with the unfiltered light from 
a low pressure mercury lamp with neon (3 mm) 
as carrier gas. Since the incident light traveled 
through a distance of about 12 inches in air, it 


2B. de B. Darwent, E. W. R. Steacie, and A. F. van 
Winckel, a Chem. Phys. 14, 551 (1946). 

3E. W. R. Steacie and B. de B. Darwent, unpublished. 

‘For opinions on the structure of ethylene oxide see, 
however, P. by Zimakov, J. Phys. Chem. (U.S.S.R.) 20, 
133 (1946); A. D. Walsh, Nature 159, 165, 712 (1947): 
R. Robinson, ibid. 159, 400 (1947), 160, 162 (1947). 

5]. E. Puddington, Ind. Eng. Chem. (Anal. Ed.) 17, 
592 (1945). 
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is reasonably certain that only traces of \1849 
could have reached the reactor. 

The products of the reaction were separated 
by pumping off a trap kept successively at 
—183°, —126°, and —80°C. The products 
volatile at —183°C (Hz, CO, and CHy,) and at 
—126°C (C. hydrocarbons) were analyzed by 
conventional methods. The residual ethylene 
oxide (volatile at —80°C) was estimated by 
P—V relations and by titration as recommended 
by Lubatti.*® Specific tests showed formaldehyde? 
to be absent from the Cz (—183° to —126°C) 
fraction but to be present in the “‘ethylene oxide’”’ 
fraction. Acetaldehyde was also found® in the 
“ethylene oxide” fraction in which it was 
estimated by the sodium _bisulfite—iodine— 
sodium thiosulfate titration? after correction 
for formaldehyde which was estimated colori- 
metrically.!° A small blank correction was neces- 
sary when acetaldehyde was estimated after the 
titration for ethylene oxide. The correction was 
negligible and the reproducibility better when 
the ethylene oxide titration was omitted; in 
these cases the residual ethylene oxide was deter- 
mined by P—YV relation by correcting for the 
total aldehyde present. No detectable pressure 
developed on allowing the trap to warm up, 
after removal of ethylene oxide, from —80° to 
room temperature. 

A non-volatile white solid was visible on the 
walls of the reaction vessel at the end of each 
experiment at room temperature. It was also 
found that under these conditions a large propor- 
tion of the ethylene oxide reacted could not be 
accounted for by the gaseous products and alde- 
hydes. Accordingly, experiments were conducted 
in a single-pass system to obtain enough of the 
heavy products for identification. In these experi- 
ments a much more intense light source was used 
and the flow rate of the ethylene oxide adjusted to 
permit 10-20 percent reaction per pass. No white 
solid was deposited on the walls, but droplets of 
a relatively high boiling liquid appeared early 


in these experiments. The products obtained 


6 QO. F. Lubatti, J. Soc. Chem. Ind. 54, 424 (1935). 
( 910) B. Schryver, Proc. Roy. Soc. (London) 82B, 226 
1910). 

8 A. Ley, J. pharm. chim. (6) 22, 107 (1897). 

97. M. Kolkoff and N. H. Furman, Volumetric Analysis 
_— Wiley and Sons, Inc., New York, 1929), Vol. II, p- 


TD. Matukawa, J. Biochem. (Tokyo) 30, 386 (1939). 
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from one experiment separated on cooling into 
two layers, the upper of which boiled at about 
165°C and the lower at about 99°C. Fractional 
distillation of the products of further experi- 
ments gave three fractions, the identifications of 
which were unsuccessful, but tests showed them 
to be aldehydic in nature and showed the absence 
of such plausible compounds as dioxanes, dioxo- 
lanes, 1,3-diepoxybutane, and alcohols. The 
white polymer produced in the low intensity cir- 
culating runs was found to be water soluble and 
aldehydic in nature; it was easily removable 
from the walls by warming or by flaming in an 
oxygen atmosphere, the latter being the pro- 
cedure adopted for cleaning the cell before each 
run. It is worthy of note that this white polymer 
was not formed on the incident face, but rather 
on the cylindrical side walls of the cell, and 
started about 1 mm from the incident face. 

Some experiments were carried out at higher 
temperatures. In these runs the quartz reactor 
was enclosed in a furnace fitted with a quartz 
window. All other conditions were unaltered. No 
visible polymer was formed in these experiments, 
and the material balance was satisfactory at 
200°C and higher temperatures. 

The light intensity was measured by the rate 
of hydrolysis of monochloracetic acid (0.5) 
in the cell. The corrections as given by Smith, 
Leighton, and Leighton" for the dark reaction 
were applied. 


RESULTS 


Details of the experimental conditions and 
results on which subsequent discussion will be 
based are given in Table I. From these results it 
is evident that the rates of production and com- 
position of the gaseous products at any given 
pressure are independent of time. The data for 
the amount of ethylene oxide reacted are much 
more variable, being small differences of two 
large quantities, but the rate here may also be 
taken as independent of time within the limits 
of accuracy of the determination. The rate of 
production of aldehydes is not independent of 
time (Fig. 1); however, it seems likely that the 
initial rate of production of aldehydes is approxi- 





"R, N. Smith, P. A. Leighton, and W. G. Leighton, 
J. Am. Chem. Soc. 61, 2299 (1939). 


mately independent of the initial pressure of 
ethylene oxide. By basing the discussion on the 
initial rates the complications connected with 


TABLE I, The mercury photosensitized decomposition of 
ethylene oxide. 








Experimental conditions and results 


Volume of system—726 cc, A2537 absorbed—9.82 10-5 einstein/hr. 
Volume of reaction vessel—230 cc, Hg pressure, 1.2 X10-? mm, 
Circulation rate—400 cc/min. at prevailing pressure. 


I. Experiments at room temperature—25+2°C 





Ethylene 
oxide 
Dura- reacted Products ce at 25°C 
Expt. tion, cc Total 
No. hours 25°C He co CH«e C2 ald. 
A. Initial pressure of ethylene oxide =10+1 mm 
29 1.0 0.5 0.27 0.46 0.09 0.16 -— 
28 2.0 1.8 0.48; 0.92 0.19 0.28 -— 
33 5.1 2.4 0.98; 1.91 0.39 0.52 — 
B. Initial pressure of ethylene oxide =40 mm 
44 0.25 — - - -—- = 0.05 
45 0.5 — —- —- —- =— 0.16 
9 0.5 0.5 0.19 0.26 0.05 0.08 0.18 
3 1.0 Lf 0.38 0.58 0.10 -0.15 0.16 
46 1.0 — — — — — 0.16 
1 2.0 2.5 0.72 1.06 0.20 0.33 0.10 
2 4.0 3.9 1.36 2.04 0.37 0.55 0.07 
11 4.0 -- —_-_ =—- _- = 0.09 
C. Initial pressure of ethylene oxide = 100 mm 
43 0.25 — _-_ - —_—- — 0.10 
32 0.25 0.6 - - —- =— 
42 0.5 — —_— —_-— -_- — 0.15 
25 1.0 2.2 0.39 0.52 0.09 0.12 = 
41 2.0 — _-_ - -—- — 0.47 
23 4.0 5.2 LS2 2 0.32 0.46 — 
40 4.3 — _-_ - —_-_ — 0.61 
24 8.0 8.88 2.78 3.57 0.59 0.66 os 
D. Initial pressure of ethylene oxide =200 mm 
14 0.17 — _-_ - -—- — 0.15 
10(A) 0.5 — 0.21 0.24 0.05 0.08 -- 
12(A) 0.5 — —-_ - —_- — 0.27 
8 1.0 1.7 0.42 048 0.09 0.11 0.20 
11(A) 1.0 — —- —- —- — 0.42 
4 2.0 2.9 0.88 1.04 0.16 0.20 0.48 
6 4.0 6.3 1.72 1.94 0.31 0.31 0.81 
13 6.0 — —- -—- —- — 0.91 
5 8.0 10.4 3.12 3.66 0.57 0.55 0.75 
E. Initial pressure of ethylene oxide =600 mm 
18 1.0 2.2 0.45 0.43 0.09 0.10 — 
21 1.0 2.6 —_-_ —- —- — — 
15 2.0 _— 0.89 0.91 0.14 0.18 0.42 
20 3.0 5.1 —- -—- -—- = — 
16 4.0 — 1.83 1.78 0.28 0.24 — 
16(A) 6.0 9.7 2.57 2.42 0.37 0.42 0.99 
22 10.0 12.9 4.48 3.94 0.67 0.64 — 


II. Experiments at 200°C and an ethylene oxide pressure 
of 100 mm 


36 1.0 2.5 0.70 1.29 0.30 0.37 1.66 

37 3.0 4.5 2.03 4.09 0.95 1.19 1.74 

III. Experiments at 300°C and an ethylene oxide pressure 
of 100 mm 

38 1.0 3.1 0.85 2.03 0.71 0.53 1.57 

38(A) 1.0 (Blank —lamp off —total non- 1.41 


condensible gas = 0.043 cc) 
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DURATION OF EXPERIMENT (HRS) 





Fic. 1. The production of aldehydes at 25°C; initial 
pressure of ethylene oxide as follows: squares—600 mm; 
circles—200 mm; triangles—100 mm; half-circles—40 mm. 
(In experiments at 40, 200, and 600 mm the ethylene oxide 
was estimated by titration before the aldehydes were 
determined; the aldehyde content was less accurate than 
in the experiments at 100 mm in which no titration for 
ethylene oxide was attempted. 





the disappearance of aldehydes in the later stages 
of the reaction may be eliminated. 

The effect of pressure on the initial rates is 
shown in Table II. The data given here show 
that the gaseous and aldehydic products do not 
account for all the material reacted except at the 
lowest pressure. It has been mentioned previously 
that a white solid polymer was formed during the 
experiments. The effect of pressure on the rate 
of polymerization has been calculated (Table 
Ill) from material balances, and the results 
show the increase in the extent of the polymeriza- 
tion reaction with pressure. The composition of 
the polymer (ratios H/C and C/O) is given in 
columns 11 and 12 of Table III. The H/C ratio 
is constant, within experimental error, at 2.0 
and is approximately independent of pressure, 
whereas the C/O ratio decreases, indicating in- 


TABLE II. Effect of pressure on the mercury photosensi- 
tized decomposition of ethylene oxide at 25°C. 








A. Experiments at room temperature 
Initial 
pressure Rates—cc/hr. (25°) 





mm —(CH2)20 He co CHa Ce ald. 
10 0.56 0.21 0.42 0.082 0.118 0.45 
40 1.02 035 33 .094 145 45 

100 1.30 .36 475 .076 .093 45 

200 1.44 411 47 .074 .081 45 

600 1.50 445 42 .070 .067 45 
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creasing oxygenation, with increasing pressure. 
The above results and the fact that two immis- 
cible liquids were found in the single-pass ex- 
periments lends plausibility to the suggestion 
that the total polymer consisted of two types of 
substance ((CH2)2O), and (CH,),, in which X 
lies between 1.2 and 2.1 with an average value 
of 1.8, as shown in Table IV. 

The effect of pressure is shown in Fig. 2 in 
which the rates given in Tables II and IV are 
plotted against the reciprocal of the pressure. 
The relationship between pressure and the rates 
shown in Fig. 2 for the experiments at room 
temperature is as follows: 

3.5 


d 
dt p 





d 4.5 
—(CO) =0.425+—, 
dt p 
d 25.0 
—[((CH2).0),,] = 0.660 -—, 
dt p 
d 3.1 
—(C;) =0.060+-—, 
dt p 
d 1.0 
—(CH,) =0.064+—, 
dt p 


where the rates are in cc (25°C) per hour and 
p=pressure (mm Hg) of ethylene oxide. The 
lowest pressure investigated (10 mm) has not 
been included in Fig. 2 because the results indi- 
cate that quenching was incomplete below about 
40 mm. 


DISCUSSION 


The following processes have been shown to 
occur in the mercury photosensitized reactions of 
ethylene oxide at room temperature: 


(a) Isomerisation to acetaldehyde: the initial 
rate of production of acetaldehyde being 
independent of the ethylene oxide pressure. 

(b) Decomposition to Hz, CO, CH4, C2 hydro- 
carbons and formaldehyde. The rates of 
production of the above substances (ex- 
cept formaldehyde) are linear functions of 
the reciprocal pressure such that the rates 
are finite even as p—, the pressure de- 
pendence in all cases being small. 
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TABLE III. Effect of pressure on the mercury photosensitized polymerisation of ethylene oxide at 25°C. 

















Reaction rates Gases and aldehydes* Polymer (by diff.) Comp. of 

Pressure cc atom/hr. (25°C) ce atom /hr. (25°C) cc atom/hr. (25°C) polymer 
mm - H Oo te H Oo * H Oo H/C C/O 
10 1.12 2.24 0.56 1.52 2.92 0.87 —0.40 0.00 —0.31 ms a 
40 2.04 4.08 1.02 1.70 3.38 0.98 0.34 0.70 0.04 ook 8.5 
100 2.60 5.20 1.30 1.$2 3.04 0.92 1.08 2.16 0.38 2.0 2.8 
200 2.88 5.76 1.44 1.48 3.08 0.92 1.40 2.68 0.52 1.9 om 
600 3.00 6.00 1.50 1.40 3.09 0.87 1.60 2.91 0.63 1.8 y 








+ The average composition of the aldehydes was 75 percent CH3;-CHO, 25 percent HCHO. 


(c) Polymerization to yield substances of the ethylene oxide, would lead to kinetics such that 
types ((CH2).0), and (CHe),—the rate of _ the rate of formation of acetaldehyde would be 
production of the former being quite inversely proportional to the concentration of 
strongly pressure dependent, whereas the ethylene oxide. 
latter is approximately independent of The effect of pressure on the decomposition to 
pressure. hydrogen and carbon monoxide is also in agree- 

ment with the supposition that an active mole- 

cule plays the predominant role in the initial step 


lack iil ae Se saneicer as a of this reaction. If the decomposition were 
I » pressure in > 2 is ‘ 
dahl ee tie ail maypinlas nage largely due to a C—H split followed by ‘‘atomic 


reaction, indicate that an active ethylene oxide nae 
Se a Ree a al ; é cracking,” the type of process frequently pos- 
> > 1 ) 1 c 7 c- . . . 
ee ee er ee ee tulated to account for C—C splitting in the 
Hg(?So) +hv2537Hg(*P1), (0) photosensitized reactions of the paraffin hydro- 


Hg(P1) +(CH:2),0-C.H,O*+Hg('So), (1) carbons, it would be expected that, since H 
atoms react with ethylene oxide™ at about the 


The occurrence of an isomerisation reaction, 
by which acetaldehyde is formed as a major 


from which acetaldehyde is produced by either 



































a first or second order process: 0.8 
C,.H,O*—CH;CHO, (2) 
0.6 
- o Be 
C,H,O*+(CH2)2.0—>CH;CHO+(CHz2),0. (2) : oe on oe - 
. rA_ o> 
The above reactions (2) and (2!) are energetically « ” dt a 
possible since the active molecule may possess ‘j - | me 
as much energy as 112 kcal. per mole, and thermal : NI] POLYMER 
investigations” indicate that only 52 kcal. per E ‘il | yx CeCH2d20l4 


mole are necessary for isomerization. Any process 
by which acetaldehyde could be formed by 
atomic and/or radical recombination, providing 
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TaBLE IV. Effect of pressure on the composition of 5 “a ae J CHa 
the polymer. a = i 
eae es . a o ” 
— ——— — —— 7 
Pressure ce ((CH2)20)n ec (CHz)n El 0.04 
mm per hour per hour X « 
10 es rae ria) 0.00! 
a oa me 7 “- 5 Tr ry a a 
ow oS be x 103 
200 0.52 0.36 Pe habe 
600 0.63 0.34 Lz 


: Fic. 2. The effect of pressure on the decomposition and 
said cuties polymerization reactions at 25°C. 

2 W. W. Heckert and E. Mack, J. Am. Chem. Soc. 51, 
2706 (1929). - 18 W. R. Trost and E. W. R. Steacie, unpublished. 
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same rate at room temperature as they do with 
the paraffins, the rate of decomposition would 
approach zero as the concentration of ethylene 
oxide approaches infinity. Since such a rela- 
tionship has not been found, it must be con- 
cluded that a C—H split cannot play an im- 
portant role in the mechanism. The decomposi- 
tion of the excited molecule may be first or 
second order: 


C.H,O*—H2+CO+ CHa, (3) 


C.H,O*F+ (CHg) 200O—-H.+CO+CH, 
+(CH2).0. 3") 


An alternative process by which the decomposi- 
tion proceeds by virtue of acetaldehyde in the 
stable or activated state, such as was found for 
the thermal and H-atom reactions" of ethylene 
oxide, does not agree with the results of the 
present investigation. That type of process 
would lead either to an induction period or to a 
pressure dependence such that the rate of de- 
composition would approach zero at higher 
pressures. 

If the excited molecule acquires almost all of 
the energy of the Hg(*P1) atom (112 kcal. per 
mole), it will be possible for reactions (2!) and 
(3') to occur in a single step leading to the 
isomerization of one molecule and to the decom. 
position of the other in the collision between one 
excited and one ground-state ethylene oxide 
molecule. 

The results obtained do not allow of a choice 
between the first- and second-order processes 
outlined above. The limiting high pressure rates, 
obtained by extrapolation in Fig. 2, give the 
following quantum yields: 


@ Decomp. = 0.18, 
@ Isom. = 0.12, 

@ Polym. =0.36, 
Total =0.66. 
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Hence only 0.30 mole of ethylene oxide isomerizes 
or decomposes per einstein absorbed. This inef- 
ficiency may be due to either a first- or second- 
order process: 


C.H,O*—(CH2).0+hr (4) 
C.H,O*+ (CH2)20—2(CH2).0. (4") 


Since no investigation of fluorescence was made, 
it is not possible to decide between the alter- 
natives since either a short-lived excited mole- 
cule with fluorescence or a long-lived one which 
on collision always either reacts or is deactivated 
will explain the results. 

The very small effect of pressure on the rate 
of production of the minor products (CH4, C2 
hydrocarbons and formaldehyde) indicates that 
they also arise in some way from the decom- 
position of the excited molecule. 

The above mechanism accounts for the major 
characteristics of the isomerization and decom- 
position reactions but not for the pressure de- 
pendence, small though it is, of the decom- 
position. Also the fact that polymerization is 
important at high pressures, and that more 
hydrogen than carbon monoxide is formed at 
high pressures show that a small portion of the 
reaction probably proceeds by a free radical 
mechanism. The following type of process 


Hg(®P1) + (CH2)20—-C.H;,0+H+Hg('Sp), (1") 


is well known in photosensitized reactions and is 
probably responsible for this aspect of the reac- 
tion. The nature of the polymer is consistent with 
the initiation of polymerization by the C.H;O 
radical, and the increase in ‘‘atomic cracking”’ 
with the decreasing pressure explains the small 
pressure dependence, providing that the major 
portion of the initial act is reaction (1) rather 
than (1'). The effect of temperature can be 
explained by the decomposition of the C,.H;O 
radical leading to increasing CO and decreasing 
polymer at high temperatures. 
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The Solubility of Water in Benzene 
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A new method for determination of the solubility of water in hydrocarbons is based on the 
use of radioactive hydrogen isotope, tritium, as tracer for the dissolved water. This method 
is applicable to the determination of the solubility of water in compounds in which it is 


sparingly soluble. 





LTHOUGH numerous methods have been 
used or proposed for the determination of 
water solubility in hydrocarbons, the number of 
hydrocarbons studied has been singularly limited. 
Solubility of water has been measured only in 
benzene,!~> toluene,” > cyclohexane,” various gas- 
olines,’**7 and oils.5~!° The synthetic or sealed- 
tube method employed by Groschuff,! and by 
Tarassenkow and Poloshinzewa? consisted in 
observing the temperature at which a measured 
amount of water is visually soluble in a measured 
weight of hydrocarbon. This is done by super- 
heating the mixture sufficiently to bring all the 
water into solution and then cooling until a fine 
mist is formed; from the temperature of appear- 
ance and disappearance of this fine mist, solu- 
bility data may be obtained. The saturated 
solutions were, in another method, volatilized 
and absorbed in a dehydrating agent which does 
not absorb or react with organic solvent (calcium 
chloride).*® The water dissolved in the hydro- 
carbons has also been determined by converting 
it to hydrogen through addition of calcium 


* Present address: Allied Chemical and Dye Central 
Research Laboratory, Morristown, New Jersey. 

' Groschuff, Zeits. f. Elektrochemie 17, 348 (1911). 

2N. D. Tarassenkow and E. N. Poloshinzewa, Zhur. 
obshchei Khim., Khim. Ser. 1, 71 (1931); Chem. Abs. 25, 
4762 (1931); Chem. Abs. 26, 2363 (1932); Ber. d. d. Chem. 
Ges. 65B, 184 (1932). 

3A. E. Hill, J. Am. Chem. Soc. 45, 1143 (1923). 

4C. W. Clifford, Ind. Eng. Chem. 13, 631 (1921). 

5C. K. Rosenbaum and J. H. Walton, J. Am. Chem. 
Soc. 52, 3568 (1930). 

6 Seidell, Solubilities of Inorganic and Metal Organic 
Compounds (D. Van Nostrand Company, New York, 1940), 
3rd edition, Vol. I, p. 587. 

7E. W. Aldrich, Ind. Eng. Chem., Anal. Ed., 3, 348 
(1931). 

8R. N. Evans and J. E. Davenport, Ind. Eng. Chem., 
Anal. Ed., 14, 732 (1942). 

®G. Wernimant and F. J. Hopkinson, Ind. Eng. Chem., 
Anal. Ed., 15, 272 (1943). 

10M. M. Acker and H. A. Frediani, Ind. Eng. Chem., 
Anal. Ed., 17, 793 (1945). 
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hydride® or a sodium-potassium alloy’ to the 
solution. The method used by Hill* is based upon 
the property of silver perchlorate of being 
exceptionally soluble in water and somewhat 
soluble in various organic solvents; the amount 
of water present in the hydrocarbon is estimated 
from the solubility of silver perchlorate in the 
solution, reference curves having been previously 
obtained using synthetic mixtures of water and 
hydrocarbon. The Karl Fischer method" has 
also been applied for the measurement of water 
solubility in oils.8-! 8 
The most serious disadvantage of all these 

methods is that they require large amounts of 
hydrocarbons if any accuracy is to be obtained. 
There was need for a method, accurate as well 
as general, which would require small amounts of 
solvent and therefore permit the study of sol- 
vents only available in the pure state in small 
quantities. The recent availability of water con- 
taining radioactive hydrogen rendered possible 
the development of a new method for the 
measurement of water solubility based on satu- 
rating the solvent with radioactive water and 
measuring the radioactivity of the solution. For 
radioactive hydrogen or tritium, having a half- 
life of 3148 years,“ there is no appreciable 
change in the radioactivity of the water during 
the course of the solubility measurement. The 
radioactive water used in this study was prepared 
in the Berkeley cyclotron by bombardment of 
heavy water with high speed deuterons; it was 
therefore necessary to correct the experimental . 


11 Brouckere and Gillet, Bull. Soc. Chem. Belg. 44, 473 
(1935). 
2D, M. Smith, W. N. D. Bryant, and J. Mitchell, Jr., 
J. Am. Chem. Soc. 61, 240 (1939). 
13 K, Fischer, Angew. Chemie 48, 394 (1935). 
( 4 R, D. O’Neal and H. Goldhaber, Phys. Rev. 58, 574 
1940). 

























results to obtain the solubility data of ordinary 
water in the solvents. Since the concentration of 
tritium in the water used was only of the order 
of 10-” mole tritium per mole of heavy water, 
the radioactive isotope is essentially all in mole- 
cules of TDO. The isotope effect due to tritium 
is correspondingly reduced. The mode of calcu- 
lation from this isotopic species to H2O solubility 
will be detailed later. 


EXPERIMENTAL METHOD AND TECHNIQUE 


The experimental technique consists in a series 
of successive operations: (a) the saturation of 
the hydrocarbon sample with radioactive water 
at the desired temperature by bubbling air, 
saturated with radioactive water vapor, through 
the hydrocarbon sample until saturation is 
attained ; (b) the taking of an aliquot portion of 
the saturated solution; (c) the measurement of 
the radioactivity present in the sample. 


APPARATUS AND OPERATING PROCEDURE 


1. Saturation and Sampling 


A unidirectional flow of air, saturated with 
radioattive water vapor and hydrocarbon vapor, 
is effected in the closed circulatory system 
represented in Fig. 1 by the combined action of 
a mercury-filled reciprocating pump and the two 
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mercury valves, V and V’. The removable 
chamber A contains the liquid hydrocarbon 
solvent; its temperature is kept constant during 
saturation by means of the water-filled Dewar 
flask, N. The circulating air has to go through 
trap C before entering the saturator A; the 
temperature of trap C is held at about 0.1°C 
below the temperature of saturator A by means 
of the water-filled Dewar flask, M. This temper- 
ature differential is important and has to be 
carefully maintained since it prevents the con- 
densation of water on the inlet tube walls of the 
saturator and on the hydrocarbon itself. In order 
also to avoid condensation of liquid in parts of 
the apparatus other than the trap C the whole 
circulatory system is kept at a temperature 
higher than that of the Dewar flasks by enclosure 
in an electrically heated air thermostat. In 
addition, the vertical tube, L, is electrically 
heated to prevent condensation of the vapors 
over the cool mercury. 

After complete evacuation (10-° mm Hg) and 
flaming of the apparatus, dry air is introduced, 
trap C being at liquid air temperature. When 
the system is at atmospheric pressure, it is 
isolated and the circulating pump is started. A 
small amount (about 0.1 cc) of radioactive water 
is then distilled from its storage chamber into 
the apparatus where it is trapped in cylinder C. 
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The circulation is stopped and about five cubic 
centimeters of liquid hydrocarbon solvent are 
added to the removable saturator A. Tempera- 
tures of the trap C and saturator A are adjusted 
to the desired level by means of the water-filled 
Dewar flasks, stopcock 2 being opened to the 
atmosphere to permit the elimination of the air 
in excess. When the conditions of temperature 
are attained, stopcock 2 is closed and the air is 
circulated. Results indicate that about four 
hours of circulation at the rate of approximately 
300 cc per minute will effect complete saturation. 

When saturation has been completed and 
circulation stopped, part of the saturated sample 
is drawn up into the inlet tube S by opening 
stopcock 2 to the vacuum system. Saturator A 
is then detached from the apparatus at Y. The 
outside of the tube S is carefully wiped with 
absorbent cotton moistened with~ alcohol to 
remove all traces of radioactive water and then 
dried. The calibrated sample tube A’ (see Fig. 1) 
is attached at the joint Y and by turning stop- 
cock 2 to the atmosphere the saturated sample 
is delivered into tube A’. The sample tube A’ is 
then removed and transferred to the other 
section of the apparatus where its radioactivity 
can be measured. 


2. Measurement of the Radioactivity of the 
Saturated Sample 


Tritium emits soft beta-radiations'® whose 
intensity is best measured by direct introduction 
of the radioactive sample into the Geiger-Miiller 
counter. This was most conveniently done by 
separating the radioactive water in solution from 
the hydrocarbon solvent through absorption on 
calcium oxide and by subsequently transferring 
the tritium in the counter through equilibration 
with ethanol vapor. The analytical apparatus is 
schematized in Fig. 2. 

When the volume of the saturated hydro- 
carbon sample, taken in the calibrated sample 
tube A’, has been measured at an observed 
temperature, the sample tube A’ is connected to 
the analytical apparatus by means of the ground- 
glass joint at W. The sample is solidified by 
cooling A’ with liquid air and the apparatus is 
evacuated in its entirety to remove the non- 
condensable gases. The sample is then transferred 
into the U-tube E which contains freshly de- 
hydrated calcium oxide and is cooled with 
liquid air while the temperature of A’ is warmed 
slowly up to room temperature. 

When the transfer is visually completed the 


18 R, D. O'Neal, Phys. Rev. 60, 359 (1941). 
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U-tube is isolated from the rest of the apparatus 
and brought back to room temperature; water 
is allowed to react with calcium oxide. The 
radioactive water left adsorbed on the walls of 
the sample tube is transferred into the Geiger 
counter for counting through equilibration with 
ethanol vapor and introduction of the latter into 
the counter. This is most conveniently done by 
evaporating into the sample tube section small 
amounts of ethanol vapor and expanding them 
into the counter until the total pressure reaches 
5 mm of mercury. Small amounts of argon gas 
are then similarly introduced into the sample 
tube section and expanded into the counter 
until the total pressure is 2 cm of mercury. The 
radioactivity of the mixture is then measured. 
The same operation is repeated until the radio- 
activity becomes negligible. 

After absorption of the water by the calcium 
oxide the hydrocarbon is removed by vacuum 
distillation. The first distillation fractions are 
introduced into the Geiger counter to ascertain 
the amount of possibly unreacted water they 
may carry. When the distillation fractions are 
exempt from radioactivity the hydrocarbons are 
condensed again into the sample tube A’. The 
bulk of the water in the hydrocarbon sample is 
now present as calcium hydroxide from which 
tritium is extracted for counting, again through 
exchange reaction with ethanol vapor. This is 
effected by introducing into the U-tube a quan- 
tity of ethanol such that when expanded into 
the evacuated Geiger counter it has a pressure 
of 5 mm; the U-tube is heated to 250°C and 
kept at that temperature for 30 minutes. The 
radioactive ethanol vapors are then expanded 
into the Geiger counter; small amounts of argon 
gas are successively introduced into the U-tube 
and expanded into the Geiger counter until the 
total pressure in the latter reaches 20 mm. The 
radioactivity of the ethanol vapors is determined. 
Since the equilibrium constant for the isotopic 
hydrogen exchange between calcium hydroxide 
and ethanol is of the order of unity at 250°C, 
the process of extraction of tritium from the 
hydroxide has to be repeated until the radio- 
activity of the ethanol vapor washings is small. 


3. Detecting Unit 


The Geiger-Miiller counter employed consisted 
of a brass cylinder 25 cm in length and 5 cm in 
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diameter with glass ends waxed into place with 
Piceine cement. A five-mil tungsten wire, ce- 
mented at the center of each glass end, was used 
as anode. The standard counter gas mixture was 
constituted of one part of ethanol and three 
parts of argon at a total pressure of 2 cm of 
mercury. The actual counting was operated by 
a “high speed Geiger-Miiller counter Set Model 
G.L. 532” manufactured by Herbach and Rade- 
man. This instrument permits the counting of 0 
to 300 discharges per second. The amplification 
system consists of a five-stage high gain short 
time constant amplifier. The scaling system 
contains circuits employing five pairs of thyra- 
tron tubes yielding scales of 2, 4, 8, 16, and 32. 
An electromechanical watch recording every 
second impulse was used as registering instru- 
ment. 

Under the actual conditions of radioactivity 
measurement, a plateau of about 100 volts was 
observed where the number of counts is practi- 
cally constant with increasing voltages. Meas- 
urements were carried out 50 volts above the 
threshold voltage. 


4. Estimation of Corrections 


(a) Coincidence correction. The Slay cut-off 
coincidence correction formula has been used. 

(b) Background correction. The background 
correction was determined by counting a gas 
mixture containing inactive ethanol vapor and 
argon. 

(c) Sensitivity correction. The sensitivity of the 
counter was checked after each series of counts 
by measurement of the added radioactivity 
caused by the introduction of a radioactive 
material, a bottle of uranium acetate, in a 
definite position near the Geiger-Miiller counter. 


5. Standardization of the Radioactive Water 


Air was saturated with radioactive water at a 
definite temperature (20°C) and, employing 
pressure-volume relationships, a measured quan- 
tity of the saturated air was introduced into the 
Geiger-Miiller counter. Ethanol vapor and argon 
were then successively used, as before, as flushing 
agent until the total pressure in the counter was 
2 cm of Hg. The radioactivity of the standardized 
water was 8.84 10° counts per minute per mole. 








Qo fF = RD oo A 


§ 


wi 


pr 
pr 





th 


ed 
as 
ee 


by 
lel 
le- 


on 


off 
ed. 
nd 
ras 
nd 


ta 
ing 
an- 
the 
ron 
ing 
vas 
zed 


le. 








EXPERIMENTAL RESULTS 


1. As an example illustrating the method, the 
experimental data leading to the calculation of 
the solubility of water in benzene saturated at 
20°C are reported below. The sample of satu- 
rated benzene taken for analysis had a volume 
of 1.10 cc. All the following counts have been 
adjusted to a standard sensitivity of 1000. 

(a) Background correction 383 counts per min- 
ute. 

(b) Ethanol vapor washing of the sample tube 
section. 743 — 383 = 360 c/m. 

(c) Radioactivity of first portions of hydrocarbon 
volatilized from CaO U-tube (0.5 cm Hg in G.M. 
counter each). 1st portion : 1915 —383 = 1532 c/m. 
2nd portion : 843 — 383 =460 c/m. 

(d) Radioactivity of washings of CaO oxide with 
ethanol vapor. ist washing: 8360—383=7977 
c/m; 2nd washing: 4040—383=3657 c/m; 3rd 
washing: 2770—383=2387 c/m; 4th washing: 
1172—383=789 c/m; 5th washing: 900—383 
=517 c/m. Estimated radioactivity remaining 
520 c/m. Total radioactivity 18200 c/m. 

Thus the 1.10 cm* of benzene saturated at 
20°C contain an amount of water corresponding 
to 18200 c/m or 16550 c/m/cc of benzene. 

Measurement of the water radioactivity. The 
amount of water introduced in counters: 2.62 
X10-* moles. Number of discharges: 2697 — 383 
= 2314 c/m. Number of discharges per mole of 
D.O water: 8.84 10-8 c/m. 

Calculation of solubility of ordinary water at 
20°C. This calculation actually requires the 
assumption of the validity of Henry’s Law, a 
knowledge of the Henry’s Law constants and 
also the vapor pressures of the isotopic species 
TDO and H.O. We have assumed the validity 
of Henry’s Law because of the extreme dilutions 
that are here involved. In the absence of data 
we further assume that the Henry’s Law constant 
k= p/n, is the same for TDO and HO. We can 
proceed in the absence of data on the vapor 
pressure of TDO in the following manner. 

The molar concentration of TDO present in 
the sample (volume = 33.4 cc; p=7 cm; t=22°C) 
used in the preceding section for measurement 
of the radioactivity is 


33.4X7X273Xpi 


22400 X 76 X 298 X 760 
= 1.685 X 10-7; moles TDO/mole, 





(1) 
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where /; is the pressure in mm of the TDO in the 
gas sample measured. This is also the pressure 
of TDO in the gas phase in equilibrium with the 
solution when the benzene is saturated with 
water. 

Let A be the number of counts/min./mole of 
TDO which is used. The number of moles in 
question is defined by Eq. (1). The observed 
activity of the sample was 2314 c/m. Hence, 


2314 2314 
——=1.685X10-"p; or pi= 
A 1.685.A 





10’, 


1.375 





10° mm. 


p= 


In 1 cc of benzene, saturated at 20°C, we found 
16550 c/m. The molar concentration of TDO in 
the benzene solution is therefore 


16550 78.05 1.445 
A 0.895 A 





10° mole TDO/mole benzene. 


The Henry’s Law constant, k, for TDO becomes 
therefore 
pi 1.375A-!X101° 


—= = 0.952 X 104. 
m, 1.445A-'X108 





If the assumption is now made that the 
Henry’s Law constant is the same for H.O as 
for TDO the solubility of H:O at 20°C is, 
therefore, 


pe 


te —————— 
0.952 X 104 


where 2 is the vapor pressure of H.O at 20°C, 
i.e., 2 =17.535 mm.!617 
Inserting this value for p2 we receive 


17.535 
ny, = ——_—_——_ 
0.952 x 104 
= 18.4X10-4 mole H2O/mole CgHg. 
The corresponding solubilities in other units are: 


(a) solubility per 100 cc benzene = 2.105 « 10-% 
mole or 3.80 10-*g H:0; 
(b) solubility per 100g benzene = 4.25 X 10-g 
H.O. 
16F, T, Miles and A. W. C. Menzies, J. Am. Chem. Soc. 
58, 1065 (1936). 


7G. N. Lewis and R. T. McDonald, J. Am. Chem. Soc. 
55, 3057 (1933). 
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TABLE I. 
Temper- References 
ature°C This work 3 1 2 4 5 
10 0.030 0.043 0.040 0.035 0.032 0.045 


0.0425 
20 (O04: 0.057 0.056 0.055 0.044 0.057 














26 0.054 0.070 0.070 0.069 0.057 0.067 
TABLE II. 
b/Po 0.60 0.77 0.92 0.95 1.00 
t(monolayers) 1 7.2 46 70 184 








2. Solubilities at 10, 20, and 26°C 


The available literature data on the solubility 
of benzene at 10, 20, and 26°C, extrapolated or 
interpolated when necessary, are reported in 
Table I along with our actual results. 

Examination of Table I reveals that our 
results are in good agreement with those of 
Clifford‘ who measured the water of solution by 
volatilization and absorption by calcium chloride, 
but that those obtained using the other methods, 
synthetic or chemical, are somewhat higher. It 
had been suggested that Clifford’s results were 
possibly too low because the volatilization of 
the water may have been incomplete. Such a 
possibility has been eliminated in our method 
since we have measured the residual radioactivity 
of the sample tube section after distillation of 
the liquid sample. It should also be noted that 
the chemical and synthetic methods tend to 
yield high solubility results. First, any water, 
from atmosphere or other sources, introduced 
during the analytical procedure will increase the 
measured solubility since physical adsorption of 
water on glass surfaces is far from negligible at 
ordinary temperatures. McHaffie and Lenher!® 
among others, have shown that the adsorbed 
water film on glass surfaces is multimolecular. 
Table II summarizes their results; p is the vapor 
pressure of the water in equilibrium with the 
surface film of thickness ¢ at the adsorption 
temperature, and p is the vapor pressure of the 
water at the same temperature. The thickness ¢ 
is expressed in number of monolayers calculated 
assuming that the surface of adsorption is the 
geometrical surface of glass. 


18 T, R. McHaffie and S. Lenher, J. Chem. Soc. 127, 
1559 (1925). 
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These measurements were carried out in the 
temperature interval 25-37°C and it was ob- 
served that the thickness of the film is only a 
function of the water partial pressure. Recent 
investigations!®*° have shown that the real 
surface area of glass varies with chemical, 
thermal, and other treatments of the glass. 
Brunauer suggests that the actual surface area 
of a glass may be from three to ten times greater 
than the apparent or geometrical area. There- 
fore, according to Brunauer, the water film 
thickness reported by McHaffie and Lenher 
would be 18 to 62 monolayers thick at p:po=1. 
It may be assumed that desorption and adsorp- 
tion of water on glass are reversible except for 
the first unimolecular layer which is chemi- 
sorbed and can only be desorbed at high temper- 
ature. The washings of the sample tube with 
ethanol vapor indicate that there was always a 
minimum of residual radioactivity of the order 
of 360 to 400 counts per minute even after 
transfer as complete as possible of the hydro- 
carbon sample into the calcium oxide tube. 
This corresponds to about 4.3X10-7 moles of 
water or 2.59 X 10!" molecules of water. Assuming 
11A? as the area occupied by one water molecule, 
the unimolecular layer would be 285 cm?. The 
estimated geometrical surface area of the sample 
tube section of the apparatus is about 60 cm? 
which gives a roughness factor of 4.75, in good 
agreement with Brunauer’s suggestion. 

The amount of water adsorbed on the walls 
of the glass vessel at saturation pressure will be 
especially important when the solubility is small 
and the ratio of volume of the sample to surface 
area of containing vessel is small. 

When the synthetic method is used, colloidally 
dispersed water may be included with the dis- 
solved water. In the method based on the 
hydrogen liberated by addition of calcium hy- 
dride, Rosenbaum and Walton appear to have 
neglected to consider the amount of gas dissolved 
in the hydrocarbon sample. 

The purity of the benzene employed is another 
factor which may influence the solubility ob- 


19P, Emmett, Advances in Colloid Science (Interscience 
Publishers, Inc., New York, 1942) edited by E. O. Kremer, 
Volume 1, Chapter 1. 

20 Brunauer, The Adsorption of Gases and Vapors. 
Physical Adsorption (Princeton University Press, Prince- 
ton, New Jersey, 1943), p. 323. 
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tained. The benzene which was used in our 
determination was prepared by the Pure Hydro- 
carbon Research Program of the Ohio University 
and had the following physical characteristics :*! 


Density at 20°C, d2°=0.87893. 
Refraction index, np*”® = 1.5012. 
Dielectric constant, E?°=2.2933. 


Preliminary results obtained by use of com- 
mercial benzene during the development of the 
technique indicated solubility results somewhat 
higher than those reported in Table I. The 
increase of solubility due to the presence of 
impurities has been noted by other investi- 
gators®” and willjbe more fully discussed in a 
subsequent paper. 

The reproducibility of the results is indicated 
by the two determinations at 20°C. The loga- 
rithm of the solubilities when plotted against the 
reciprocal of the absolute temperature shows a 
quite satisfactory linear relationship. The heat 
of solution thus obtained is 6.1 Kcal. per mole 
of water. 


Limitations of the Method 


The results are subject to certain sources of 
error. The probable error in determining a 
sample size of the order of about 1 cc may reach 
several percent. The exchange procedure used 
in the determination of the radioactivity is 
never theoretically complete and estimation of 
the residual radioactivity is necessary if an 
exceedingly time-consuming operation is to be 
avoided. The probable error arising in the 


*t Private communication from Dr. C. P. Smyth. 





counting procedure itself may also be of the 
order of a few percent being the function of the 
total radioactivity ; this error can be minimized 
by longer periods of counting or the use of water 
with higher tritium content, or both. 


CONCLUSION 


Fractions of one milligram of water containing 
tritium oxide at a molar concentration of the 
order of 10- can be measured with accuracy. 
This is of particular interest in the determination 
of solubility of water in hydrocarbons available 
only in small quantities in the pure state. This 
method is generally applicable to a wide variety 
of materials but is obviously limited to solvents 
in which no exchange takes place between the 
hydrogen atoms of the water and the solvent 
under the conditions of saturation and analysis. 
It may be applied to determinations over an 
extremely varied solubility range by employing 
water containing adequate tritium contents. It 
has the additional advantage of being easily 
adaptable to high pressure measurements when 
required, as in the case of the lower hydrocarbons 
at room temperature. 
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The author’s data on the properties of chromium electro- 
deposited from a standard chromic acid electrolyte are 
examined in the light of a theory first proposed by Bozorth. 
It is shown that the data fit well the assumption that the 
preferred orientation found in these chromium electro- 
deposits is the result of a slipping process, analogous to 
cold-working, which occurs during deposition when internal 
contractile stress reaches a certain critical value. It is 
further shown that the type of orientation found in 
chromium deposits is in accord with the theory and it is 
stated that (111) orientation is to be expected generally 


for body-centered cubic metals electro-deposited in a state 
of contractile stress, while (110) or (100) is to be anticipated 
for face-centered cubic metals deposited with the same 
type of stress. The orientations of a number of electro- 
deposited metals are shown to be in agreement with this 
expectation. Certain anomalies in the stress and orientation 
results for silver, nickel, and aluminum deposits have also 
been studied experimentally, and the results explained on 
the basis of the above theory. In conclusion, the possible 
mechanism of internal stress formation in electro-deposits 
is examined. 





I. INTRODUCTION 


T was first shown in 1924 by Glocker and 
Kaupp,’ as the result of x-ray diffraction 
analysis, that metals formed by electro-deposi- 
tion could under certain conditions possess an 
orientated structure. The orientations observed 
were of a type known as fiber textures, the metal 
crystals of the deposits being aligned with a 
particular crystallographic axis perpendicular to 
the surface of the cathode but possessing freedom 
of rotation through 360° about this axis. Subse- 
quent investigations by other workers have 
extended the range of electro-deposited metals 
examined by x-ray diffraction and have shown 
that the majority of the commonly deposited 
metals may possess a fiber texture of the type 
first noted by Glocker and Kaupp. These investi- 
gations have also proved that the extent of the 
preferred orientation and to some extent the 
nature of the orientation are both dependent 
upon such variables as composition of electrolyte, 
temperature and current density at which the 
electro-deposition took place, and nature of the 
cathode material. ‘The effect of the cathode 
metal, however, and its influence on the orienta- 
tion of deposits can frequently be ignored. 
While under certain favorable conditions a de- 
posit is able to reproduce the structure and 
lattice spacing of the basis metal on which it is 
produced for some distance, this effect is gener- 
ally very limited, particularly if the deposit 


1R. Glocker and E. Kaupp, Zeits. f. Physik 24, 121 
(1924). 
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itself has a strong orientating tendency. Finch 
and Sun? have shown that where the deposition 
conditions themselves give rise to orientations 
the effect of the substrate does not extend to a 
distance of more than 1000A. Orientation results 
considered in this discussion are limited to those 
found at thicknesses much greater than 1000A, 
results which are dependent solely on the depo- 
sition conditions and do not result from any 
substrate influence. The orientation results for 
chromium, given below, were found to be entirely 
independent of the nature of the substrate, 
whether the substrate was itself orientated or 
unorientated or whether it was a face-centered 
or body-centered cubic crystal structure being 
quite immaterial. 

Fiber textures have been generally ascribed to 
a crystal growth process occurring during elec- 
trolysis, and the general conception of the 
mechanism was admirably summed up by Blum, 
Beckman, and Meyer* when they stated “One 
might expect that preferred orientation would 
tend to occur generally in electro-deposits be- 
cause of the highly anisotropic nature of the 
process—the direction of current flow and 
diffusion restrictions.” 

The existence in electro-deposited metals of 
residual stresses, often of very considerable 


magnitude, was first reported in 1877 by Mills’ 


2G. E, Finch and C. H. Sun, Trans. Faraday Soc. 32, 
852 (1936). 

’W. Blum, A. O. Beckman, and W. R. Meyer, Trans. 
Electrochem. Soc. 80, 409 (1941). 

4E. J. Mills, Proc. Roy. Soc. 26, 504 (1877). 
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who showed, by preparing deposits on the bulb 
of a thermometer and noting the movement of 
the mercury thread caused by the distortion of 
the bulb, that deposits of iron, copper, nickel, 
and silver were formed in a state of contractile 
stress, i.e., were apparently endeavoring to shrink 
and thereby causing the mercury thread to rise, 
while the exact reverse was the case with 
deposits of cadmium and zinc, an expansive 
stress being shown and the mercury thread 
tending to sink. Mills showed that the forces 
involved were considerable, the contractile effect 
of copper on the bulb being as great as a uni- 
formly applied hydrostatic pressure of over one 
hundred atmospheres. Stoney investigated quan- 
titatively the contractile stresses in certain nickel 
deposits and showed them to be as much as 
19.2 tons per square inch parallel to the plane 
of the cathode. More recently, Hume-Rothery 
and Wryllie*’ reported that the contractile 
stresses in chromium deposits may reach values 
of nearly 30 tons per square inch. Much other 
valuable work on the effect of electro-deposition 
conditions on stress has been carried out." 
Relatively little attention has been paid to 
the possible effects of these internal stresses on 
the properties of electro-deposits. It was found 
by Hume-Rothery and Wyllie, however, that 
in the case of chromium deposits there appeared 
to be a clear connection between the internal 
stress and the extent of preferred orientation, 
and between stress and such physical properties 
as hardness and ductility. This paper will seek 
to explain the preferred orientation and prop- 
erties of chromium electro-deposits in terms of 
a plastic deformation process caused by internal 
stress and will endeavor to show that the theory 
may be extended to other electro-deposited 
metals. An attempt will also be made to explain 
the mechanism of the formation of both con- 
tractile and expansive stresses in electro-deposits. 


5G. G. Stoney, Proc. Roy. Soc. A82, 172 (1909). 

®W. Hume-Rothery and M. R. J. Wyllie, Proc. Roy. 
Soc. A181, 331 (1943). 

7™W. Hume-Rothery and M. R. J. Wyllie, Proc. Roy. 
Soc. A182, 415 (1943). 

8V. Kohlschutter and E. A. Vuilleumier, Zeits. f. 
Elektrochemie 24, 300 (1918). 

® V. Kohlschutter and F. Jakober, Zeits. f. Elektrochemie 
33, 290 (1927). 

10E. A. Vuilleumier, Trans. Am. Electrochem. Soc. 42, 
99 (1922). 
(1935) Marie and N. Thon, Comptes rendus 193, 31 
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II. THE DATA FOR ELECTRO-DEPOSITED 
CHROMIUM 


Hume-Rothery and Wyllie*® showed that the 
hardness of chromium electro-deposits produced 
from an electrolyte comprising 250 grams per 
liter chromic acid and with a CrO; to sulphate 
ion ratio of 100:1 was a property solely de- 
pendent on the structures of the deposits as 
revealed by x-rays and was independent of the 
exact conditions under which the structure of a 
given type was produced. Further work has 
shown that the principle holds good for the 
two other physical properties of ductility and 
internal stress. Thus with increasing temperature 
of deposition a sequence of structures is found at 
constant current density, each structure being 
characterized by its own value for hardness, 
ductility, and internal stress. With increasing 
current density the temperature at which a given 
structure occurs increases, and it appears that 
the relationship between the temperature and 
current density producing a structure of a 
particular type is an exponential one. This 
relationship, in the case of the most highly 
preferred type of chromium deposit, has already 
been pointed out.® 
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TEMPERATVRE °C 


Fic. 1. Properties of chromium electro-deposits; electro- 
lyte 250 g per liter CrOs: 2.5 g per liter [SO,’’]; current 
density: 1000 amp. per square foot. 


















The structure and properties of chromium 
deposits produced over a wide range of temper- 
ature at a constant current density may thus be 
taken as representative of all deposits obtainable 
from the particular chromic acid electrolyte used 
in the experiments. Figure 1, which has been 
adapted from the papers of Hume-Rothery and 
Wyllie, and Wyllie,” shows on the same diagram 
curves for the hardness, ductility, and internal 
stress of electro-deposits produced on a basis 
metal of carbon steel with increasing temperature 
of deposition at the constant current density of 
1000 amp. per square foot. The diagram has 
also been divided into four zones to represent 
the changes in the orientation of deposits 10 
microns thick which, as revealed by x-rays, 
occur with increasing temperature. In Zone | 
the deposits have a random orientation, but at 
the border of Zone I and Zone II a trace of 
orientation of the (111) type is first noticed. 
With increasing temperature in Zone II two 
effects are apparent: (a) the number of particles 
with purely random orientation decreases and 
(b) the perfection of the alignment of crystals 
with the (111) orientation increases. In the 
shaded zone lying between 623°C and 673°C 
there are no crystals with random orientation. 
In this region the appearance of the crystals is 
very bright and they all have a very strong 
(111) orientation, the deviation from a perfect 
(111) orientation not exceeding -7}°. 

In Zone III effects precisely opposite from 
those operating in Zone II are evident. Here 
increase of temperature brings about a progres- 
sive decrease in the exactness of alignment of 
the (111) orientation and an increase in the 
number of particles with a purely random 
orientation. At the border Zone III/Zone IV all 
orientation effects finally disappear, and the 
high temperature Zone IV deposits have a 
random structure. 

The changes of structure revealed by x-rays 
are followed, particularly noticeably, by the 
curve giving the values for the residual con- 
tractile internal stress in deposits of the same 
10-microns thickness. In Zone I stress goes up 
rapidly to a peak of nearly 30 tons per square 
inch at about the Zone I/Zone II border. In 


2M, R. J. Wyllie, Trans Electrochem. Soc. Preprint 
92-5 (1947). 
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Zone II it decreases to a minimum and almost 
zero value for the deposits of maximum (111) 
orientation. In Zone III the stress increases 
again to a smaller maximum of about one-third 
the previous value at the Zone III/Zone IV 
border, and then decreases again rapidly in 
Zone lV. 

Ductility values are also at a minimum in the 
case of the most orientated deposits, but ductility 
increases rapidly in the case of deposits prepared 
on both the high and low temperature sides of 
the region of maximum orientation. Hardness, 
too, has a maximum value coincident with the 
region of maximum orientation and falls off in 
value at both higher and lower temperatures— 
the more rapidly, however, to the high tempera- 
ture side. 

These data appear to indicate that there is 
some definite correlation between the physical 
properties of chromium electro-deposits and their 
preferred orientation and particularly between 
the internal stress and orientation. The correla- 
tion appeared to the writer to exceed any effects 
which might be attributable merely to an 
inherent tendency for the deposit to grow 
preferentially under certain deposition conditions 
in a (111) direction, and an attempt is therefore 
made to explain the correlation on the basis of 
internal stress. 


III. THE RELATIONSHIP BETWEEN STRESS AND 
TYPE OF ORIENTATION IN ELECTRO- 
DEPOSITS 


A theory attributing the formation of preferred 
orientation in electro-deposits to internal stress 
was first put forward by Bozorth.” A similar 
idea had earlier been incorporated in a paper by 
Greaves,“ but Greaves had claimed that pre- 
ferred orientations would be expected only in 
deposits produced at low current densities while 
the fine-grained hard deposits resulting from 
high current density deposition would be entirely 
random in structure. Greaves attributed internal 
stress to the effect of the electric field operating 
during deposition and claimed that only at low 
current densities, when the atoms had had time 
to arrange themselves, would orientation effects 
ensue. The data for chromium are entirely out 


13 R, M. Bozorth, Phys. Rev. 26, 390 (1925). 
4 R. H. Greaves, Metallurgist 1, 141 (1925). 
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of accord with the views of Greaves, since in the 
case of chromium it is the hardest deposits of 
minimum grain size which exhibit the maximum 
preferred orientation. Bozorth considered that 
since deformation produced in metals by cold- 
working gave rise to special orientations, it was 
natural to suppose that the tension known to be 
present in many electro-deposited metals pro- 
duced orientations in them in a manner similar 
to that in which orientations were produced in 
hard-drawn wires or rolled sheets. Bozorth 
considered that the stress was produced by 
co-deposited hydrogen which diffused out of the 
metal as the deposit thickened, thus leaving the 
metal in a state of tension. He supposed that as 
the process of deposition went on the stress 
difference between the innermost and outermost 
layers of the deposit increased and that ulti- 
mately the compressive stress in the layers next 
to the cathode became so great that the crystals 
were crushed, thus forming new crystals of very 
small dimensions. In the outer layers, however, 
where the crystals were formed under greatest 
tension, the orientations would be most marked. 
Bozorth’s idea thus also accounts for the experi- 
mental observation that preferred orientation 
tends to increase with increasing deposit thick- 
ness. Experimentally, Bozorth examined some 
deposits of iron, nickel, copper, zinc, and cad- 
mium and found marked orientation effects only 
in iron and nickel. These two metals were found 
to be much harder than the deposits of random 
orientation, and this fact Bozorth interpreted as 
lending support to his view that the orientations 
observed in the iron and nickel were produced by 
some sort of cold-working mechanism. 

The theory of Bozorth, although clear in its 
general implications, i.e., that stress in electro- 
deposits is produced by co-deposited hydrogen 
and that this stress is responsible for the pro- 
duction of preferred orientation in the deposits 
by a mechanism analogous to the cold-working 
of metals, was put forward in a somewhat 
qualitative manner and was based on rather 
scanty experimental data. The more compre- 
hensive results available for chromium deposition 
enable the writer to carry Bozorth’s theory a step 
further. 

The deposition of chromium is an extremely 
inefficient process, and in no case were efficiencies 


greater than 35 percent noted, the most orien- 
tated deposits being characterized by an effici- 
ency of about 9 percent.* The bulk of the current 
is thus consumed in the liberation of hydrogen, 
and for every atom of chromium laid down at 
the cathode numerous atoms of hydrogen are 
simultaneously liberated. It may be assumed 
that some of these hydrogen atoms are initially 
entangled in the growing crystal lattice of 
chromium, causing the lattice to be abnormally 
expanded. As deposition proceeds the major 
portion of the hydrogen in the lattice diffuses 
out, thus leaving the lattice in an expanded 
and unstable condition. The chromium atoms, 
therefore, attempt to contract together to their 
normal interatomic positions. This contraction, 
however, is opposed by the adhesion of the 
deposited atoms to the layers below, and a state 
of strain is thereby set up. It is not suggested 
that all the hydrogen initially contained in the 
lattice is lost since a residual expansion of the 
lattice of electro-deposited chromium has been 
reported'® and Hume-Rothery and Wyllie, al- 
though unable to determine the lattice parameter 
of chromium in the ‘‘as-deposited” condition 
because of the extreme diffuseness of the lines on 
a Debye-Scherrer film, nevertheless found a 
parameter of 2.8861A at 23°C, as against the 
normal value of 2.8786A, in the case of a highly 
orientated deposit annealed in vacuum for 4 
hours at 750°C. The suggestion is put forward, 
however, that immediately after deposition, 
when the bulk of the hydrogen escapes, the 
lattice is considerably strained and that this 
strain is relieved, when it reaches a certain 
critical value, by an actual slipping or yielding 
of the metal crystals which are thus left in a 
condition somewhat analogous to that of a cold- 
worked metal. On this assumption it may be 
understood why the temperature and current 
densities leading to deposits with the most highly 
preferred orientations are those which lead also 
to the lowest residual stress, and why such 
deposits are hard and brittle, since it is known 
that cold-work increases the hardness and 
decreases the ductility of a metal. It has been 
shown that the stresses involved in chromium 
deposition are very high, and at the moment of 


18 W. A. Wood, Phil. Mag. 23, 984 (1937). 





















deposition these may well be sufficient to produce 
plastic deformation. While the stresses involved 
are thus of sufficient magnitude to justify the 
concept of plastic distortion occurring, the 
amount of actual strain involved in the process 
is more difficult to assess. The distortion neces- 
sary to account for stresses of the magnitude 
actually observed experimentally in deposits 
would result from an enforced straining of the 
deposit amounting to perhaps 0.7 percent in di- 
rections parallel to the plane of the cathode. This 
order of magnitude of straining would not 
normally account for very appreciable internal 
slip in the crystals of the deposits, and certainly 
would not be expected to yield the slip into 
preferred orientations which is postulated below. 
Reduction in volume of the order of 60 percent 
is necessary if a block of massive metal is to be 
so distorted by slip as to show evidence of a 
preferred orientation. While it is realized that 
this fact constitutes a potent argument against 
the slip theory of Bozorth and the writer, it may 
be argued that the exact mechanism of slip is 
still obscure and that in the case of a massive 
block of polycrystalline metal considerable re- 
duction in volume may be essential before the 
influence of the strain on small crystals becomes 
sufficiently appreciable and directed so as to 
produce orientation effects. It is possible that 
the separating action of co-deposited hydrogen 
produces during deposition a very much greater 
strain throughout very small crystals than can 
ever be detected in a deposit completed and 
removed from the electrolyte. If this be so, it is 
suggested that during the contraction of these 
expanded crystals to approximately normal 
interatomic distances after loss of hydrogen slip 
occurs over shorter ranges and very much more 
freely than in a massive metal, this slip resulting 
both in a harder and more orientated state of 
the material than can be achieved by any ordi- 
nary cold-working process. On this theory slip- 
bands would not normally be visible in an 
electro-deposited metal, and the so-called slip 
bands observed by Hughes'* and Tamman!’ in 
iron and copper deposits, as well as the rhythmic 


mf E. Hughes, J. Iron and Steel Inst. 103(1), 356 
1921). 

17G, Tamman, Lehrbuch der Metallkunde (Leopold Voss, 
Leipzig, 1932). 
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banding observed in nickel,!® would be suscep- 
tible to explanations other than as confirmation 
of the internal slip theory suggested. 

On the theory outlined above, the data shown 
in Fig. 1 are readily explicable. Across Zone I 
the residual stress resulting from the lattice 
strain caused during deposition is increasing, and 
the general distortion of the crystal lattice causes 
both an increase in deposit hardness and a 
diminution in ductility, but not until the Zone 
I/Zone I! boundary is reached does the lattice 
strain during deposition attain the critical value 
for sufficient slip to occur to manifest itself in 
preferred orientation. In Zone II it may be 
presumed that the strain set up during deposition 
steadily increases to reach a maximum value in 
the shaded region. In this zone, therefore, 
orientation increases and the residual stress, 
indicative of strain unrelieved by slip, decreases. 
Hardness correspondingly increases and ductility 
further decreases. In Zone III it would appear 
that the strain resulting immediately after de- 
position begins to decrease, possibly because of 
the effects of increased thermal vibration ena- 
bling the co-deposited hydrogen to escape more 
readily. As the strain decreases in magnitude 
less slip, and therefore preferred orientation, 
results and more unrelieved strain is left in the 
lattice to show itself as residual stress. In Zone 
IV the strain never reaches a value sufficient to 
cause any orientation and appears to fall off 
rapidly with increasing temperature. In Zones III 
and IV hardness also decreases and ductility 
increases as would be anticipated on the basis 
of the strain effects postulated. 

It remains to be seen what effect the slipping 
process suggested would have on the crystal 
orientation, and here an apparent correlation 
has been found between the textures of electro- 
deposited metals and of the same metals in the 
drawn state. 

It is well known that when a metal is drawn 
the resulting texture is such that one or more 
crystallographic directions lie in the direction of 
drawing. It has been found that if for a particular 
metal the drawn texture is such that one or 
more crystallographic directions lie in the direc- 
tion of drawing, then the texture of the electro- 


18 See, for example, W. A. Wesley and E. J. Roehl, 
Trans. Am. Electrochem. Soc. 82, 37 (1942). 
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deposit, if the latter involves formation with 
contractile stress, is such that the maximum 
number of the same directions lie in the plane of 
the deposit. If, however, the electro-deposit 
involves expansive stress, the texture is related 
to the compression texture of the metal. 

In the case of body-centered cubic metals it 
is known that the drawing texture is such that a 
[110] direction lines in the drawing direction, 
and the above principles therefore require the 
maximum number of [110] directions to lie in 
the plane of the deposit. This condition is 
satisfied by a (111) orientation since the (111) 
planes contain three [110] directions and this is 
the maximum possible. The preferred orienta- 
tions of body-centered cubic chromium and of 
iron agree with this principle. That the orienta- 
tion in iron is (111) has been shown,"!9 and that 
the stress is contractile has also been demon- 
strated.* 9 

With face-centered cubic metals, drawing pro- 
duces one of two different textures, or both, with 
the [100] or [111] directions parallel to the 
direction of drawing.* The above principle leads 
to the expectation that electro-deposited face- 
centered cubic metals will possess preferred 
orientations with the greatest number of [100] 
and [111] directions in the plane of the cathode 
if the deposits are formed with contractile 
stress. Taking the different planes in order it is 
found: 


No. of [100] direc- 
tions in cathode 


No. of [111] direc- 


Plane parallel to tions in cathode 


cathode surface surface surface 
(100) 2 0 
(110) 1 2 
(111) 0 0 
(210) 1 0 
(211) 0 1 


The planes of the higher indices are all less 
favorable as regards the number of [100] and 
[111] directions they contain. 

It should be expected that the preferred 
orientations of electro-deposited face-centered 
cubic metals will be (100) or (110) if the deposits 


_%G, I. Finch, A. G. Quarrell, and H. Wilman, Trans. 
Faraday Soc. 31, 1051 (1935). 

* With aluminum the texture is almost entirely [111] 
whilst copper gives a mixed texture 40 percent [100] and 
60 percent [111] and silver a mixed texture 75 percent 
[100] and only 25 percent [111]. Nickel resembles copper. 


involve contractile stress. In agreement with 
this principle electro-deposits of cobalt formed 
under conditions likely to give rise to contractile 
stress® have been found with a (110) orientation.” 
Similarly deposits of copper have a (110) or 
(100) orientation’?° when formed under condi- 
tions giving rise to contractile stresses. !° De- 
posits of nickel exhibit very strong contractile 
stress,*5 and the orientations generally reported 
are (110) or (100), although (211) has been 
noted." It has been stated*! that a (100) orienta- 
tion is always obtained if the electrolyte is 
stirred by rotating the cathode. The orientation 
of gold electro-deposits has been shown to be 
(110),!° but no reference to the nature of the 
internal. stress appears to exist. The data for 
silver are more complicated. Glocker and Kaupp 
stated! that both [100] and [111] orientations 
were observed. On the theory above a [111] 
orientation would not be expected from a silver 
deposit possessing contractile stress. However, 
it has been demonstrated® that, depending on 
the current density of deposition, silver deposits 
formed from the same solution may have either 
contractile or expansive stress. Since the com- 
pression texture of silver is such that the [110] 
direction is the direction of compression, a silver 
deposit formed under conditions of expansive 
stress might be expected to possess a (111) 
orientation since the (111) plane, as previously 
observed, contains three [110] directions. 

Experiments, described below, were therefore 
undertaken to investigate the connection be- 
tween the direction of stress and orientation in 
silver deposits. Some experiments involving the 
deposition of nickel from a nickel chloride 
solution under various conditions, as well as 
aluminum from an organic electrolyte, are also 
noted. 


IV. STRESS AND ORIENTATION IN SOME SILVER, 
NICKEL, AND ALUMINUM DEPOSITS 


Silver 


Kohlschutter and Jakober® showed that silver 
deposits produced on a platinum cathode from a 


20H. Hirata and Y. Tanaka, Mem. Coll. Sci. Kyoto 
Imp. Univ. Al5, 9 (1932). 

21'W. G. Burgers and H. E. Quanjel, Rec. trav. chim. 
61, 353 (1942). 
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(a) 


(b) 


Fic. 2. (a) Silver deposit showing (111) orientation ; expansive stress; current density 5 M.A./cm?; temperature 18°C 
electrolyte: 0.2 M AgCN, 0.4 M KCN. (b) Silver deposit showing (110) orientation; contractile stress; current density 
5 M.A./cm?; temperature 18°C; electrolyte: 0.2 M AgCN, 0.2 Af KNOs, 0.4 M KCN. 


solution comprising 0.2 M AgCN and 0.4 M KCN 
possessed strong expansive stress at a current den- 
sity of 5 milliamp. per square centimeter and con- 
siderable contractile stress at a current density of 
15 milliamp. per square centimeter. It was decided 
to reproduce their results and to examine the de- 
posits by the glancing beam x-ray method.” One 
liter of electrolyte was made up 0.2 M AgCN and 
0.4 M KCN from British Drug Houses “‘Analar”’ 
silver cyanide and potassium cyanide. The con- 
tainer used was a liter glass beaker. The anode, 
which closely fitted the whole interior wall of the 
beaker, was silver foil of 99.99 percent purity. 
The cathode was silver foil of the same quality 
as the anode and was 0.05 mm thick and 5 cm 
wide. This foil had been delivered by its manu- 
facturers in the form of a tight roll and was 
naturally smooth and springy. When a length of 
11 cm was cut from the roll, it curled into a 
circular segment of about 3 cm diameter, there 
being a gap of just over 1 cm between the edges 
of the split cylinder thus formed. Advantage 
was taken of this phenomenon to measure con- 
veniently the type of stress in deposits. A copper 
lead was soldered to the inside of the cylinder 
at a point diametrically opposite the split, and 
the entire inside of the silver foil as well as the 
copper wire was stopped off with cellulose 
varnish. The cylinder was then suspended in the 
electrolyte by means of the wire so that it was 
roughly concentric with the silver anode around 


2M. R. J. Wyllie, Rev. Sci, Insts. 18, 425 (1947). 


the wall of the beaker. With this arrangement, 
since deposition can take place only on the 
outside of the split cylinder, it is clear that any 
expansive or contractile stresses in silver deposits 
formed will manifest themselves in a closing or 
opening of the gap in the cylinder. 

At a current density of 5 milliamp. per square 
centimeter and at 18°C a definite closing of the 
gap was noted. The effect was not very large, 
amounting to about 1-cm movement for a 
deposit thickness of 10 microns. At 10 milliamp. 
per square centimeter no stress was perceptible, 
and this was also the case of 15 milliamp. per 
square centimeter. The latter deposit was a 
dirty yellow in color, nodular, and rather non- 
adherent. No evidence of the contractile stress 
claimed by Kohlschutter and Jakober for a 
deposit formed at this current density could be 
detected. X-ray evidence revealed that the 
deposit produced at 5 milliamp. per square 
centimeter had a marked (111) orientation (Fig. 
2a) but those formed at 10 and 15 milliamp. per 
square centimeter were random in orientation. 

Further investigation showed that contractile 
stress in silver deposits on a silver cathode could 
be formed from a solution containing 0.2 M 
AgCN, 0.2 M KNOs, and 0.4 M KCN. This 
solution was made from the previous one by 
addition of ‘‘Analar’’ brand potassium nitrate. 
Using the same arrangement of anode and 
cathode as described previously the following 
results were obtained at a current density of 5 
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milliamp. per square centimeter and an electro- 


lyte temperature of 18°C: 


Distance between the sides of the 


Time gap in the split cylinder Expansion 


(minutes) (millimeters) (millimeters) 
0 16.21 — 
6 18.54 2.32 
12 23.70 5.16 
18 29.38 5.68 
31 38.27 8.89 


Total expansion 22.06 


X-ray investigation of this deposit (Fig. 2b) 
showed that it possessed a marked (110) orienta- 
tion. These data for silver thus tend to sub- 
stantiate the theory outlined above, insomuch 
as deposits with contractile stress show a (110) 
orientation and those with expansive stress a 
(111) orientation. 

Further experiments in silver deposition were 
conducted to ascertain whether or not the lattice 
parameter and crystal form of deposits were the 
same during deposition, i.e., while the current 
was flowing, as they were after deposition was 
complete. A possible cause of internal stress may 
be the deposition of one allotropic form of a 
metal, this form reverting to a more stable 
modification of different volume either during 
deposition or after removal of the specimen from 
the electrolyte. Such changes would normally 
escape detection as deposits are examined some 
considerable time after deposition has ceased. 
The allotropic modification theory for stress will 
be further discussed below. 

The apparatus used is shown diagramatically 
in Fig. 3. The electrolyte was 0.2 M AgCN and 
0.4 M KCN contained in a rectangular glass 
vessel 10 cm X6 cmX8 cm. The anode was of 
99.99 percent silver sheet 10 cm X8 cm in area 
arranged against the side of the vessel as shown 
in the figure. The cathode was a 99.99 percent 
silver disk 0.2 cm thick and 6 cm in diameter 
and was screwed to the shaft of an electric 
motor. The disk was rotated at i r.p.m. The 
bottom of the disk was 2 cm below the surface 
of the electrolyte. With this arrangement, depo- 
sition of silver at 18°C was carried on at a 
constant current density of 5 milliamp. per 
square centimeter on the area of the cathode 
immersed. The portion of disk plated, of course, 
changed all the time as the disk revolved, but 
since the portion of the disk not in the electrolyte 








remained covered with a film of electrolyte as it 
revolved, it is possible that some slight deposition 
always occurred on this portion as well as on the 
immersed portion. It was, in fact, hoped that 
this was the case. Normal back reflection photo- 
graphs, using a revolving circular cassette holder 
with 120° masks and unscreened cobalt radiation, 
were taken of (a) the pure silver cathode, (b) the 
deposit while depositing, and (c) the deposit 
after deposition ceased. The disk was not touched 
during this sequence of photographs nor was the 
electrolyte removed. The beam was arranged to 
impinge 0.5 cm from the edge of the disk and 
1 cm above the electrolyte surface, and the 
direction of disk rotation was such that the 
portion of deposit diffracting was that which 
had just emerged from the electrolyte. Disk to 
film distance was approximately 3 cm. 

Figure 4 shows the result obtained. It is 
evident that under these conditions there is no 
difference between the lines given by the deposit 
during deposition or after deposition has ceased. 
The lines from the deposit are, however, more 
diffuse than those from the basis metal, and 
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Fic. 3. Schematic tence he of Lae for taking back 
reflection x-ray photograp silver electro-deposits 
immediately after deposition. 





















there is possibly a very slight general increase in 
the lattice parameter of the deposited silver. 
The effect is not large. It is suggested that the 
method outlined above deserves wider applica- 
tion to the study of electro-deposition processes 
generally. 


Nickel 


The solution chosen was one which has been 
stated'* to show particularly large values of 
contractile stress. It comprised 205 g per liter of 
NiCl,-6H.O, 8 g per liter of KCl and 17 g per 
liter of NH,Cl. The salts used were ‘“‘Baker’s 
Analysed” reagents, the cobalt content being 
0.10 percent, the copper 0.003 percent, and the 
iron 0.003 percent. Measurement of stress was 
carried out by depositing nickel on one side of a 
flexible high carbon steel strip and measuring the 
direction and magnitude of the curvature of the 
strip after the deposition was complete. The 
steel used was of Swedish manufacture and 
remarkably uniform in its properties. It was 
0.25 in. wide and 0.002 in. thick. Bending was 
measured by determining the distance between 
the free end of a 14 cm long strip (12 cm being 
immersed in the electrolyte) relative to a fixed 
glass pointer when the other end of the strip was 
clamped in a brass holder. Before deposition the 
strip was cleaned by cathodic etching in a 
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Fic. 4. Back reflection x-ray photograph of silver while 
depositing; current density 5 M.A./cm* temp. 18°C; 


electrolyte: 0.2 M AgCN, 0.4 M KCN; segment 1: 
deposition anaewryes current on; segment 2: deposition 
concluded, current off; i 


segment 3: pure silver basis metal. 
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solution of sodium carbonate followed by an 
anodic etch of 30-seconds duration at a current 
density of 200 amp. per square foot in a 30 
percent sulfuric acid solution. The reverse side 
of the strip was insulated with cellulose varnish. 
Deposits were prepared in all cases at a thickness 
of 10 microns, this figure being checked by 
weighing. A pH of 5.0, measured with a com- 
mercial glass electrode, was standardized. 

The electrolyte was contained in a glass vessel 
10 cmX5 cm X15 cm, the temperature being 
controlled to 0.1°C by a constant temperature 
bath. The nickel anodes used were 14 cm X5 cm 
in area and were of high purity, the material 
being donated by the International Nickel Com- 
pany. 

Initial experiments were carried out at a 
current density of 5 amp. per square decimeter 
and were made to ascertain the effect of temper- 
ature on stress. Figure 5 shows the results 
obtained. It is obvious that, while the stress tends 
to fall off steadily with increasing temperature 
of deposition about 25°C, there is a very rapid 
and immediate decrease below 25°C. Figures of 
stress are given in millimeters of deflection of the 
tip of the strip, but it may be pointed out that a 
deflection of 1 mm represents a contractile stress 
of about 6 tons per square inch parallel to the 
plane of the cathode. X-ray investigations 
showed that in all cases the orientation of these 
nickel deposits was random. 

In terms of the theory outlined above the 
sudden change of stress below 25°C would be 
attributable to an internal relief of strain by slip 
into a preferred orientation. No orientation was 
found. However, under the microscope the de- 
posits below 25°C all showed numerous cracks, 
and it was evident that stress was being relieved 
in these deposits by rupturing of the metal film 
and not by any possible internal mechanism. 
The cracking of deposits, with consequent relief 
of stress, has been pointed out” to be a source of 
error in the measurement of stress in nickel 
deposits. It was a complication not noticed in 
the case of chromium. 

The hardness of nickel deposits has been 
attributed largely to co-deposited basic nickel 


%D. J. MacNaughtan and A. W. Hothersall, Trans. 
Faraday Soc. 31, 1168 (1935). 
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salts,%?5 and considerable evidence to justify 
this belief has been presented. At an operating 
pH of 5.0 the H of the cathode film is un- 
doubtedly much higher, as has been experi- 
mentally observed,”* and co-deposition of nickel 
hydroxide, formed as a sol adjacent to the 
cathode surface, is almost inevitable. It was 
considered that the effect of co-deposited basic 
material highly dispersed in the deposit might 
be to weaken, to some extent, intergranular 
cohesion and therefore to promote cracking as 
well as to inhibit plastic deformation. It is known 
that the effect of a finely divided second phase 
in a matrix is to make plastic deformation, and 
therefore the occurrence of orientation, more 
difficult. 

The only logical method of cutting down the 
deposition of basic material appeared to be by 
decreasing the thickness of the cathode film, and 
this was achieved by affixing anode and cathode, 
facing each other, at the end of an arm 5 cm in 
length and rotating this arm about its center in 
the electrolyte at 25 r.p.m. Under these condi- 
tions at 5 amp. per square decimeter, crack-free 
deposits of low stress values could be fairly 
reproducibly obtained below 25°C, and these 
uncracked, fine-grained deposits showed evidence 
of a (110) orientation. Deposits obtained by the 
rotating anode and cathode method at tempera- 
tures above 25°C did not differ appreciably either 
in appearance, internal stress, or orientation 
from those prepared under the same conditions 
of current density and electrolyte in the absence 
of rotation. 

A point of interest noted when determining 
stress was the fact that in many cases when the 
deposit was being soaked in an alcohol/acetone 
mixture to remove stop-off varnish prior to the 
bending of the strip being measured, a visible 
evolution of gas from the deposit occurred, while 
at the same time further bending of the strip, 
i.e., contraction of the deposit, took place. It 
was often necessary to wait several hours before 
the strip bending reached a final equilibrium 
value. 

The results for nickel deposited from a nickel 


*D. J. MacNaughtan and R. A. F. Hammond, Trans. 
Faraday Soc. 27, 633 (1931). 

25D. J. MacNaughtan, G. E. Gardam, and R. A. F. 
Hammond, Trans. Faraday Soc, 29, 729 (1933). 
*6 A, Brenner, Proc. Am. Electroplaters Soc. 28 (1941). 
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Fic. 5. Stress in nickel deposits from nickel chloride 
electrolyte; current density 5 amp. per square decimeter 
pH 5.0. 


chloride electrolyte under these conditions indi- 
cate that very high values of contractile stress 
can occur and that a peak stress is reached at 
about 25°C whether the electrolyte is stationary 
or flowing past the cathode. Stress relief below 
25°C is generally caused by cracking of the 
deposit, but if rotation is used the stress is still 
low and an orientation is observed. This result 
tends to confirm the statement™ that nickel 
deposits from a nickel acetate electrolyte cannot 
be produced with a definite (100) orientation on 
a circular cathode of 6-mm diameter unless the 
cathode is rotated at 100 r.p.m. or more. It was 
found impossible to obtain satisfactory x-ray 
photographs of nickel deposits while deposition 
was being carried on. 


Aluminum 


Aluminum was deposited on steel strips 
cleaned in precisely the same way as that used 
previously for nickel deposition. The non-aque- 
ous electrolyte, essentially aluminum bromide 
and ethylene bromide, was made up in the 
manner suggested by Blue and Mathers.?’ Pure 
aluminum was used as anode and the electrolyte, 
anode, and steel cathode were enclosed in an 
airtight cell 4 cm in diameter and 16 cm high. 
Using a current density of 1 amp. per square 
decimeter at room temperature, good adherent 


27R. D. Blue and F. C. Mathers, Trans. Am. Electro- 
chem. Soc. 65, 339 (1934). 
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and bright aluminum deposits were obtained 
after the electrolyte had been worked for some 
time. No perceptible stress, either expansive or 
contractile, was recorded, and x-ray diffraction 
showed that the deposits were of a grain size of 
10-? to 10-* cm and were random in orientation. 
This appears to be the first occasion on which a 
determination of the stress and orientation of 
deposits from non-aqueous solvents has been 
attempted. 


V. THE ORIGIN OF STRESS IN ELECTRO- 
DEPOSITS 


A theory which attributes the orientation of 
many electro-deposits to stress has already been 
given, and a mechanism for the possible forma- 
tion of the contractile stress in chromium electro- 
deposits has been outlined. The deposition of 
chromium from an aqueous chromic acid electro- 
lyte is, however, unique among the commonly 
deposited metals insomuch as the deposition 
efficiency is extraordinarily low and in all cases 
liberation of chromium is accompanied by the 
evolution of very considerable quantities of 
hydrogen. It remains to be considered whether 
the co-deposited hydrogen theory is capable of 
explaining contractile stresses generally, and 
whether any analogous theory can account for 
the formation of expansive stresses in certain 
electro-deposited metals. 

Although a number of papers have been 
published on the experimental determination of 
stress in deposits, stress data are still scanty. 
Certain generalizations on the basis of existing 
results may, however, be made. 


(1) Only contractile stresses are exhibited by metals of 
the transitional groups unless very considerable quan- 
tities of organic addition agents are added. E.g., the addi- 
tion of sufficient sodium napthalene trisulfonate to a nickel 
bath may eventually produce a slight expansive stress. 

(2) Metals such as copper and silver, which are closely 
related to transitional metals, may exhibit either con- 
tractile or expansive stresses depending on the deposition 
conditions. 

(3) Non-transitional metals such as lead, zinc, and 
cadmium, if stressed at all, show only expansive stresses. 


It would also appear that contractile stresses 
tend to be of considerably greater magnitude 
than expansive stresses. 

The fact that hydrogen is associated with the 
formation of contractile stresses appears to 
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derive support from a number of considerations. 
It is well known that transitional metals, and 
under certain conditions semitransitional metals 
such as copper and silver, tend most readily to 
form stable or metastable alloys with hydrogen. 
Expansion of the lattices of such metals during 
deposition and their subsequent contraction is 
therefore quite feasible. Residual expansions in 
electro-deposited chromium have been found in 
this work, and the same is certainly true for 
palladium.”*?° Depolarizers and alternating cur- 
rent superimposed on the direct plating current 
both serve markedly to reduce contractile stress, ° 
and their action is almost certainly the oxidation 
of hydrogen shortly after it is liberated cathodi- 
cally and before it can exert any considerable 
effect on the growing metal lattice. The work of 
Marie and Thon" shows that copper deposited 
from a neutral solution of copper sulphate 
exhibits expansive stress on bright platinum, 
silver, and gold substrates but a contraction on 
a dull platinum base. When the addition of 1 
percent sulphuric acid is made contractions occur 
on both dull and bright platinum, and in the 
presence of 2 percent H2SO, contractions are 
observed on all the substrates. These data point 
both to an over-voltage effect and the relative 
ease of liberation of hydrogen and metal. Con- 
tractile effects increase with time after removal 
of the deposit from the electrolyte and are 
accompanied by the liberation of a gas which is 
almost certainly hydrogen. The annealing of 
highly stressed deposits at low temperatures 
frequently increases contractile stress, and gas 
is evolved. There is, nevertheless, no direct con- 
nection between magnitude of stress and amount 
of co-deposited hydrogen. Thus in Fig. 1 the 
efficiency of deposition at 1000 amp. per square 
foot is about 30 percent at 30°C falling to only 
5 percent at 85°C, while the values for stress 
vary as shown. A similar lack of obvious rela- 
tionship exists in the case of nickel.*° It has 
recently been shown*! that the stress in nickel 
deposits increases with increase in the pH, and 
on these data it was claimed that co-deposited 


28. McKeehan, Phys. Rev. (2) 21, 334 (1923). 

29G. Chaudron, J. Bérnard, and A. Michel, Comptes 
rendus 218, 913 (1944). 

30D. J. MacNaughtan and A. W. Hothersall, Trans. 
Faraday Soc. 24, 387 (1928). 
31 B. Martin, Proc. Am. Electroplaters Soc. 207 (1944). 
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basic material is the cause of stress. Nevertheless, 
it would appear that some hydrogen is always 
evolved if contractile stress is to result, and 
nickel deposits prepared at a pH of as much as 
5.7 have been found by direct analysis to contain 
hydrogen.” A tentative conclusion appears to 
be that hydrogen must be evolved before con- 
tractile stress manifests itself, but that there is 
an optimum amount of hydrogen necessary to 
produce maximum stress, the amount varying 
with temperature, the rate of deposition, and 
other factors not readily capable of quantitative 
evaluation. 

The formation of expansive stresses in silver, 
copper, zinc, cadmium, and lead deposits under 
appropriate conditions is more difficult to ex- 
plain. It may be that both contractile and 
expansive stresses result from the deposition of 
one allotropic modification of a metal, this 
modification reverting to a more stable form, 
with a consequent volume change, while deposi- 
tion proceeds. No direct experimental evidence 
has yet been found to justify this theory, but 
it is perhaps noteworthy that metals such as 
nickel, chromium, iron, and cobalt, which all 
show heavy contractile stresses, can all exist in 
another crystalline form. The beta-modification 
of chromium, for example, can be formed from 
a chromic acid electrolyte if more than 18 percent 
trivalent chromium exists in the solution, and 
this modification reverts readily to the normal 
alpha-form. It is possible that the concentration 
of trivalent chromium may reach 18 percent in the 
cathode film by reduction of the hexavalent ion by 
liberated hydrogen. The beta-form is, however, 
more close packed than the alpha-modification, 
and a transformation should cause an expansive 
rather than a contractile stress. However, the 
same is true for nickel where the normal beta- 
nickel is a face-centered cube and the alpha-form 
a close-packed hexagonal structure. Nevertheless 
the density of the close-packed nickel when 
formed by sputtering in hydrogen was estimated 
as 7.04 against the normal value of 8.9.% Here 





1937) Brenner, J. Research Nat. Bur. Stand. 18, 565 
a L. Wright, H. Hirst, and J. Riley, Trans. Faraday Soc. 
31, 1253 (1935). 
*G. Bredig and R. Allolio, Zeits. f. physik. Chemie 
126, 41 (1927). 








a change from alpha- to beta-form would cause 
a contraction of a metal film. 

The data of Marie and Thon, given above for 
copper deposition, show that the stress in the 
deposits changes from expansive to contractile as 
the over-voltage of the cathode is lowered or the 
acidity of the electrolyte is increased. The change 
from expansive to contractile stress found in 
silver deposits on platinum® formed from the 
same solution of potassium argento-cyanide 
takes place with increase of current density. It 
would appear that in these cases there are two 
processes at work, one endeavoring to give 
expansive stress and one contractile. Under 
conditions giving rise to the liberation of a large 
quantity of hydrogen the second effect prevails 
over the first, and contractile stresses result. 

It is now suggested that in almost all deposited 
metals there is a tendency to form expansive 
stress, but in the case of transitional’ metals this 
tendency is always overcome, and in the case of 
borderline transitional metals it may be over- 
come by the formation of contractile stress 
resulting from interference by hydrogen in the 
growing lattice. Only in the case of non-transi- 
tional metals such as lead, zinc, and cadmium is 
the expansive tendency always greater than the 
contractile. The proposed mechanism for the 
formation of the expansive stress is based on the 
observation that in every stressed deposit investi- 
gated the AF° value for the reduction to metal of 
the metal ion in solution by monatomic hydrogen 
is favorable to the reaction taking place. Thus 
the AF° values are —61524 calories for zinc, 
—90820 for lead, and —112610 for copper for 
reduction from their divalent ions, and —42120 
calories for the reduction of silver from the 
Ag(CN)2- complex.*® Iron, nickel, cobalt, and 
palladium theoretically may also be reduced. 
These AF° values, of course, only indicate that 
the reduction is possible, since the precise con- 
centration of ions and hydrogen in the cathode 
film cannot be calculated. Assuming, however, 
that reduction of the metal ion by monatomic 
hydrogen takes place, metal atoms will be pro- 
duced at the cathode by two separate and 
distinct electron-transfer mechanisms. It is sug- 
gested that the metal formed by the hydrogen 


% Values calculated from W. M. Latimer, Oxidation 
Potentials (Prentice-Hall, Inc., New York, 1938). 











reduction has little tendency to continue the 
lattice structure of the nuclei resulting from the 
deposition of metal by electron gain from the 
cathode, and that it tends to grow between 
these nuclei. What may be envisaged crudely as 
a shouldering apart of the growing deposit occurs, 
and expansive stress ensues. Under conditions 
where more monatomic hydrogen than can 
reduce the metallic ions is being formed, con- 
tractile stresses will occur if the metal being 
deposited is capable fof absorbing hydrogen. 
Thus with silver at 5 milliamp. per square 
centimeter the stress is expansive ; at 15 milliamp. 
per square centimeter it is contractile. With 
copper contractile stress occurs only when hydro- 
gen is more readily evolved. The inability of 
non-transitional metals to form systems with 
hydrogen causes the stress in electro-deposits of 
these metals always to be expansive. 


VI. LIMITATIONS OF THE STRESS-ORIENTA- 
TION THEORY 


The theory proposed satisfactorily fits the 
data for chromium. Experimental findings are 
limited in the case of other metals, and insuffi- 
cient data exist to permit a detailed correlation 
of stress, orientation, and physical properties. 
For nickel deposits, however, this lack of data is 
not entirely true. Many papers on nickel deposits 
have been published, but because of the wide 
number of electrolyte compositions used and the 
complicating variable of pH in addition to 
current density and temperature, it is not easy 
to generalize the results. It would appear, how- 
ever, that the internal cold-working mechanism 
suggested to account largely for the hardness of 
chromium (supplemented by extremely small 
grain size and possibly by the effects of a residual 
interstitial solid solution of hydrogen), can 
account only for part of the hardness of many 
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deposits of nickel. It is possible that the major 
portion of the hardness of these deposits results 
from very finely dispersed basic material con- 
tained in the lattice. The hardness of all very 
hard electro-deposits is greater than that which 
could be engendered by cold-working the same 
massive metal or by the known effects of precipi- 
tation hardening, alloy formation, or diminution 
in particle size. It is possible that internal cold- 
working taking place within small grains immedi- 
ately after their deposition might produce greater 
effects than those resulting from ordinary cold- 
working processes, but on the whole it is felt 
that the anomalous hardness figures for electro- 
deposits result from a combination of all the 
hardening factors outlined above. The degree to 
which each contributes may vary from metal to 
metal, i.e., primarily cold-work in chromium and 
precipitation hardening by co-deposited basic 
material in nickel deposits. 

It is not suggested that growth mechanisms 
can never be responsible for orientations in elec- 
tro-deposits since some experiments***’ clearly 
show that, particularly at low current densities, 
preferential growth of distinct crystal planes 
may take place. It is, however, suggested that 
orientations in hard, fine-grained electro-deposits 
are the result of an internal cold-working mecha- 
nism occurring during deposition. 
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O the general reader the rule concluding 
this article will be of more interest than 
its derivation. It is slightly different from 
Koenigs original statement! and is perhaps more 
easily grasped. The major part of this article is, 
however, concerned with the group-theoretical 


and physical foundation of this rule.? 

By way of introduction, consider Maxwell’s 
equations for the electromagnetic field. This set 
of well-known equations is invariant for: 

a. simultaneous changes in sign of E, H, D, 

and B, 
b. simultaneous interchanges (E, H), (D, B), 
and (€, u), 

c. other transformations, not to be considered 

here. 


The physical interpretation of a is the rather 
trivial fact that the positive sense of all vectors 
may be reversed; of more importance is the 
meaning of 6 that nature is symmetrical with 
regard to electricity and magnetism, as far as the 
general field equations for insulators go (not 
meaning, of course, that for a definite substance 
e(Z) should always equal yu(A)). 

In much the same way a set of transformations 
exists that leaves the set of all general thermo- 
dynamic relations invariant. If for succinctness 
the argument is confined to a single homogeneous 
substance, and the Nernst law is disregarded, 
all general formulae may be deduced from the 
fundamental set: 


energy =U, dU=TdS—Pdv, 
enthalpy U+PV=W, dW=TdS+ VaP, 
freeenergy U-—-TS =F, dF =—SdT-—Pdv, 
fr.enth. U+PV—TS=G, dG =—SdT+VdP. 
(1) 


The symmetry of this set with respect to changes 
of sign of the “primary” variables P, V, 7, S and 


1F, O. Koenig, J. Chem. Phys. 3, 2° (1935); Professor M. 
Born informs me that the subject was also treated in his 
lectures at Géttingen and that a poem was made for memo- 
rizing the substitutions. 

?J. A. Prins, Physica 13, 417 (1947). 


(Received September 2, 1947) 





the “secondary” variables or characteristic 
functions U, W, F, G is summarized in (2) and 
(3). Each horizontal line of simultaneous changes 
is allowed: 


PVTS UG WF 
E ++ ++ ++ ++ 
A -- +4 ++ +41 wy 
Bo o++-- ++ +4’ 
c -- -- ++ ++ 
a -+ —+ -- -- 
b +- 4- -- -- 
a -+ + I al 
6B +--+ -- -- 


Physically, this table expresses little more than 
that the characteristic functions depend on 
products of “intensive” variables (P,7) and 
“extensive” variables (V, S). Another symmetry 
with respect to the interchanging of variables is 
again much more important. It is approximately 
summarized by the next set of eight substi- 
tutions (including unity E=no change): 





E P V T S UG WF) 

I -V P T S WF UG 
IP V-S T FWGU 
HI -V P -S T GU FUW\ 
IV S T V PGU WFf- 
VS T P-V FW UG 

VI T-S V PWF GU 
Vil T-S P-V UG FW 


It should be noted that “conjugated pairs” 
of variables keep together, not only (P, V) and 
(T, S), but also (U, G) and (W, F). Physically, 
this set expresses that in general thermodynamic 
relations (disregarding signs) no distinction is 
present between mechanical and thermal vari- 
ables nor between intensive and extensive vari- 
ables, a rather unexpected fact. 

A not very important but unpleasant com- 
plication is offered by the fact that the set (4) 
does not form a group and can only be made to 
do so either by disregarding the signs (which is 

































Fic. 1, Modification of Koenig’s thermodynamic square. 
The seven principal thermodynamic substitutions are 
found by rotating the square over 180° around its digonal 
axis and over 90°, 180°, and 270° around its tetragonal axis. 


inadmissible) or extending the number of trans- 
formations from eight to 32 (which is cumber- 
some) by multiplication of (4) by (2). One may 
even raise it to 64 by multiplying (4) by (2) and 
(3), or by (2) and a. 

In the total group of 64 transformations the 
32 form a maximal invariant sub-group. Its 
factor group is C2. In the group of 32 the only 
maximal invariant sub-group is (2), which is 
isomorphous with Klein’s ‘‘Vierergruppe,” V4. 
Its factor group is found to be Klein’s dihedral 
group for the square, which is usually called D, 
but will here be called Dg to indicate the number 
of its elements. So the composition of the total 
thermodynamic substitution group* of 64 ele- 


3 The group considered in this article, Gs, may be repre- 
sented in a totally different way by the 64 linear trans- 
formations (matrices), working on the two complex vari- 
ables P+iV and T+7S: 


a 0 0 a 
0 3 and Cl, ol 


where c stands for either taking the complex conjugate or 
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ments is Ggs=C2V4Ds. The unpleasant fact 
remains that D, itself does not occur as a sub- 
group in the 64 or 32 substitutions, though, on 
the other hand, its hidden presence lends the 
most important feature to the group, viz. the 
interchanges of variables. 

To sail between the scylla of non-group 
character and the charybdis of too many ele- 
ments, Koenig! and McKay,‘ although they do 
not state the dilemma so explicitly, have resorted 
to different devices. With a view to simplicity | 
suggest, however, the following 

Rule: Jn Fig. 1 a square 1s in its normal position 
with T at the top and P to the right, the conjugate 
variables being opposite and the characteristic 
functions occurring between those variables wherein 
they are normally expressed, as shown for instance 
by the differentials in (1). The operations indicated 
by the four diagonal axes and the one tetragonal 
axts transform the square into itself and in doing 
so substitute the letters in such a way, (4), that all 
general thermodynamic relations remain true after 
transformation. For the extensive variables V and 
S the additional rule must be observed that they 
change sign in passing the dotted line. 

As an example the reader may apply the 
transformations to one of Maxwell's thermo- 
dynamic equations, getting the other ones, or to 
the formula for the Joule effect, getting that for 
the Joule-Kelvin effect: 

P—T(dP/dT)y 
(@T/dV)v= - 


V 





—V+T(dV/dT)p 
(0 T/dP) ww C ’ 
P 





C standing for T0.S/dT. 


not, and independently a and 6 stand for either 1, —1, 2, 
or —i. 

The composition given in the text, Ges=C2V4Dsz, is not 
the only one possible. By taking (2)+(3) as first maximal 
invariant sub-group the composition DsC2V, is found. By 
Hélders theorem the two compositions differ, of course, 
only in the sequence of the factors. 
4H. A. C. McKay, J. Chem. Phys. 3, 715 (1935). 
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The infra-red absorption of gaseous hexafluoroethane between 2 and 25y has been inves- 
tigated with prism spectrometers of high resolution. On the basis of normal vibration calcula- 
tions and the known Raman frequencies it has been possible to interpret all but 3 of the 51 
observed absorption maxima. Only combination bands permitted for both the staggered (Dsa) 
and the eclipsed (Ds,) model were used in the assignment of bands. Five fundamental fre- 
quencies have been determined in addition to the six obtained from the Raman spectrum. No 
combination band involving the remaining frequency, associated with the torsional oscillation 
of one CF; group relative to the other, could be identified. An attempt to determine the 
moments of inertia of the C2Fs molecule from the envelope of the strongest parallel band gave 


unsatisfactory results. 





INTRODUCTION 


NUMBER of fluorocarbons have recently 

been prepared.' Their unique physical and 
chemical properties make the investigation of 
their infra-red and Raman spectra of consider- 
able theoretical interest, as well as of practical 
value in providing data for analytical applica- 
tions. Because of the structural analogy between 
fluorocarbons and hydrocarbons it may be ex- 
pected, moreover, that a study of the spectra of 
the former will throw light upon the vibrational 
problems presented by compounds of the latter 
important class. 

The Raman and infra-red spectra of tetra- 
fluoromethane, which has been available for some 
time, have been studied by several workers.*? 
P. Torkington and H. W. Thompson‘ have inves- 
tigated the infra-red spectrum of tetrafluoro- 
ethylene. Recently, Rank and Pace® have ob- 
served the Raman spectrum of liquid hexa- 
fluoroethane. In the present paper a study of the 

* Presented at the Symposium on Molecular Structure 
and Spectroscopy held at Ohio State University, Columbus, 
Ohio, June 9 to 14, 1947. 

1Compare March 1947 issues of Ind. and Eng. Chem. 
and Anal. Chem. 

2D. M. Yost, E. N. Lassettre, and S. T. Grosse, J. 
Chem. Phys. 4, 325 (1936). 

3C. R. Bailey, J. B. Hale, and J. W. Thompson, Proc. 
Roy. Soc. London 167A, 555 (1938). 

*P. Torkington and H. W. Thompson, Trans. Faraday 
Soc. 41, 236 (1945). 


5D. H. Rank and E. L. Pace, J. Chem. Phys. 15, 39 
(1947). 
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infra-red absorption spectrum of gaseous hexa- 
fluoroethane will be reported. 


EXPERIMENTAL 


The sample of hexafluoroethane was prepared 
and purified in the Fluorine Laboratory of the 
Pennsylvania State College and was kindly 
made available by Professor J. H. Simons. It 
was part of a sample on which thermodynamic 
measurements have been made by Aston and 
Pace.® Its purity, as estimated by these workers, 
was 99.6 mole percent. The small amount of the 
impurity was found by Aston and Pace to be 
low-boiling, and a subsequent distillation raised 
the purity to 99.9 mole percent. A part of this 
latter sample was used by Rank and Pace for 
obtaining the Raman spectrum. 

The infra-red absorption spectrum was re- 
corded, from 2 to 154, by means of a spectrometer 
similar to one described by Wright.? The 60° 
rocksalt prism, in Littrow mounting, has faces 
10 by 8 cm and the collimating parabolic mirror 
a focal length of 91.4 cm. A Pfund-type com- 
pensated vacuum thermocouple serves as re- 
ceiver, and the galvanometer deflections are 
amplified by a balanced barrier-layer photo-cell 
device and recorded on bromide paper whose 


6]. G. Aston and E. L. Pace, to appear shortly in J. 
Am. Chem. Soc. 
7™N. Wright, Ind. Eng. Chem., Anal. Ed. 13, 1 (1941). 

















68 J. R. NIELSEN, C. M. 





motion is synchronized with that of the Littrow 
mirror. The cell length was 43 cm. Several 
pressures up to 194 mm of mercury were used, the 
temperature of the gas being about 30°C. 

Through the courtesy of Mr. D. C. Smith of 
the Naval Research Laboratory, a few records 
were made with a spectrometer of higher resolv- 
ing power.’ In particular, the absorption in the 
region from 15 to 25u was determined with this 
instrument. A 10-cm cell was used with the 
NRL spectrometer. The pressure ranged from 
less than one millimeter to 800 mm. 


RESULTS 


The infra-red absorption of hexafluoroethane 
gas from 2 to 15y is shown in Fig. 1 which gives 
complete percent transmission curves on a uni- 
form wave-length scale for three different 
pressures, as well as partial transmission curves 
for the two strongest bands obtained at a 
pressure of less than one millimeter. The curves 
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are based on the records made with the Univer- 
sity of Oklahoma spectrometer but have been 
checked against the data obtained with the 
NRL instrument and have been slightly im- 
proved in a few places. The scale of the drawing 
is not quite large enough to show all of the 
details revealed by the NRL spectrometer. 
Faint bands, at 2890, 2967, and 3060 cm™, 
appeared on the first records obtained with 
hexafluoroethane at pressures of 6 and 48 mm 
but were absent, or at least greatly reduced in 
intensity, on the last records obtained at a 
pressure of 194 mm. Before the last records were 
made the material was frozen several times and 
the remaining vapor removed by pumping. It 
appears, therefore, that the bands mentioned are 
caused by a volatile impurity. Their wave 
numbers lie close to the Raman frequencies, 
2885, 2974, and 3040 cm~, observed by Hatcher 
and Yost? for liquid CH3CF3. The band of highest 
frequency nearly coincides with the intense 
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Fic. 1. Infra-red absorption spectrum of hexafluoroethane gas from 2 to 154. Percent transmission is plotted against 
wave-length on a uniform wave-length scale. Cell length: 43 cm. Pressures: 6, 48, and 194 mm. Temperature 30°C. 


8 J. Rud Nielsen, F. W. Crawford, and D. C. Smith, J. Opt. Soc. Am. 37, 296 (1947). 
®W.H. Hatcher and D. M. Yost, J. Chem. Phys. 5, 992 (1937). 
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Raman band, 3062.0 cm found by Glockler and 
Edgell?® for liquid fluoroform. When the cell was 
being evacuated after the first filling, a weak band 
was observed at 1267 cm~, close to the most 
intense hexafluoroethane band, 1250.5 cm. 
After the material had been frozen several times, 
the ratio between the intensities of 1267 cm— and 
1250.5 cm~ had fallen off very greatly. It was 
concluded, therefore, that the former band 
belonged to a volatile impurity. Since Bailey, 
Hale, and Thompson* have found the most 
intense band of CF, to lie at 1265 cm, the 
impurity causing the 1267 cm band is evidently 
tetrafluoromethane. 

Beyond 15y, only three bands were observed: 
a strong diffuse. band with maximum at 522.5 
cm~!, a weak band with Q-branch at 632.5 cm—, 
and P- and R-branches at 613 and 649 cm=!, and 
a very faint band at 599 cm—. The band at 632.5 
cm was at first not regarded as an impurity 
band (but was interpreted as v7—vu, although 
with misgivings because of what seemed to be an 
abnormally large separation of the rotational 
branches). However, when it was noted that this 
band agrees closely, both as regards frequency 
and intensity distribution, with the second 
strongest infra-red absorption band found by 
Bailey, Hale, and Thompson’ for tetrafluoro- 
methane, and when the 1267 cm~ band had been 
interpreted as caused by this impurity, it 
became clear that the 632.5 cm~! band must also 
be ascribed to the presence of a small amount of 
CF. 

The wave-lengths of the observed absorption 
maxima (or “shoulders’’), after elimination of 
the impurity bands mentioned, are listed in the 
first column of Table I and the corresponding 
wave numbers in the fourth column. Rough 
indications of the intensities and shapes of the 
bands are given in the second and third columns 
under the joint heading ‘‘Description.”’ The fol- 
lowing abbreviations are used: vs extremely 
strong, vs very strong, s strong, m medium, w 
weak, vw very weak, vw extremely weak, vw 
barely observable, d diffuse, and sh sharp. In 
many cases, especially at shorter wave-lengths, 
it is not possible to designate bands as either 
diffuse or sharp. 


10 G,. Glockler and W. F. Edgell, J. Chem. Phys. 9, 224 
(1941). 
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TABLE I. Infra-red absorption spectrum of 
hexafluoroethane gas. 




















Wave- Descrip- Calculated 
length tion Frequency Assignment frequency 
19.14 s d 522.55cm> | vs — 
16.7 vw d 599 |] 1 vo vue 596 cm"! 
14.13 vs d 707.3 P 
14.005 vs sh 714.03 | vs — 
13.89 vs d 720.2 R 
13.04 veew 766.9 | Va 768 
12.78 vovw 782.5 
11.88 vow d 842 in veri 836 
11.52 vow d 868 | Lo ovr—rie 871 
11.38 vovw d 879 L sonatos 872 
11.19 vw a 894 Va vi—vs 897 
11.08 m 903 |} Ls sovs+ri2 903 
10.56 vw 947 
9.756 m d 1025.0 . a vetvs os 
9.50 vow 1053 P 
9.410 vw sh 1062.7 || vitve 1063 
9.35 vow d 1070 R 
9.01 s dad 1110 P 
8.953 ws sh 1116.9 | vs ao 
8.89 ts 1125 R 
8.771 m sh 1140.1 ie oT 1143 
8.75 vw da 1143 R 
8.714 m 1147.6 || 1 v6 +29 1146 
8.298 m 1205.1 LL sn—s 1204 
7.997 ws d 1250.5 Low - = 
7.799 om 1282.2 | L (va+-2012 +02) —vs 
7.787 m 1284.2 | Lvs +2012 1283 
7.528 m 4d 1328.4 Ls vets 1332 
6.92 vow 1445 P 
6.882 m sh 1453.1 | L vo+y10 1453 
6.83 vw 1464 vs+vor R 1466 
6.57 ww sh 1522 | vetve 1523 
6.246 m 1601.0 i. v3+r7 1600 
6.238 m 1603.1 Lo (ator tey) —ry 
6.130 m sh 1631.3 |] 1 vr +oie 1631 
6.10 vuvw 1639 5. vitvs or R 1636 
5.770 m d_ 1733.1 Louse 1737 
5.64 vrw 1773 | L vs+vi0 1760 
5.357 m 1866.7 |} Lo wrt 1871 
5.198 vow 1924 | vatve 1926 
5.151 m 1941.4 vitvs 1943 
5.115 ww d 1955 L vet 1951 
4.867 w 2055 a voter 2060 
4.69 vow 2132 | ntvs 2134 
4.572 vw 2187 
4.44 vere 2252 |] Lorton tris 2251 
4.252 m 2352 Lovet 2354 
4.033 m 2480 | Ls ovr-+tvw 2488 
4.01 m 2494 Lovrt2eu 2491 
3.734 om 2678 Loner 2671 
3.64 vw d 2750 2a vitvatvs 2752 
INTERPRETATION 


Two equilibrium configurations must be con- 
sidered for the hexafluoroethane molecule, the 
“staggered” model of symmetry D34¢ and the 
“eclipsed” form of symmetry Ds,. The sym- 
metry species and selection rules for the twelve 
fundamental vibration frequencies are shown, for 
either of these models, in Table II. The number- 
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FABLE II. Normal vibrations of ethane-like molecules. 











Dsa D3 

Designation Species Activity Species Activity 

m1 Aig R-p Ai’ R-p 

v2 Aig R-p Ay’ R-p 

V3 Aig R-p Ay’ R-p 

V4 Aw inactive A,” inactive 

V5 Aou I | A," I | 

V6 Ax I A,” I 

vy ae E’ I 1, R-dp 

vs E, IL E’ I 1, R-dp 

V9 me £2 E’ I 1, R-dp 

v0 E,  R-dp E"” ~— R-dp 

ru E,  R-dp E” R-dp 

Vi2 E, R-dp ze” R-dp 
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TABLE III. Calculated values of fundamental frequencies 
of the hexafluoroethane molecule. 











Nota- No F-F F-F repulsion; F-F repulsion; 
tion repulsions r-c-c =kp-c-F ky-c-c >ky-c-F 

V1 1734 cm™ 1392 cm 1510 cm™ 

v2 824 799 813 

V3 346 209 301 

V4 — — oe 

V5 1298 875 940 

V6 722 505 - 742 

Vi 1710 1214 1223 

Vs 566 515 ES | 

V9 316 134 316 

V10 anee ae jo 


Vil wn a 


Vie ee a 








ing of the fundamentals is that used by Herzberg" 
for ethane. J means active in infra-red absorp- 
tion, R active in the Raman effect, || denotes a 
parallel and | a perpendicular infra-red band; 
p indicates a polarized and dp a depolarized 
Raman band. 

As an aid in identifying the observed infra-red 
frequencies, we have computed eight of the 
fundamentals, following essentially Howard’s 
procedure for ethane,” except that repulsive 
forces between F atoms within each CFs; group 
were assumed in setting up two of the three 
potential functions used. The C—C stretching 
force constant was that given by Howard. The 
three remaining force constants were obtained by 
means of Urey and Bradley’s equations" which 
give the CX, frequencies when X—X repulsion 
exist. The symmetry D3, was assumed but the 
results would have been the same for D3, or, 
more generally, for Howard’s D3,’. The results 
are given in Table III. In the first column are 
listed the frequencies obtained when no repulsive 
forces are assumed to act between the fluorine 
atoms, in the second column are given the fre- 
quencies obtained when such repulsion is 
assumed to exist and all angle constants are 
taken to be equal, and in the last column are 
listed the frequencies computed under the as- 
sumption that repulsion between the fluorine 
atoms exists and that the C-C-F angle constant 
is different from the constant for the F-C-F 


1G, Herzberg, Infrared and Raman Spectra of Poly- 
atomic Molecules (D. Van Nostrand Company, Inc., New 
York, 1945). 

12 J. B. Howard, J. Chem. Phys. 5, 442 (1937). 
193 Dy, C. Urey and C. A, Bradley, Phys. Rev. 38, 1969 





angles. Since these calculations were made we 
have solved the normal vibration problem more 
completely. The results will be published else- 
where. 

In spite of the uncertainty as to which of 
three sets of computed frequencies is most nearly 
correct, comparison with the observed infra-red 
and Raman frequencies gives with certainty all 
but two of the twelve fundamentals. The inter- 
pretation of the Raman spectrum has already 
been given by Rank and Pace.® The infra-red 
spectrum of hexafluoroethane contains two very 
intense bands, at 714.0 and 1116.9 cm~, which 
clearly exhibit parallel structure. They must be 
vg and ys, respectively. Another very intense 
band, at 1250.5 cm™, appears from the record 
obtained at the lowest pressure to be a per- 
pendicular band. It must be interpreted as 7. 
The rather strong diffuse band at 522.5 cm 
must be 7s. 

Great effort was spent in determining v9. The 
fact that no appreciable absorption is observed 
near the 25y limit of the potassium bromide 
region shows that vg must lie well below 400 
cm-!. The value of vy which leads to the most 
satisfactory interpretation of the infra-red spec- 
trum was obtained in the following manner. The 
observed spectrum was first interpreted as com- 
pletely as possible in terms of the ten funda- 
mental frequencies, 1, v2, v3, ¥5, V6, V7, Vs, V10 
11, and v32. When that was done only four bands 
of medium intensity remained unassigned. If one 
of these, 1025 cm=, is assumed to be v2+ 79, 
making vy=216 cm™, each of the other bands 
finds either an entirely or a fairly satisfactory 
interpretation, Several other assumptions, lead- 
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ing to different values for vs, have been con- 
sidered and discarded as less plausible. 

All attempts to identify combination bands 
involving vs, and thus to determine this fre- 
quency, failed. The values for the other funda- 
mental frequencies are listed in Table IV. The 
spectrum from which the frequency is deter- 
mined is listed in the third column. The funda- 
mental frequencies obtained from the Raman 
spectrum of liquid hexafluoroethane are probably 
somewhat less accurate than the frequencies 
determined from the infra-red absorption spec- 
trum of the gas. The value for v9 is considerably 
less certain than those found for the other funda- 
mentals. 

On the basis of the fundamental frequencies 
listed in Table IV it has been possible to give a 
satisfactory interpretation of all of the 51 ob- 
served infra-red absorption maxima with the 
exception of three very faint bands. The assign- 
ment of each band is given in the fifth column 
of Table I. The calculated frequency is given 
in the last column. The agreement between the 
observed and the calculated frequency is very 
close for practically all of the bands. Only for 
four bands is the difference greater than 5 cm 
and only for one band does it exceed 10 cm™. 

The designations ||, 1, and ||1, for parallel, 
perpendicular, and parallel-perpendicular, bands 
given in Table I, assume the molecular equilib- 
rium configuration to have the symmetry Dg3a. 
If the molecular symmetry were D3, some com- 
bination bands, such as V7, 8,92 V10, 11, 12 and 
¥5,6-2v7,8,9, would be || instead of ||. Other 
bands, such as v7,s,9-+2vs,9,7, would be instead 
of || 1. 


TABLE IV. Fundamental vibration frequencies 
of hexafluoroethane. 
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DISCUSSION 


Rank and Pace‘ observed a faint Raman band 
at 655 cm and suggested that it might be 2g, 
making »vy=328 cm~. This assumption explains 
the very weak infra-red band at 947 cm= as 
vgt+vy; and the extremely weak band at 2187 
cm~! as vg+vio+711. However, it leads to the not 
very plausible interpretation of the band 1025 
cm as 2v3+v9 and leaves the bands 1147.6, 
1205.1, and 1453.1 cm unassigned. For these 
reasons, Rank and Pace’s interpretation of the 
faint Raman band at 655 cm=! seems unac- 
ceptable. However, we have not been able to 
interpret this Raman band in any other plausible 
manner. 

Four bands are interpreted as difference bands, 
Ya—v. In all but one case the corresponding 
band, »v,+™, has been observed with about the 
intensity to be expected from the intensity of 
¥a—v and the Boltzmann factor for the excited 
initial state ». The exception is 1205.1 cm 
which is interpreted as v;—¥v9. The absorption 
maximum at 1639 cm™ is probably »1+ 9, but 
its intensity is abnormally low. It is possible 
that the intensity of the band 1205.1 cm has 
been enhanced because of the proximity of 7. 

The splitting of the bands at 1282-1284 and 
1601-1603 cm is interpreted by the assumption 
that the weaker component in either case is 
caused by a transition from an excited vibra- 
tional level to an upper level in which this 
vibration remains excited.“ It is not possible to 
determine the initial excited levels v, and vy. 

Two cases where the assigned band may be 
merged with a rotational branch of a stronger 
band are indicated in Table. I. A number of 


TABLE V. Observed combinations va+ 7. 

















Spectrum in which Eu*Eg =A. 

Notation Frequency observ: Aou*Aig=A2u Aou*Eg =Eu Eu*Aig=Eu +Aou+Eu 

V1 1420 cm R (liq) Vetve view Vstvio m vrtv1 m vitvie m 

v2 809 R (liq) Vet V3 VW vst Vil om vit Yo W vr+ Vil ™m 

V3 349 R (liq) Va—V3 VUW vstvig VW vytvs mM vitvie m 

V%4 —_ —_— vet vvw Vet vio vw vet m VI—Vi2e VIUW 

V5 1116.9 I (g) vetve vw vetvis vi-vg UW vetvig vevw 

V6 714.0 I (g) Vetvs vw vetve m vetvi mM 

v7 1250.5 I (g) vetv3s vovw vetrie m 

Vs 522.5 I (g) vitveg vow votvio m 

Vg 216 I (gz) Vi Vg.—_—_ M vot 2 ee es 

Vi0 1237 R (liq) votv, m votrie vw 

Vil 620 R (liq) 

V2 380 R (liq) 








14 Compare reference 11, p. 267. 
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NIELSEN, C. M. RICHARDS, AND H. L. McMURRY 


Temperature 28° 
re Pressure 8mm 


Cell Length 43 cm 


Fic. 2. Intensity distribution 
in the parallel band at 714.0 
cm~, The extinction, log(Jo/Z), 
obtained with hexafluoroethane 
at 6 mm and 30°C in the 43-cm 
cell is plotted against angular 
setting of Littrow mirror. 





695700 +705. 710 715 720 725 


Wave Numbers (cm™) 


other bands, especially some of high wave 
numbers, can be interpreted in more than one 
way. Thus, the band at 902.5 cm may also be 
interpreted as v7—v3; 1328.4 cm~! may be inter- 
preted as vg+vi1, 1866.7 as vo+v3+6, and 2132 
as 2v2+ vs. At least two alternative assignments, 
Yvst2v11 and vs+ve+yvs, are possible for 2352 
cm~, 2480 cm! may also be interpreted as 
vitvstre and 2494 cm- as v3+3y¢. In general, 
the interpretation chosen as the most plausible 
and given in Table I is that for which the total 
change in vibrational quantum numbers is least. 

The combinations used in the assignment are 
all such which are active in infra-red absorption 
for both equilibrium configurations, D3g and D3,, 
of the hexafluoroethane molecule. The observed 
combination bands »,+», may be divided into 
four groups according to the symmetry species 
of vy, and ». This is shown in Table V which also 
gives the intensities of the observed bands, 
except in the case of vs+v11 which overlaps with 
vet+vs. It will be noticed that all of the com- 
bination bands of species Aey (for Da) are very 
faint, whereas most of the other combination 
bands are of higher intensity. 

It is a remarkable fact that only three com- 
bination bands of the types listed in Table V 
were not observed. The band »3+¥7;, which may 
be expected to be very faint, should lie in the 
short-wave region where identification is dif- 
ficult. The band v9+ v3, on the other hand, should 
lie"in the potassium bromide region where the 
presence of a faint band may have escaped ob- 
servation, since only one record was made. The 





730 735 


third missing band, vg+v12, may be masked by 
vs. It is surprising also that it was not found 
necessary to assume the existence of Fermi 
coupling between any of the excited vibrational 
states. 

The fact that no overtone band of the type 
2v, was observed may be taken as evidence in 
favor of the staggered model of symmetry D3a; 
for, while all such overtones are forbidden for 
Dsa, overtones of type 2», belonging to the 
species E’, i.e., 2vz, 2vs, and 2y9, can occur if the 
equilibrium configuration has the symmetry Ds;. 
There is no indication of a band in the long-wave 
region where 2v9 should occur. The absence of 
2vg cannot be claimed with certainty, since it 
would nearly coincide with the P-branch of the 
parallel band v3+ 5. 2v7 would lie approximately 
at 2501 cm—. A band found at 2494 cm is 
interpreted in Table I as v7+2v1:. The possi- 
bility that this band is actually 2v7 cannot be 
entirely excluded ; for, if the frequency v4, cor- 
responding to a torsional oscillation of one CF; 
group relative to the other, is less than RkT/h, 
a hindered rotation may occur and destroy the 
rigorous distinction between the molecular sym- 
metries Dsz and D3,. Six combination bands 
V7,8,9+ V3, 9,7 and V0, 11,12 +11, 12,19 are permitted 
if the symmetry of the molecule is D3, but are 
forbidden for Dza. The band v7+ gs, if present, 
should lie at 1773 cm—". An extremely faint band 
found at this wave number has been interpreted 
in Table I as vg+yv12, which should occur at 
1760. vz-+v9 should coincide with the parallel 
band vs+v3. The band vg+y9 should lie far 
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enough from the R-branch of vg to be observable. 
However, there is no indication of its presence. 
Of the bands of the second type, vio+711 should 
lie at 1857 cm=', which is close to »7+v7, at 
1866.7 cm. There is no indication of its 
presence. The bands y+ 712 and »31:+ 712 should 
both lie sufficiently far from other bands to be 
observable, but neither was found. Thus, the 
evidence favors the conclusion that the sym- 
metry of the hexafluoroethane molecule is repre- 
sented better by the staggered than by the 
eclipsed model. 

The bands v4 V0, 11,12 and 42740, 11,12 
are active in infra-red absorption for both 
molecular symmetries. The bands +»4+27,¢9 
are permitted for the staggered form but for- 
bidden for D3,. For either symmetry of the 
equilibrium configuration, +24 added to any 
active fundamental or combination band will 
give a pair of permitted bands. A thorough 
search was made for bands of the types men- 
tioned, but no band involving v4 could be iden- 
tified with any degree of certainty. Evidently, 
transitions in which the quantum number V; 
changes have a very low probability. It may be 
of interest to recall that a few bands involving 
v, are believed to occur in the infra-red spectrum 
of ethane." However, in no case have both 
bands of a pair, corresponding to both signs for 
v4, been found. Hence, the assignment of these 
ethane bands seems open to doubt. 


MOMENTS OF INERTIA 


According to Gerhard and Dennison" it is 
possible to determine the moments of inertia of 
fairly heavy symmetrical top molecules from 
measurements of the ratio of the intensity of the 
Q-branch to the total band intensity and the PR 
separation in a parallel band. We have tried to 
apply their theory to the 714 cm hexafluoro- 
ethane band. The extinction, logI)/J, as obtained 
with hexafluoroethane at a pressure of 6 mm and 
a temperature of 30°C in the 43-cm cell, is 


18S. L. Gerhard and D. M. Dennison, Phys. Rev. 43, 
197 (1933). 


plotted against the angular setting of the 
Littrow mirror in Fig. 2. No correction was 
made for stray radiation. The ratio between the 
intensity of the Q-branch and the total intensity 
of the band was estimated to be 0.17 by meas- 
uring with a planimeter the areas under the 
envelopes shown. The separation between the 
extinction maxima in the P- and R-branches 
is 13.1 cm. From Gerhard and Dennison’s 
curves B=A/C-—1 is found to be 3.1 and 
S=rAv/(kT/A)* is found to be 1.2. Solving for 
the molecular moment of inertia, A, about an 
axis perpendicular to the threefold symmetry 
axis and the moment of inertia, C, about the 
axis, we find A =401X10-” and C=98X10-*”. 

These results, especially the low value for C, 
are not very reasonable. If the distance between 
the carbon atoms is taken to be 1.54A, as in 
ethane, if the distance C—F is 1.36A as in tetra- 
fluoromethane, and if tetrahedral angles are 
assumed, the values A = 462.4X10-”, C=311.2 
X10-“, and 6=0.486 are obtained. If the C—C 
distance is smaller in hexafluoroethane than in 
ethane, as claimed by Brockway,'* A and 6 
would be somewhat smaller. Although no great 
accuracy can be claimed for the extinction curve 
given in Fig. 2, nor perhaps for Gerhard and 
Dennison’s procedure, the high value found for 
8 is nevertheless surprising. On the other hand, 
the perpendicular hexafluoroethane bands appear 
more different from the parallel bands than one 
would expect on the basis of Gerhard and Den- 
nison’s theory if 8 is as small as 0.486. 
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Vibrational potential functions of the quadratic type 
for CF, and C.F., which assume no interaction between 
the CF; groups in C2F¢ and include only the interaction 
terms involving bond-bending with bond-stretching 
motions, are used with the observed normal frequencies to 
evaluate a set of force constants for the molecules. The 
force constants associated with the CF; group are used 
interchangeably between the two molecules. The repro- 
duction of the observed values of the normal frequencies 


resulting from the use of the force constants is of the 
order of five percent or better except for three of the 
frequencies. The omission of terms in the potential function 
for the interaction of the CF; groups is largely responsible 
for the relatively poor reproduction of the type EZ, fre- 
quencies of C:F.. Additional interaction terms, e.g., 
between bond-bending motions, are necessary for the 
improvement of the agreement in. the case of the type F, 
frequencies of CF,. 





INTRODUCTION 
HE nature of the forces holding atoms 
together in molecules can be determined 
from a consideration of the force constants in a 
given molecule or in molecules with structurally 
similar groups. The force constants are calcu- 
lated from the observed normal frequencies and 
the assumed model for each molecule. The cal- 
culation depends fundamentally on the proper 
choice of a vibrational potential function. At 
present, the most promising potential function 
is the general quadratic type, limited to small 
vibrational amplitudes, which contains the bond- 
stretching and the bond-bending force constants 
together with the interaction force constants 
that are considered important for the particular 
molecule. In general, the validity and usefulness 
of the potential function rests on the repro- 
duction of observed normal frequencies in a 
greater number than the potential constants in- 
volved in it as well as in the agreement obtained 
between observed and calculated frequencies 
when the force constants for a given group in 
one molecule are used for the same group in a 
different molecule. For several molecules com- 
posed of carbon, hydrogen, deuterium, oxygen, 
and other atoms, excellent reproduction or pre- 
diction of observed frequencies has resulted from 
the transference between molecules of stretching, 
bending, and interaction force constants for the 
same group.'~* 
1H. W. Thompson and J. W. Linnett, J. Chem. Soc. 
1376 (1937). 
2B. L. Crawford, Jr. and S. R. Brinkley, Jr., J. Chem. 
Phys. 9, 69 (1941). 


3J. E. Kilpatrick and K. S. Pitzer, J. Research Nat. 
Bur. Stand. 38, 191 (1947). 
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The assignment of the normal frequencies of 
C.F, from the infra-red‘ and Raman spectra,° 
when combined with that for CF,,* now makes 
possible the consideration of the force field for 
molecules consisting solely of carbon and fluorine 
atoms. Because of the more complicated force 
field associated with the fluorine atom, the vibra- 
tional potential functions cannot be expected to 
be as simple as for the corresponding hydrogen 
or deuterium analog. Therefore, whereas a simple 
valence-force field is found to be adequate for 
CH, and CD,, it is unsatisfactory for CF4.’ 
Likewise, the potential functions employed by 
Howard? and Stitt® for C2H¢ and C2D¢ do not 
give real force constants when used with the 
observed normal frequencies for C.F. 

Some information concerning the significant 
interaction terms in the potential function for 
molecules containing fluorine can be derived 
from a consideration of the structurally analogous 
hydrocarbons. Hence, in C2Hg and C2Dg, the 
fact that the frequencies for either the two de- 
generate or the two non-degenerate types are 
essentially the same indicates that the C—H and 
C—D stretching motions are not coupled with 
any other motion. Furthermore, the comparison 
of the corresponding symmetrical frequencies in 


‘From the paper by J. R. Nielsen, C. M. Richards, 
and H. L. McMurry, presented at the Symposium on 
Molecular Structure and Spectroscopy, Columbus, Ohio, 
June, 1947. 

( 6 D H. Rank and E. L. Pace, J. Chem. Phys. 15, 39 
1947). 

6 G. Herzberg, Infrared and Raman Spectra of Polyatomu 
Molecules (D. Van Nostrand Company, Inc., New York, 
1945), p. 167. 

7 See reference 6, p. 181. 

8 J. B. Howard, J. Chem. Phys. 5, 442 (1937). 

°F, Stitt, J. Chem. Phys, 7, 297 (1939). 
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C:H¢ and C2Dg points to the absence of coupling 
between the H—C—H bending and the C—C 
stretching motions in C,H, but to the presence of 
considerable interaction between the D—C—D 
bending and the C—C stretching motions in 
CoDeg. 

Stitt® also suggests that the weak intensity or 
entire absence of the line for the symmetrical 
hydrogen bending frequency in the Raman spec- 
tra of several molecules, e.g., C2Hs, CHCl, and 
CH,;OH, may be caused by the polarizability of 
the molecule undergoing little or no change in 
the uncoupled hydrogen bending vibration. Con- 
sequently, the appearance of the line for this 
vibration with considerable intensity in the 
Raman spectra of molecules such as C2D¢. and 
CoH, may be attributed to the interaction of 
the symmetrical bending motion with another 
having a pronounced change in polarizability 
such as the C—C stretching motion. 

An examination of the data from the Raman 
and infra-red spectra of C.F, in the light of 
that for C2Hs and C2D¢ reveals two pertinent 
facts. First, the frequencies for both the C—F 
stretching and F—C—F bending vibrations are 
quite different for the two non-degenerate types. 
Secondly, the line in the Raman spectrum asso- 
ciated with the symmetrical bending frequency, 
vo, has a higher intensity than any other line.® 
Therefore, strong coupling probably exists be- 
tween the F—C—F bending and both the C—F 
and C—C stretching motions. 

In order to test the validity of the preceding 
conclusion, the succeeding calculation of force 
constants for CF, and C.F, from observed 
normal frequencies is performed with potential 
functions containing bond-stretching and bond- 
bending force constants, and all interaction force 
constants involving bond-bending with bond- 
stretching motions exclusive of those for the 
interaction of the two CF; groups in C2F5. The 
development and solution of the secular equation 
from the potential function and geometry of the 
molecule follows the method devised by Wilson." 


CALCULATIONS FOR CF, 


The form of the potential function for CF,, 
obtained from the preceding considerations, is 





_E. B. Wilson, Jr., J. Chem. Phys. 7, 1047 (1939); 
ibid. 9, 76 (1941). 
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im4 i= 
2V=fa dD (Adi)? +fa SY (dAai;)? 
i=1 i>j=1 


i=4 

+2 faa x (Ad;) (dAa;;), 
in which V is the vibrational potential energy, 
d, the C—F; equilibrium distance, a;;, the 
F;—C—F; equilibrium bond angle, and fu, fa, 
faa are, respectively, the C—F stretching con- 
stant, the F—C—F bending constant, and the 
constant for the interaction of these two motions. 

The derivation of the secular equation for the 
vibrations of a molecule of this type has been 
done by Meister and Cleveland." If their solution 
is modified by setting faa=0, the resulting 
equations become applicable to CF, with the 
preceding potential function. 

The CF, molecule has the symmetry 7,4. 
Therefore, the secular equation factors into one 
equation for the non-degenerate type A, vibra- 
tion, a second equation for the doubly degenerate 
type E vibration, and a third equation for the 
two triply degenerate type F2 vibrations. The 
two equations for the type A: and type E vibra- 
tions give values for fz and fa of 9.15 10° and 
0.71105 dynes per centimeter, respectively, if 
the frequency assignment given by Herzberg’ is 
used. However, in the equation for the F»2 fre- 
quencies, it is impossible to evaluate a real force 
constant, faa, which is consistent with the ob- 
served frequencies and the values of fa and fa 
calculated from the other type frequencies. 


CALCULATIONS FOR C.F, 


Following the procedure for CF,, the potential 
function of the quadratic type containing all 
interaction terms involving bond-stretching with 
bond-bending motions excepting those for the 
opposing CF; groups for an assumed Dzq con- 
figuration is 

i=3 


2V=fp(AD)?+fa > [(Ad,’)?+ (Ad,’’)? ] 


i=1 


i=3 


+fa D2 [(dAa;;')?+(dda;;’)*] 


i>j=1 


+fa © [(das,’)*+ (48,")*) 


ot (continued on next page) 


1 A. G. Meister and F. F. Cleveland, Am. J. Phys. 14 
13 (1946). 

















































t=3 


+ 2foa ; (AD) (dAa;;' +dAa;;'’) 


i>j=1 


t=3 


+2fop u (AD) (dAB,’ +d48,’’) 


i=3 
+2 faa >> ((Ad,’) (dAay, 41’ +dAai;, »-1’) 


in 
+ (Ad,’’) (da, 41 +dAc;, -1’’) ] 


t=3 


+2 fas > [(Ad;’) (d4B;’) + (Ad,”’) (d4;’) ], 


in which V is the vibrational potential energy ; 
D and d;, the equilibrium C—C and C—F;, bond 
distances; a;; and §;, the equilibrium F;—C—F; 
and F;—C—C bond angles. The single and 
double primes refer to different CF; groups. The 
force constants fp, fa, fa; fs, foa, fos, faa, and fag 
are defined by the bond and angle displacements 
associated with them in the above expression. 








PACE 


By making use of the symmetry of the mole- 
cule, the secular equation for the normal vibra- 
tions can be factored into equations for the E,, 
E., Aig, and As, symmetry types. Since these 
can be derived simply by the Wilson method” 
from the corresponding F and G matrices, the 
matrices will be presented in the following 
treatment rather than the expanded equations. 
The matrices are symmetrical with respect to 
the main diagonal. Hence, only the off-diagonal 
elements on the upper side of the main diagonal 
are shown. 

The G matrices, assuming a tetrahedral CF; 
group, are, for the type Ai, and type Ao, 
vibrations 


2uc —(6'/3)ue (4/3) (6')uce 
urtpuc/3 — (4/3)uce ; 
[2ur+(16/3)uc Je 


and for the £, and E, type vibrations 


ur+(4/3)uc  —(23/3)ucSe (4/3) (24)uce 
[urt+(uc/6)S Je  [ur/2—(2/3)ucS]e |, 
L(5/2)ur+(8/3)uc je 


in which yr and ye are, respectively, the recip- 
rocal masses of the fluorine and carbon atom; e, 
the reciprocal C—F equilibrium bond distance 
and S=1 for type E, vibrations and S=6(r/e)+1 
for the type E, vibrations, 7 being the reciprocal 
of the C—C bond distance. 

The F matrix for the type A1, vibrations is 


fo 9 33d(foa— fog) 
fa 2'd(faa—fag/2) |, 
(d?/2)(fa+fs) 


while that for the A2, vibrations is obtained 
from the preceding by setting fo=fpa—fps=9. 

The F matrix for both the type E, and type 
E,, frequencies is 


fa fag —Adfaa 
Pf, 0 
df 


The potential function and, hence, the secular 
equation make use of a total of eight different 
force constants, of which seven are independent 
since fpa and fps appear only in the combination 


foa— fp. 


In evaluating the most consistent set of force 








constants for C.F., the values of two of the 
force constants, f2=9.15X10° dynes per cen- 
timeter and f,=0.71 105 dynes per centimeter, 
are taken from the type A; and the type E 
frequencies for CFy. The remaining force con- 
stants are derived by a process of adjustment to 
the best fit between observed and calculated 
values for the type A1, and type E, frequencies 
for CFs. The values obtained, in units of 10° 
dyne per centimeter, are: fa=1.13, fag=1.35, 
f2=0.80, foa—fog= —0.28, and fo=5.A45. 

A C-—F distance of 1.35A and the observed 
value of 1.45A for the C—C distance!” were used 
in all calculations. The CF; group was assumed 
tetrahedral. 

The order of the precision maintained in the 
calculations was such that the frequencies cal- 
culated from the force constants were correct to 
within +2 cm~. In view of the approximate 
nature of the potential function and the tedious- 
ness of the calculations, it was felt that greater 
precision was not warranted. 


21. O. Brockway, J. H. Secrist, and C. M. Lucht, 
Abstracts, Buffalo Meeting of American Chemical Society, 
September, 1942. 
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DISCUSSION 


A comparison of the observed normal fre- 
quencies with those calculated from the set of 
force constants is summarized in Table I for 
CF, and in Table II for CoF¢. 

In the case of CF,, the test of the accuracy of 
the potential function is found in the results for 
the type F. vibrations. The agreement between 
the observed and calculated normal frequencies, 
although poor, represents a considerable im- 
provement over the simple valence-force field. 
The addition of at least another interaction 
term, e.g., angle-angle, to the potential function 
is necessary. 

The observed frequencies for C2F¢, except for 
the type E, frequencies, are reproduced to the 
order of five percent or better. The relatively 
larger discrepancies noted for the E, frequencies 
stem largely from the assumption of no inter- 
action between the CF; groups. From the nature 
of the equations for the type E,, frequencies, it is 
this particular interaction which would be 
effective in reducing the discrepancies between 
the observed and calculated values. The existence 
of this interaction is supported by heat capacity 
measurements which have assigned a barrier of 
the order of 4350 calories per mole, restricting the 
free rotation of the CF; groups.” 

The value of 5.45105 dynes per centimeter 
for the C—C stretching force constant in CeF¢. 
is significantly greater than that observed for 
saturated hydrocarbon molecules. This is in 
apparent agreement with the observation by 
Brockway” that the presence of fluorine atoms, 
adjacent to a C—C bond, has a stabilizing 
influence on it. Hence, a C—C distance of 1.45A 
is observed as compared to the value of 1/54A 
for many saturated hydrocarbons. However, 
because of approximations introduced in the 
potential function, the interpretation that the 


13 E. L. Pace and J. G. Aston, to be published. 
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TABLE I. Comparison of calculated and observed normal 
frequencies for CF,. 











Obs. Calc. 
Designation Type (cm~!) (cm~) 
V1 Ai 904 904* 
V2 E 437 436* 

V3 F, 1265 1316 

v4 F, 630 516 


Force constants in 105 dynes per centimeter: fa=9.15, 
fa=0.71, faa=1.13. 








* Frequencies used to evaluate fa and fg. 


TABLE II. Comparison of calculated and observed normal 
frequencies for CF ¢. 











Obs. Calc. 
Designation Type (cm~) (cm=) 
V1 A lg 1420 1417 
v2 Ag 809 828 
V3 A lg 349 330 
V5 A Qu 1 1 17 1 107 
ve Av 714 697 
V7 Eu 1251 1406 
Vs Ey 523 488 
V9 Ey 216 238 
Vio Ey, 1237 1230 
Vint Eg 620 631 
Vi2 Eo 380 399 


Force constants in 10° dynes per centimeter: fa=9.15, 
fan0.71. faa =1.13, fag=1.35, fg=0.80, fo=5.45, foa—fng 








high value of the force constant is the result of 
a shortening of the C—C bond is not altogether 
justifiable. A more rigorous treatment is required 
to resolve this point. 

From the preceding results, it may be con- 
cluded that the potential functions assumed for 
CF, and C.F, reproduce the observed normal 
frequencies to a good approximation. Moreover, 
an improvement in the agreement between cal- 
culated and observed values would result from 
the addition to the potential function of terms 
for the interaction between bond-bending mo- 
tions, and between the CF; groups in CoF¢. 
However, a normal coordinate treatment in- 
volving such a potential function would be so 
time consuming as to be impracticable. 
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The energy matrix of the rigid asymmetric rotor, 
evaluated in terms of symmetric rotor wave functions, has 
been examined under conditions where J is increased 
indefinitely. Asymptotically, the energy matrix assumes a 
form that is very much like the matrix that may be ob- 
tained from the characteristic value problem of Mathieu's 
differential equation. The characteristic values of Mathieu’s 
equation serve as the basis for a good approximation to the 
energy values of those asymmetric rotor levels which, for 


a given, large value of J, correspond to a small value of K 
in the limiting symmetric cases. The differences which 
exist between the two matrices are accounted for by a 
perturbation technique which permits an accurate deter- 
mination of the energy values. 

When the characteristic values of Mathieu's equation 
lead to a successful approximation to the energy values of 
the asymmetric rotor, an estimate may be made of the 
asymmetric rotor wave functions. 





INTRODUCTION 


LTHOUGH the equations defining the 
energy levels of the rigid asymmetric rotor 
have been derived by various authors,' the 
numerical solution of the equations is not avail- 
able for all values of the rotational quantum 
numbers. Recently, however, King, Hainer, and 
Cross? have published tables of solutions for 
J<10. There have been three methods for 
extension of the solutions to higher J values. 
Numerical solution of the equations becomes 
increasingly laborious for large J’s. The approxi- 
mations resulting from the usual perturbation 
theory yield rather slowly convergent series. 
Finally, there is the application of the corre- 
spondence principle.* 
The result of the correspondence principle 


*The support given by the Navy Department is 
gratefully acknowledged for some of the computational 
work herein reported. It was carried out under Task 
Order V of Contract N5ori-76, Office of Naval Research, 
by Mrs. Dorothy A. H. Brown and Mr. F. C. Merriam. 

** National Research Council Predoctoral Fellow. 

1E. E. Witmer, Proc. Nat. Acad. Sci. 13, 60 (1927); 
S. C. Wang, Phys. Rev. 34, 243 (1929); H. A. Kramer 
and G. P. Ittman, Zeits. f. Physik 53, 533 (1929); 58, 217 
(1929); 60, 633 (1930); O. Klein, Zeits. f. Physik 58, 730 
(1929); H. B. Casimir, Rotation of a Rigid Body in Quantum 
Mechanics (J. B. Wolter’s, The Hague, 1931). 

2G. W. King, R. M. Hainer, and P. C. Cross, J. Chem. 
Phys. 11, 27 (1943). These tables recently have been 
extended to J=11 and 12, R. M. Hainer, Thesis, Brown 
University, 1947. 

3G. W. King, Phys. Rev. 70, 108 (1946); J. Chem. Phys. 
15, 820 (1947); E.E. Witmer, Monthly Progress Reports of 
the University of Pennsylvania, Thermodynamics Research 
Laboratory, Contract NObs-2477, Navy Department 
Bureau of Ships. 
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arguments is to give an asymptotic expression 
for the energy levels which fails to reveal one of 
the most important properties of the asymmetric 
rotor spectrum ; viz., the removal with increasing 
asymmetry of the twofold degeneracy in K (K #0) 
which exists in the limiting case of the symmetric 
rotor. It would appear that the applicability of 
the correspondence principle result is restricted 
to those asymmetric rotor energy levels which 
are essentially twofold degenerate, i.e., under 
conditions of slight asymmetry and large J and 
K (for the limiting symmetric case). 

In order to avoid this restriction, the energy 
matrix of the rigid asymmetric rotor, evaluated 
in terms of a basis of symmetric rotor wave 
functions, has been examined, and the limiting 
form that it assumes when J is increased indefi- 
nitely has been determined. This limiting form 
differs but slightly from that arising in the 
characteristic value problem involving Mathieu’s 
differential equation. As a consequence of this 
asymptotic similarity, the characteristic values 
of Mathieu’s equation have been utilized for a 
first approximation to the energy levels of the 
asymmetric rotor. It has been found possible to 
modify the Mathieu equation by the addition of 
a suitable perturbation operator that permits a 
better approximation to the energy values. In 
principle, a series of suitable perturbation oper- 
ators may be developed which, together with the 
matrix equations derived below, will permit the 
energy values to be determined to as great a 
degree of accuracy as may be desired. 
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ASYMPTOTIC FORM OF THE ENERGY MATRIX‘ 


The energy levels of any rigid asymmetric 
rotor have the form 


a+c a-c 
E,’ (a, b, c) “re (1) 


where J is the quantum number associated with 
the squared total angular momentum of the 
rotor, f?J(J+1); J=0, 1, 2, ---; 7 is a pseudo 
quantum number which indicates the order of 
the 2/+1 energy levels having a given value of 
J; 7 takes on integral values between —J for 
the lowest energy level, and +J for the highest ; 
a=h?/2Iq, b=h?/2I,, c=h?/2I,, and Ia<In< Ii, 
the principal moments of inertia ; x = [(2b—a—c)/ 
(a—c) ], an asymmetry parameter ;> E,/(x) is a 
characteristic value of the energy matrix E(x) 
discussed below. 

The energy matrix E(x), when evaluated in 
terms of symmetric rotor wave functions 
v(J, K, M) as basis functions, does not depend 
upon M in the absence of an external field. 
Furthermore, it is diagonal in J, each J; J block 
being of order 2/+1. In any J; J block it has 
the following non-vanishing elements: 


Ex; x= FJ(J+1)+(G—F)R?, (2) 





RIGID ASYMMETRIC ROTOR 





79 


F, G, and H depend upon the manner in which 
a set of Cartesian axes within the rotor is 
identified with the principal axes,* and upon the 
asymmetry parameter kx. 


f(J, n) =f(J, —n) =3[J(J+1) —n(n+1)] 
X[J(J+1)—n(n—1)]. (4) 
The energy matrix may be written as 


E(x) = FJ(J+1)I+ (G—F)E’(«), (5) 


where I is the unit matrix of appropriate order 
(2J+1), and E’(x) now has the following non- 
vanishing elements in each J; J block: 
B'x, e=K°, (6) 
and 
H 
E'x; x42 =E' +2; Kk=——Lf(J, K+1)]}}, (7) 
G-—F 
provided G— F#0. 
Now for finite , the following asymptotic 
relationship (J) may be obtained: 
rh, n)} eh o( n? (8) 
S|) steers * es pageseg 
4 2 2J(J+1) 





Thus, if the symmetric rotor wave functions 


are ordered: ¥(J, —J, M), --- ¥(J, K, M), 




















and -++¥(J, J, M), the energy matrix E’(x) assumes 
Ex; x+2=Exy2, x=H[f(J, K+1)}. (3) _ the following limiting form. 
H J(J+1) 
2? (0) 0) 
G-F 2 
H J(J+1) 
0) 1° 0 , 
G-F 2 
H J(J+1) H J(J+1) 
E’(k)~| °°: 0 0 0 vee] (9) 
G-F 2 G-—F 2 
H J(J+1) 
1? 0 
G-F 2 
H J(J+1) 
0 22 
G-F 2 








* The notation used by King, Hainer, and Cross, reference 2, will be adopted here. 
° The use of x as an asymmetry parameter is due to B. S. Ray, Zeits. f. Physik 78, 74 (1932), and results from 
his choice of diagonalizing the angular momentum about the intermediate axis of inertia. 


® Reference 2, Table III. 
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neglecting terms O({?/J(J+1)]) compared to 
unity for very large values of J. Equation (9) is 
similar in structure to the corresponding matrix 
that arises in the characteristic value problem 
that involves Mathieu’s equation. It differs from 
the latter in that the former is of order (27+1), 
while the latter is of infinite order; also, the 
latter gives rise to off-diagonal elements which 
are truly constant, while the off-diagonal ele- 
ments of the former are only asymptotically and 
approximately constant for J>K. 


THE MATHIEU EQUATION MATRIX 


To obtain the matrix arising from the char- 
acteristic value problem of Mathieu’s equation, 
consider the equation itself.’ 


d*y 
—+(a—26 cos2x)y=0, (10) 
dx* 


a and @ are parameters which are independent 
of x. 

This equation has both periodic and non- 
periodic solutions, depending upon the values of 
a and 6. However, for those solutions which are 
periodic in 7 or 27 it is found that a may assume 
only special values, which depend upon @. The 
characteristic value problem, then, is to deter- 
mine those values of a, for each value of @, 
which permit solutions periodic in m or 27. This 
may be accomplished as follows. Let #=0, 
obtaining 

d*yo 
—+aoyo=0. (11) 
dx? 


This equation will have solutions periodic in a 
or 27 if and only if aj=K? so that 


1 
(2m)! 





yo= eX, K=O, +1,+42---, (12) 


which form a complete orthogonal and normal- 
ized set of functions in the interval 0<x<2r. 
The quantity (—26cos2x)y in Eq. (10) may be 


7Other forms of Mathieu’s equation are frequently 
used, viz., Eq. (29) below. The particular form here is 
that used by E. L. Ince, Proc. Roy. Soc. Edin. 52, 355 
(1931). An extensive bibliography relating to Mathieu’s 
equation, with summaries of the contents of the most 
important papers, is given by W. G. Bickley, Mathematical 
Tables and Other Aids to Computation I, No. 11 (1945); 
II, No. 13 (1946). 
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regarded as a perturbation term that is added to 
Eq. (11). This perturbation term, as well as 
those solutions y which are periodic in zw or 27, 
may now be expressed in terms of a linear 
combination of the unperturbed functions yp. 
If the unperturbed functions, yo, are ordered: 

















1 1 1 a: 1 
kee eikz e~HK-1)z e* ree, 
(27)} (27)} (2)! 
1 
ei(K-1)z eiKz... 
(27)! (27)# 


application of the usual perturbation procedure 
leads to the matrix 


0 80 0 
0 20 06 0 
M= 6000 6 , (13) 
0 @ 0 1 0 
0 0 60 2 








\ 


the characteristic values of which determine 
those values of a for which solutions of Eq. (10) 
having the required periodicity are obtained. 
Except for being of infinite order, Eq. (13) is 
clearly identical with Eq. (9), if 


H J(J+1) 
G-F 2 


With the aid of Eq. (14) and tables of the 
characteristic values of Mathieu’s equation, an 
estimate may be made of the energy levels of 
the asymmetric rotor. How this may be done 
will be considered in greater detail in a later 
section. 


= 6. (14) 





THE ELLIPTIC CYLINDER FUNCTIONS 


The elliptic cylinder functions are those solu- 
tions of Mathieu’s equation which are periodic 
in + or 2x. Although their properties are well 
known,’ a summary will be given of those 
properties which pertain to the present treat- 
ment. 


8 See W. G. Bickley, reference 7, for numerous references. 
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If Mathieu’s equation is written as 


da? 

(-—+20 cos2e)y=ay, (15) 
dx? 

or 


Ly=ay, 


then it is readily verified that the operator £ is 
left invariant by the transformations 


E: x—-2r+x, 
Ri: x—-2r—x, 
Re: Xx 7T—X, 
R3: Xr +x. 


(16) 


The complete set of orthonormal functions 
generated by Eq. (15) can, at most, experience 
a change in sign when subjected to each of the 
transformations of Eq. (16). Since the set of 
transformations (16) form a group isomorphic 
with the point group Dy», the functions generated 
by Eq. (15) must belong to each of the four 
irreducible representations of De. In this respect 
the elliptic cylinder functions manifest a close 
analogy to the asymmetric rotor wave functions, 
which also belong to the irreducible representa- 
tions of De. 

The elliptic cylinder functions have been 
defined® as 


C€on (x, 0) = >> Ao,” (8) cos2rx, 
r=0 
S€2n41(X, 6) = es Bop 2"* (8) sin(2r+ 1)x, 
r=0 
Cons1(X, 8) = >> A op 1 2"t) (8) cos(2r+1)x, 
r=0 


on 
S€on4.2(X, 6) = : Bo, 2"*?) (6) sin2rx, 
r=0 





where ” takes on all non-negative integral values. 
The characteristic values corresponding to each 
of these species of functions are, respectively, 
2n(8), Don+1(8), don+1(8), Don+2(8). The limiting 
values for @ equal to zero are 
2n(0) = (2n)?, 
bon+1(0) = (2n+1)’, 
G2n41(0) = (2n+1)?, 
Don+2(0) = (2n+-2)?. 


Moreover, from the tabulated values it is found 


(18) 


* Ince, see reference 7. 
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that for any n 


siting <d2n(4) <bon41(8) <Gen41(8) 
Sdan42(8) Sang) (8)-++, (19) 


so that a distinct ordering of the characteristic 
values is possible that is independent of 0. 

In view of the isomorphism existing between 
the elliptic cylinder functions and the asym- 
metric rotor functions it would seem desirable to 
find a one-one correspondence between them, 
based upon symmetry considerations. This would 
allow the unique identification of each character- 
istic value approximated by solution of the 
Mathieu equation with an energy value of the 
asymmetric rotor. This, unfortunately, is im- 
possible because the transformations in Eq. (16) 
may not be identified in a unique fashicn with 
the rotations forming the elements of the point 
group D» to which the asymmetric rotor wave 
functions belong. Depending upon the evenness 
or oddness of J, a representation may be found 
which permits of the appropriate identification 
of elliptic cylinder functions with the wave 
functions of the rotor. 

It is possible, however, to identify the two 
kinds of functions through their corresponding 
characteristic values. In the notation of reference 
2, an asymmetric rotor energy level is denoted 
by the two numbers J,=JK_1,Ky1 with r=K_, 
—K4:1; K 1 is the |K| corresponding to the 
limiting prolate symmetric rotor and K4, is the 
|K| corresponding to the limiting oblate sym- 
metric rotor. 

Observe that either 


K.1+Kii1=J+1 case (a),or (20) 
a J case (b), 
so that 
Kyi=J+1—K-1 
=J—K_; 


case (a), 
case (b). 
Hence 
t=2K_,—J—1 case (a), 
=2K_4—J case (b). 
For the same value of K_, (corresponding to the 
limiting case of K_, degeneracy) 
2K_1—J>2K_1—J—-1. 
Therefore, levels corresponding to case (a) are 


lower than those corresponding to case (b), for 
the same value of K_1. 
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If J is added to both sides of Eq. (20), 


J+K1t+K41=(2J+1) 
=2/J 


case (a), 
case (b), 
so that 


case (a) corresponds to odd (J+K_1+K41) 
case (b) corresponds to even (J+K_1+Ky41). 


Equation (19) may now be applied, permitting 
the following correspondences 


epee for odd (J+K1+K41), (21) 
aK_10JK_1,K41 for even (J+K1+K,41). 

Insofar as the characteristic values of Mathieu’s 
equation approximate the asymmetric rotor 
energy levels, it should be possible by Eq. (21) 
to arrive at an estimate of the asymmetric rotor 
wave functions for J large and |K]| small (for 
the limiting symmetric case). 

The asymmetric rotor wave functions are 
commonly expressed as a linear combination of 
a basis of the Wang! functions, which consist of 
the sums and differences of those pairs of sym- 
metric rotor wave functions which have the 
same J and M but have K’s of opposite sign. 
All these basis functions are normalized to unity, 
and the Wang function for which K is zero is 
identical with the symmetric rotor wave function 
for which K is zero. 

This combination of the symmetric rotor wave 
functions corresponds precisely to the combina- 
tion of the exponential functions generated by 
Eq. (11) to form either sines or cosines. It is in 
terms of these latter functions that the elliptic 
cylinder functions are commonly expressed. As 
a consequence, the coefficients employed in 
expressing the asymmetric rotor wave functions 
in terms of the Wang functions are simply 
related to the Fourier coefficients of the corre- 
sponding elliptic cylinder function. 

In order to employ the tabulated Fourier 
coefficients of Ince, a slight modification is 
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necessary. This modification arises because of 
Ince’s choice of the basis function for which K is 
zero. When allowances are made for this, it 
turns out that all that is needed is to multiply 
the quantities A» by the quantity v2. With this 
modification the transformation coefficients from 
Wang functions to asymmetric rotor wave func- 
tions may be taken directly from the tabulated 
values, or interpolated when necessary. Some 
care must be exercised to maintain the proper 
sign of @ (i.e., as it appears in the energy matrix) 
if the correct characteristic vectors are to be 
obtained. 

Since the characteristic vectors obtained for 
the elliptic cylinder functions are of infinite 
dimensionality, while those for the asymmetric 
rotor are finite (for finite J), it is necessary to 
adjust for this disparity. It is suggested that the 
first m Fourier coefficients of the appropriate 
elliptic cylinder function be taken as the approxi- 
mation to the corresponding asymmetric rotor 
characteristic vector of dimensionality m. Con- 
sistency may be attained by normalizing the 
resulting vector to unity. 


MODIFICATION OF MATHIEU’S EQUATION 


To obtain a better approximation for the 
energy levels than that given simply by the 
characteristic value of Mathieu’s equation the 
higher order terms of Eq. (8) must be considered. 
The following expansion may be obtained 

J(J+1) 
Lf(J, —_ 





n? 1 n? 2 
-—(1+ ~)+0(——_) . (22) 
2 2J(J+1) J(J+1) 
For those cases where nx<J>1, the first two 
terms give an excellent approximation for the 


function. With this approximation, then, Eq. 
(7) becomes 


H jAeTh (K+1)? 1 
E'x; x42=E' x42; k= ( ) - +) 
G-F/t 2 2 2J(J+1) 
H \fJ(J+1) 1 1 
(rst) 
G-F/t 2 2 2J(+1) 








lf H 


[1+ |x++m0. (23) 
2\G-F 2J(J+1) 
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TABLE I. First-order corrections apart from factor 6’. 








\char. 





\ value 
oo beo bor bei boe be2 bes bos bes bos bes bos 
0 0.00000 —1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 J 0.00000 0.00000 0.00000 
1 0.09424 —0.32005 1.83468 1.32876 1.25074 1.19579 1.04153 1.07173 1.05792 1.04121 1.04073 1.02826 
2 0.41579 0.25006 2.84494 2.63144 2.32442 2.50407 1.91671 2.17244 2.06401 2.08152 2.07401 2.05472 
3 0.81976 0.74573 3.97885 3.88779 3.42380 3.88749 2.69191 3.32478 2.97109 3.12456 3.08712 3.07807 
4 1.22537 1.18935 5.14014 5.08654 4.66321 5.31776 3.48317 4.38195 3.74826 4.17802 4.06222 4.09782 
5 1.61321 1.59454 6.26512 6.22294 6.06460 6.77265 4.37809 5.82805 4.39930 5.25192 4.97741 5.11481 
6 1.98021 1.97009 7.33437 7.29965 7.59044 8.23415 5.42184 7.18040 4.97322 6.35700 5.80951 6.13092 
7 2.32744 2.32176 8.34847 8.31365 9.17178 9.68775 6.63650 8.59208 5.54855 7.50303 6.53867 7.15008 
8 2.65694 2.65367 9.31301 9.29165 10.74443 11.12230 8.02937 10.05370 6.20048 8.69766 7.15560 8.17707 
9 2.97075 2.96881 10.23381 10.21786 12.26998 12.52989 9.59010 11.55498 6.98025 9.81414 7.66933 9.21794 
10 3.27069 3.26953 11.11594 11.10419 13.73425 13.93145 11.28789 13.08523 7.91630 10.08346 8.11090 10.27919 








The first bracketed term is constant for a given 
J, and may be identified with the @ in Mathieu’s 
equation. It can be verified that the second term 
in brackets may be obtained from the perturba- 
tion operator 


d? d 
P(x) =— 20 costs ——2 sind] (24) 
dx? dx 


so that Mathieu’s equation, modified by the 
addition of this perturbation operator, is altered 
to 


2 


y 
—-+(a—26 cos2x)y 
dx? 


9 


2y _. dy 
—26'{ cos2x—-—2 sin2x— }=0, (25) 
dx? dx 


with 


=< wae... 
(26) 
(5 Ale aaa) 


Equation (25) will be considered in the next 
section, where it will be solved by a perturbation 
procedure. In principle, additional perturbation 
operators may be found that will permit, with 
one exception, the construction of a secular 
equation that is identical with that arising in 
the asymmetric rotor problem. The exception is 
the disparity in the magnitudes of the order of 
the two secular equations. However, this may 
be handled by a perturbation procedure also. 








PERTURBATION TECHNIQUE 


To take into account the fact that for any 
finite J, the energy matrix of the asymmetric 
rotor is finite, while the matrix generated by 
Eq. (25) (in terms of the functions of Eq. (11)) 
is of infinite order, augment the former by rows 
and columns of zeros, so that E’(x) remains 
symmetrical with respect to the main and anti- 
diagonals, i.e., equal numbers of rows and 
columns added to each side of E’(x). Then, 
symbolically, the augmented matrix may be 
represented as 


E” (x) =M+P-+R, (27) 


where M is the matrix in Eq. (13), P is the 
matrix arising from the perturbation operator in 
Eq. (24), and R is a remainder added to the 
sum of M and P to complete the equality; all 
matrix elements are evaluated in terms of the 
functions generated by Eq. (11). 

The remainder matrix R will have elements 
that are small in magnitude and vanish asymp- 
totically, connecting those diagonal terms which 
are of interest (|K|<J). Hence it will have but 
slight effect upon the characteristic values as 
determined from M+P alone when J is only 
moderately large. 

Formally, Eq. (27) may be diagonalized as 
follows. By its construction M will give char- 
acteristic roots of Mathieu’s equation. These 
may be taken as known. P may be evaluated in 
terms of the periodic solutions of Mathieu’s 
equation, the elliptic cylinder functions, and the 
perturbation operator in Eq. (24). Since a per- 
turbation operator corresponding to R has not 
been obtained it is convenient to first determine 
R in terms of the basis of functions of Eq. (12). 
Transformation to a basis of elliptic cylinder 
functions is then carried out by matrix multi- 
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TABLE II. 
Corresponding 
Asymmetric (J+K-1 Notation Notation _ elliptic cylinder 
rotor level +K41) of ref. 7 of ref. 12 function 
JK.1,K41 odd bx, box+, S€K +, (x, 4) 
JKu1,K4: even aK4, bex+, ceK +, (x, 8) 
TABLE III. 
\\Representation 
Quantity \, I III 
F xk—1 k+1 
2 2 
G—F 3—K ” 3+k 
2 2 
H 1+« 1—« 
G-F 3—k 3+k 








plication. The required transformation matrices 
are readily obtained from the Fourier coefficients 
tabulated by Ince.!° 

In view of their construction both P and R 
may be expected to have elements which have 
absolute values considerably smaller than the 
differences in the characteristic values of the 
levels of M which they connect. Consequently, 
a rapid convergence of the usual perturbation 
series expansion may be expected. 

Some comment is required concerning the sign 
of 6 and @’ in Eq. (26). Since the characteristic 
values of Mathieu’s equation are usually tabu- 





C€on(X, 0): 


P(2n, 2n') = —40 3) (r?—1) {A 22 A op 92") $A 9p OMA op 2") J 5 


r=2 


S€on41(X, 4): 


GOLDEN 


lated for positive 6, it would be convenient to 
have @ always positive. However, since the sign 
of 6 (also 6’) is determined by [H/(G— F)], it 
may be either positive or negative,® depending 
upon the representation that is chosen in which 
to express the symmetric rotor wave functions. 
It is clearly possible to find an orthogonal 
transformation that will change the sign of the 
non-diagonal elements of E(x), while the sign of 
the diagonal elements remains unchanged. Since 
the characteristic values are unaltered by orthog- 
onal transformations, the quantity [H/(G— F) ] 
may be taken positive always. The same does 
not apply to the quantity (G— F) however. 

In evaluating the matrix elements of P in 
in terms of a basis of elliptic cylinder functions, 
a certain simplification is possible based upon 
the symmetry properties of these functions. 
Since the perturbation operator ® has character 
+ 1 for all operations in Eq. (16), it is clear 
that P factors into four submatrices correspond- 
ing to the four species of elliptic cylinder func- 
tions. A similar factoring may be applied to M, 
but is unnecessary, since the characteristic values 
of M may be taken as known. This factoring is 
analogous to the factoring of the secular equation 
of the asymmetric rotor. 

The necessary integrals have been evaluated 
by the use of Eq. (17). The general matrix element 
for each of the four species of functions follows: 


(28a) 


P(2n+1, 2n’+1) = —o| BysrrnByew 4S (4r2—1) 


r=1 


x { Bory 1 2"t Bo, 2” 4 Boy Ort Bo, nt) | (28b) 


C€on41(X, 6) ° 


P(2n+1, 2n’ +1) = —o| — Anan 4B (4r2—1) 


r=1 


X {A ong O2tVA op 20D 4 A 004 Om’tDA 9, 1 Ont) | (28c) 


S€2on42(X, 6) ° 


P(2n+2, 2n' +2) = —40 ¥ (r?—r) { Bo, 2"+2) Bo, 920" +2) + Bo, (2042) Bo, (2mt2) } 


r=2 


(28d) 


10 Ince, see reference 7; J. A. Stratton, P. M. Morse, L. J. Chu, and R. A. Hunter, Elliptic Cylinder and Spheroidal 


Wave Functions (John Wiley and Sons, Inc., New York, 1941); S. Goldstein, Camb. Phil. Soc. Trans. 23, 303 (1927). 
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The diagonal elements of P, apart from the 
factor 6’, have been determined from the coeffi- 
cients given by Ince" for a number of the elliptic 
cylinder functions. These quantities are given 
in Table I. 


COMPUTATIONAL PROCEDURE 


The computation of the characteristic values 
of Mathieu’s equation, which form the basis for 
the approximation considered here, is made 
particularly simple by the recently published 
Tables of Characteristic Values of Mathieu’s 
Differential Equation.” 

Certain differences occur between the notation 
employed in the latter work and that employed 
here. These differences will be summarized. 
Mathieu’s differential equation may be written 
in the form 


d*y 
—+(b—s cos*x)y=0. (29) 
dx? 
Comparison with Eq. (10) reveals that 
titi t 
ie aS, (30) 
et 


The tables give the characteristic values } to 
eight decimal places for various values of s. 
Interpolation for intermediate values of s is 
made particularly simple by the simultaneous 
tabulation of modified second central differences. 
When Everett’s interpolation formula” is used 
with these modified differences the error is less 
than one half-unit in the last tabulated place. 


. 
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The notation used for the characteristic values 
differs from that used above. The following 
correspondences are obtained. 


Aon >beon 

Dong 19 D02n41 
Aon41 be on41 
bon42d02n+2. 


(31) 


The characteristic value of Mathieu’s equation 
corresponding to a particular asymmetric rotor 
energy level may be obtained from Table II. 

Before carrying out a computation for an 
energy value approximation it is best to choose a 
representation that will make | H/(G— F)| small. 
In general, that representation should be chosen 
for which the asymmetric rotor energy matrix is 
most nearly diagonal. Thus for asymmetries in 
the vicinity of the limiting prolate symmetric 
case (x= —1) a type I representation of reference 
2 should be used; in the case of nearly oblate 
symmetry (x=+1) a type III representation 
should be used. Table III gives the quantities 
needed in these computations for both repre- 
sentations. It should be remembered that K_, 
is used with a type I representation, while K,, is 
used with a type III representation. 

The quantities @ and 6’ may now be determined 
with the aid of Eq. (26) and Table III. The 
appropriate characteristic value to be computed 
may be determined from Table II. With the aid 
of Eq. (30) the quantity s may be obtained; the 
value of the appropriate characteristic value } 
may then be obtained from the Tables by 
interpolation, if necessary. The characteristic 


TABLE IV. Approximate second-order correction apart from factor (6’)?. 











beo bor be: bo2 bez bos bes bos 
0 — 0.00000 — 1.12500 — 1.12500 — 5.33333 — 5.33333 — 12.93750 — 12.93750 — 23.46667 
1 —0.39178 — 1.20178 — 1.37094 — 5.41561 — 5.00017 — 12.88162 — 12.71756 — 23.40270 
2 — 0.97928 — 1.41298 — 2.11591 — 5.64637 — 4.58435 — 12.72354 — 12.06048 — 23.22154 
3 — 1.40973 — 1.65907 — 3.26382 — 6.00745 — 4.42034 — 12.55062 — 11.07556 — 22.95080 
4 — 1.72170 — 1.90163 —4,48031 — 6.44982 —4,42417 — 12.39361 — 10.10610 — 22.62205 
5 — 1.96944 — 2.12919 — 5.53756 — 6.93442 — 4.46974 — 12.26193 —9.38104 —22.27251 
6 — 2.18254 — 2.34026 — 6.39651 — 7.44045 — 4.46752 — 12.16077 — 8.92747 —21.91700 
7 — 2.37602 — 2.53642 — 7.09759 — 7.94698 — 4.40606 — 12.09629 — 8.67357 —21.57118 
8 — 2.55769 — 2.71991 — 7.68968 — 8.44328 — 4.34503 — 12.07602 — 8.51668 —21.24234 
9 — 2.73176 — 2.89281 — 8.21038 — 8.92404 — 4.36265 — 12.10735 — 8.35227 —20.93472 
10 — 2.77796 — 3.05685 — 8.67038 —9.38738 —4.51675 — 12.19631 — 8.11453 — 20.65218 











_' Different authors use somewhat different conventions relating to the “normalization” of the functions. Ince 
normalizes” the function to w instead of unity. This fact was taken into account in preparing Table I. 
” A report prepared for the Applied Mathematics Panel, NDRC, by the Mathematical Tables Project, National 


Bureau of Standards, AMP Report 165.1R. 


®E. T. Whittaker and G. Robinson, The Calculus of Observations (Blackie and Son, Ltd., London, 1932). 
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TABLE V. Computation of characteristic values of 
asymmetric rotor. 








Third-order 
perturbation 
formula, ref. 2 


Mathieu 
function 
approximation 


k=—0.9 
— 106.11116 
— 106.67708 
— 100.32955 
—96.96996 
—95.36418 
— 86.74966 
— 86.59730 
—73.14556 


«= —0.8 
— 104.00146 
— 103.90131 
—94.02079 
—92.48039 
— 87.56856 
— 81.32363 
—80.20768 
—67.99363 


Tables of 


Level ref, 2 





— 106.40266 
— 105.67431 
— 100.32513 
—96.97205 
—95.06940 
— 86.75226 
— 86.60144 
—73.14326 


—106.11121 
— 105.67708 
— 100.32958 
— 96.96993 
—95.36420 
— 86.74965 
— 86.59725 
— 73.14548 


100,10 
10;,10 
10:,9 
102, 9 
102.8 
103 
103 7 
1047 


— 106.79625 
— 103.83947 
—94.04603 
—92.99250 
— 84.71850 
— 80.28408 
—80.17775 
—67.95775 


— 104.00141 
— 103.90124 
—94.02065 
—92.48017 
— 87.56942 
— 81.32404 
— 80.20810 
— 67.99406 


100, 10 
10:1, 10 
10:,9 
102, 9 
102s 
103s 
1037 
1047 








value a is then determined by Eq. (30). The 
approximate characteristic value of Eq. (5) is 
then given by 


E,1(«)=FI(J+1)+(G—Fya. (32) 


Corrections to Eq. (32) may be obtained as 
indicated above. The diagonal elements of P, 
except for the factor 6’, are given in Table I, 
from which the appropriate first-order correction 
to Eq. (32) may be determined by interpolation. 
Higher order corrections may be determined by 
evaluation of the non-diagonal elements of P, 
using Eq. (28a)—(28d), followed by the usual 
perturbation treatment. However, more non- 
diagonal elements of P are required that can be 


determined from the few elliptic cylinder func- 
tions that have been tabulated. Nevertheless a 
good approximate second-order correction may 
be obtained from the finite submatrix of M and 
P which contains those levels for which coupling 
elements may be evaluated. The second-order 
corrections so obtained will be larger than the 
complete second-order correction. Because of 
the rapidly increasing difference between succes- 
sive characteristic values of M, the neglect of 
the influence of higher levels will probably not 
prove to be serious. 

The approximate second-order correction to 
be applied to a, except for the factor (6’)?, may 
be determined from Table IV by interpolation, 
if necessary. 

As an illustration of the procedure outlined 
above, the eight lowest energy values for J = 10, 
with asymmetries of k= —0.9, —0.8, have been 
computed. In Table V they are compared with 
the accurate values obtained by King, Hainer, 
and Cross, as well as with the values obtained 
from the third-order perturbation formulas given 
by these authors. It is seen that the agreement 
between the results of the present computation 
and the accurate values is very good, while the 
results of the third-order perturbation calcula- 
tions are quite inadequate for the values of « 
considered. 
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